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Studying  the  Physics  and  Operation  of  Multi-terminal  Near- 
and  Sub-micron  Length,  Hot  Electron  Semiconductor  Devices 


SUMMARY 


This  final  report,  N00014-78-C-0269,  summarizes  recent  studies  by 
H.  L.  Grubin,  the  Principle  Investigator,  as  they  relate  to  near-  and  sub-micron 
device  physics  studies. 
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I .  INTRODUCTION 


Daring  the  last  two  or  three  years  we  have  seen  a  remarkable  resurgence  of 
interest  in  the  physics  of  nonlinear  semiconductor  devices  and  transport.  While 
there  are  probably  as  many  reasons  for  this  development  as  there  are  practitioners, 
we  can  group  them  into  several  categories. 

First  there  are  the  anticipated  fruits  of  the  VLSI  and  VHSIC  programs.  These 
programs  are  geared  primarily  to  the  development  of  conventionally  conceived  devices 
fabricated  on  a  submicron  scale,  and  provide  the  primary  motivation  for  all  other 
submicron  programs.  While  the  VLSI  and  VHSIC  development  will  provide  new, 
important  military  and  commercial  gains,  its  potential  accomplishments  are  limited 
by  our  present  understanding  of  semiconductor  device  physics.  Most  devices  fabri¬ 
cated  as  part  of  this  program  use  design  tools  generated  by  semiclassical  trans¬ 
port,  the  best  known  being  the  mobility  concept  where  the  carriers  respond  instan¬ 
taneously  to  changes  in  electric  field. 

The  inadequacy  of  the  mobility  concept  for  submicron  length  high  speed  devices 
was  raised  implicitly  by  Butcher  and  Hearn  in  1968  XRef .  1) j  and  Rees  in  1969  (Ref.  2); 
and  explicitly  by  Ruch  in  1972  _(Re^-  3^  in  his  well-known  paper  on  velocity  over¬ 
shoot.  Indeed,  Ruch's  paper  was  a  watershed  for  device  transport  studies  insofar 
as  it  focused  attention  on  the  transient  dynamics  of  carrier  transport  and  showed 
that  the  response  of  the  carrier  velocity  to  a  sudden  change  in  field  was  considerably 
different  from  that  associated  with  mobility  models.  i^In  fact,  values  of  velocity 
three  or  four  times  higher  than  steady  state  values  were  commonly  reported.  With 
the  exception  of  a  few  workers,  e.g..  Refs.  4-6,  there  was  not  much  follow-up  to 
Ruch's  study.  The  situation  changed  dramatically  in  the  late  1970's  with  a  plethora 
of  papers  on  transient  carrier  dynamics.  At  this  time,  we  also  saw  private  industry 
increasingly  able  to  fabricate  submicron  length  electron  devices  while  simul taneouslv 
expressing  concern  about  transient  effects  on  device  performance,  particularly  with 
regard  to  GaAs  FET's.  Of  significance  were  such  questions  as:  (1)  could  the  high 
transient  velocities  be  adequately  exploited  in  submicron  design  geometries,  (2) 
could  negative  differential  mobility  be  eliminated  from  sufficiently  short  sub¬ 
micron  devices.  Answers  were  rapid  in  coming,  but  they  were  of  a  patchwork  quality 
insofar  as  they  were  based  primarily  on  uniform  field  calculations.  Thus,  while 
temporal  transients  received  a  great  deal  of  attention,  spatial  transients  were 
largely  ignored.  In  spite  of  the  broad  range  of  assumptions  here,  it  has  become 
increasing  clear  that  a  physics  of  submicron  devices  has  been  developing,  one  that 
we  have  been  calling  the  numerical  physics  of  semiconductor  devices:  see  Grubin, 
et  al.  (Ref.  7). 

The  numerical  physics  of  semiconductor  devices  serves  to  determine  the  tem¬ 
poral  and  spatial  transient  dependence  of  submicron  length  electron  devices  and  to 
account  as  well  as  possible  for  the  environment  surrounding  the  device.  As  indicated 


above,  however,  most  activities  are  based  upon  uniform  field  calculations,  although 
Rees,  and  Rees,  et  al.  (Ref.  8)  have  made  major  advances  in  the  art  of  numerical 
device  analysis  when  nonuniform  fields  are  present.  While  the  need  for  nonuniform 
field  calculations  is  apparent,  Rees'  methods  are  not  finding  widespread  use,  due 
primarily  to  high  computational  costs.  Instead,  for  nonuniform  fields  and  multi¬ 
dimensional  device  simulations,  more  frequent  use  is  made  of  the  method  of  moments 
as  applied  to  the  Boltzmann  transport  equation  (BTE) .  Here  the  most  common  form 
for  the  moment  equation  is  that  summarized  by  Blotekjar  (Ref.  9)  under  the  assump¬ 
tion  that  the  carriers  are  described  by  a  displaced  Maxwellian  function.  Represen¬ 
tative  calculations  using  these  equations  appear,  for  uniform  fields  in  Ref,  1 
(Butcher  and  Hearn),  and  Ref.  10  (Grubin,  Ferry  and  Barker),  for  nonuniform  fields 
and  micron  length  variations  in  Ref.  11  (Bosch  and  Thim) ,  and  Ref.  12  (Cheung  and 
Hearn),  for  submicron  length  devices  in  Ref.  7  (Grubin,  Ferry,  Iafrate  and  Barker). 
And  while  there  is  widespread  concurrence  that  realistic  device  simulations  require 
the  incorporation  of  spatial  derivatives,  the  primary  origins  of  the  field  non¬ 
uniformities  associated  with  these  derivatives  is  less  certain.  These  must  surely 
arise  from  contact  effects  (which  are  imperfectly  understood  from  a  device  view¬ 
point),  material  variations  (can  background  doping  variations  be  adequately  treated 
by  'jelllum  distributions'),  environmental  interactions  (nearby  devices  and  cir¬ 
cuits).  In  addition,  one-dimensional  approximations  to  three-dimension  devices 
are  almost  certainly  inadequate  for  submicron  length  electron  device  simulations. 

A  variety  of  predictions  have  been  made  from  simulations  of  submicron  electron 
device  behavior.  In  particular,  as  the  device  length  shrinks,  the  peak  velocity  is 
expected  to  increase  and  the  regions  of  negative  differential  mobility  becomes  less 
pronounced.  In  addition,  a  ballistic  regime,  where  the  carriers  suffer  no,  or  a  few, 
collisions,  is  thought  by  some  to  be  conceptually  possible  (see  Ref.  13).  The 
experimental  observation  of  such  effects  occupies  the  third  area  of  active  intersest, 
with  workers  in  France,  Laval,  et  al.  (Ref.  14),  and  in  the  U.S.,  involved,  e.g., 
Eastman,  et  al.  (Ref.  15)  and  Zuleg  (Ref.  16).  While  some  of  these  experiments 
clearly  demonstrate  the  length  dependence  of  carrier  transport,  there  are  generally 
too  many  variables  involved  to  allow  an  unequivocable  explanation  of  device 
behavior. 

The  discussion  above,  with  regard  to  the  numerical  physics  of  submicron 
length  electron  devices,  has  focused  primarily  on  problems  associated  with  the 
Boltzmann  transport  equation.  The  BTE  is  the  cornerstone  of  semiclassical  device 
transport  drift  and  diffusion  semiconductor  equations.  Its  central  concept  is  the 
idea  of  a  single  carrier  distribution  function,  f(r,p,t),  which  may  be  used  to 
compute  expectation  values  for  macroscopic  current  flow.  When  we  examine  this 
idea  from  a  quantum  mechanical  viewpoint,  we  see  some  necessary  averaging  if  f  is 
to  be  regarded  as  a  simultaneous  function  of  position  and  momentum.  From  a  device 
perspective,  care  must  therefore  be  exercised  if  f  and/or  the  critical  device 
dimension,  L,  is  smaller  than  an  electron  wavelength,  A.  For  a  quasi-particle 
with  an  effective  mass  m* 
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A  =  h/m  v 


therm 


(1) 


where  A  is  the  thermal  de  Broglie  wavelength.  For  GaAs  central  valley  electrons 
A  -  270  angstroms. 

In  addition  to  spatial  consideration,  from  a  temporal  perspective,  the  BTE 
describes  irreversible  phenomena  and  it  is  assumed  that  collisions  occur  on  a  scale 
short  compared  to  an  observation  time.  Thus,  e.g.,  relaxation  times  should  be 
short  compared  to  a  transit-time.  This  assumption  will  be  routinely  violated  in 
submicron  devices  and  again  caution  must  be  exercised  in  using  the  BTE. 

It  is  clear,  that  from  a  device  modeling  viewpoint,  we  are  faced  with  serious 
problems.  On  the  one  hand,  the  BTE  overcomes  many  objections  to  the  use  of  the 
semi-conductor  drift  and  diffusion  equation.  On  the  other  hand,  the  above  argu¬ 
ments  suggest  abandoning  the  BTE  for  ultrasmall  devices  and  replacing  it  with  a 
quantum  transport  formulation.  While  the  latter  may  be  necessary,  the  current 
approach  is  to  retain  the  BTE  as  long  as  useful  quantum  effects  can  be  incorporated. 
It  is  abandoned  when  first  order  quantum  effects  occur,  as  e.g.,  transport  in 
quantized  inversion  layers. 

A  summary  of  the  above  discussion  is  represented  by  the  flow  diagram  of 
Fig.  1.  The  arrows  flow  from  the  major  motivating  discipline  to  the  new  area  of 
study. 

In  the  following,  we  summarize  recent  work  by  the  Principle  Investigator  on 
ONR  Contract  N00014-78-C-0269  as  it  relates  to  near-  and  submicron  device  physics 
studies.  Detailed  discussions  of  our  studies  are  contained  in  published  or  sub¬ 
mitted  papers.  A  listing  of  these  papers,  along  with  select  copies  of  them  is 
included . 
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II.  SUMMARY  OF  STUDIES  UNDER  ONR  CONTRACT  N000U-78-C-0269 


II. 1  Introduction 

The  studies  carried  out  under  ONR  support  involve  the  formulation  and  simulation 
of  the  transient  behavior  of  transferred  electron  multiterminal  one-  and  two-dimen¬ 
sional  electron  devices.  The  devices  studied  fall  into  two  broad  categories:  (1) 

The  first  is  for  devices  for  which  the  current  response  is  represented  by  the 
equation 


j  =  Ne  y  F+eD  grad  N  +  e 


II 

St 


(2) 


where  F  is  a  field,  y  the  mobility,  and  D  the  diffusivity.  This  equation  is  based 
on  the  Boltzmann  transport  equation  (BTE)  for  a  distribution  function  only  slightly 
modified  by  a  self-consistent  field.  For  high-field-nonlinear-hot-carrier-effects 
the  distribution  function  is  strongly  dependent  on  F  and  the  usual  approach  replaces 
yF  by  a  nonlinear  v(F)  curve,  and  D  by  a  field  dependent  diffusivity;  (2)  The 
second  category  is  for  devices  for  which  Eq.  (2)  is  inadequate  and  the  Boltzmann 
transport  equation  is  needed: 


3f 

St 


where  f(r,p,t) 
kernal. 


-*  3f  -*  S 

+  V  •  —  +  e  F 

Sr  3p 


f(p')  S(p',p)  -  f(p)  S(p.p') 


(3) 


is  the  carrier  distribution  function,  and  W(p',p)  is  a  scattering 


In  our  studies,  Eq.  (2)  is  often  solved  directly;  Eq.  (3)  is  not.  Instead, 
working  under  the  assumption  that  when  the  carrier  concentration  exceeds  a  critical 
concentration  (i.e.,  provided  that  the  rate  of  inter-carrier  energy  exchange  is 
sufficiently  large),  the  carrier  distribution  is  taken  to  be  a  displaced  Maxwellian. 


Then  the  transport  problem  for  a  given  spatially  dependent  field  reduce  to  com¬ 
puting  the  carrier  density,  the  displaced  momentum  and  the  electron  temperature, 
for  each  valley  of  a  multivalley  semiconductor,  such  as  GaAs  or  InP.  These  are 
obtained  by  self-consistently  solving  3A  balance  equations,  where  A  is  the  number 
of  nonequivalent  valleys.  The  balance  equations  are  generated  by  multiplying  the 
BTE  by  successive  powers  of  k;  i.e.,  ^(k)  *  k-' ,  where  j  =  0,1,2  and  integrating 
over  all  k  space.  The  result  is 
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9  i  3  i-*  3  i 

—  <<J>J>  +  ■  <<fc  V>  +  eF  •  < — x  = 

3t  9r  3p  v 


(p')  -  <t>  (P)  S(p.p')  > 


(5) 


where  the  brackets  represent  k-space  integrations: 


<(  )  > 


(  )f  dk 


(6) 


The  first  three  moment  equations  are  well  known  (see  Ref.  9),  and  are  summarized 
below  (Particle  Balance), 
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(Energy  Balance) 
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where  the  subscripts  identify  a  particular  valley.  Here, 


pi  '  "iNivi 


wi  ■  I  Nikii  +  4  nivi2 


(10) 

(ii) 


and  m^  is  the  effective  mass  of  electrons  in  the  iC^  valley.  is  a  measure  of 

the  average  thermal  energy  of  an  electron,  and  is  regarded  as  the  average  total 
kinetic  energy  density.  The  terms  in  brackets  on  the  right-hand  side  of  Eqs.  (7) 
through  (9)  represent  scattering  integrals. 

Equations  (8)  and  (9),  for  momentum  and  energy  balance,  represent  the  essential 
structural  differences  between  the  semiconductor,  Eq.  (4),  and  transport  as  given 
by  the  BTE.  Apart  from  intervalley  transfer,  Eq.  (7),  the  particle  balance  (or 
continuity)  equation,  is  also  required  for  the  drift  and  diffusion  equations. 
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Solutions  to  Eq.  (4),  with  continuity,  and  to  Eqs.  (7)  through  (9),  have  occupied 
most  of  the  Principle  Investigators  efforts  over  the  past  few  years.  We  will 
illustrate  results  of  this  study  below,  but  it  is  instructive  to  see  how  Eq .  (4) 
is  recovered  from  the  balance  equations.  To  do  this,  we  need  explicit  representa¬ 
tions  for  the  scattering  integrals. 


II. 2  Scattering  Integrals 


The  scattering  terms  in  Eqs.  (7)  through  (9)  may  be  approximated  analytically 
or  given  exactly  for  a  displaced  Maxwellian  distribution  function.  In  our  studies 
we  have  chosen  the  latter,  and  use  the  integrals  summarized  in  Butcher's  review 
article  (Ref.  18).  The  scattering  integrals  are  represented 


(12) 


(13) 


(14) 


We  illustrate  a  set  of  scattering  curves  for  a  T-X  orientation  in  Fig.  2.  The 
parameters  used  in  the  calculation  are  given  in  Table  1.  These  results  should  be 
compared  to  those  of  Bosch  and  Thim  (Ref.  11).  We  retain  the  T-X  orientation 
because  most  of  the  early  space  charge  work  involving  the  moment  equations  was 
performed  for  this  ordering. 


Recovery  of  the  drift  and  diffusion  equation  requires  the  judicious  neglect 

of  select  space  and  time  derivatives.  For  the  particle  balance  equation,  conduction 

within  a  single  parabolic  band,  and  neglect  of  avalanching  allows  us  to  set  (3N/3t)  = 

c 

We  then  obtain  current  continuity 


3N  -*• 

—  +  V*(v  N)  =  0 
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(13) 
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For  the  momentum  term,  we  set  (3P4/9t)  =  -P. /ip.  Then,  neglecting  all  spatial 

and  temporal  transients,  i.e.,  the  terms  3P^/3t,  and  V-(vP),  we  obtain 


*>  P  -*  u 

Nv=”=NuF-ii  V(NkT) 
m  e 


where  u  =  ei^/m.  Using  Einstein's  relation,  D  =  pkT/e,  we  obtain 


(16) 


Nv  =  -N  pF  -  V  (ND) 


(17) 


which,  when  coupled  to  the  equation  of  continuity  (Eq.  (15)),  yields  the  semi¬ 
conductor  equation  (Eq.  (2)).  Based  on  the  above  approximations,  it  is  clear  that 
within  the  framework  of  the  Boltzmann  equation,  we  see  that  even  in  the  simplest 
case  the  semiconductor  equations  must  be  used  with  caution. 

The  question  as  to  just  how  serious  the  above  approximations  are,  can  be 

examined  analytically.  However,  for  highly  nonlinear  problems,  a  more  useful 

approach  is  a  numerical  one.  In  perhaps  the  only  serious  comparison,  Cheung  and 

Hearn  (Ref.  12),  studied  Gunn  domain  propagation  using  the  semiconductor  equations 
with  that  obtained  using  the  balance  equations  (neglecting  the  temporal  derivative 
in  the  momentum  balance  equation).  Their  results  show  dramatic  differences  and  are 
displayed  in  Fig.  3. 

While  one  may  have  many  serious  reservations  about  unbridled  use  of  the  semi¬ 
conductor  equations,  its  effectiveness  for  micron  length  and  moderate  frequency 
devices  has  been  established  by  experience.  It  also  has  the  significant  advantage 
of  familiarity.  A  balanced  numerical  physics  program  should  involve  continued 
use  of  the  semiconductor  equations,  as  well  as  the  Boltzmann  transport  equation. 


II. 3  Studies  Using  the  Drift  and  Diffusion  Equation 

In  the  past  several  years  we  have  used  the  semiconductor  equations  to  describe 
time-dependent  transport  in  two  and  three  terminal  negative  differential  mobility 
devices.  The  two  terminal  work  was  concerned  primarily  with  ellucidating  the  opera¬ 
tional  physics  of  GaAs  and  InP  X-band  devices.  Most  of  this  work  has  been  summarized 
in  a  monograph  "The  Gunn-Hilsum  Effect"  (Ref.  III-l,  see  Section  III*)  and  in  a  review 
article  (Ref.  III-2) .  The  studies  described  in  these  reviews  represent  the  non¬ 
linear  semiconductor  as  a  one- dimensional  entity. 


References  preceeded  by  the  Roman  numeral  III  have  been  wholly  or  partially 
supported  by  ONR  Contract  N00014- 78-C-0269 .  They  are  listed  in  Section  III. 


Studies  in  1980,  in  collaboration  with  the  Solid  State  Device  Technology 
Branch  at  the  Naval  Research  Laboratory,  required  a  broadening  of  this  program. 

Here,  using  a  well  developed  two-dimensional  simulation,  we  sought  design  concepts 
to  obtain  repeatable  coherent  microwave  transit  time  oscillations  from  coplanar 
InP  two  terminal  devices  (see  Fig.  4).  The  basis  of  this  study  was  the  observa¬ 
tion  at  NRL  that  repeatable  and  predictable  oscillatory  behavior  from  device  to 
device  had  not  been  forthcoming.  Our  simulations  revealed  that  some  sort  of 
significant  drop  in  the  background  doping  near  the  source  contact  may  be  required 
for  repeatable  and  predictable  oscillations  from  device  to  device. 

The  two-dimensional  simulations  used  in  the  NRL  study  were  developed  earlier 
to  examine  the  time-dependent  behavior  of  GaAs  FET's.  Here  our  initial  theoretical 
work  predicted  the  presence  of  oscillations  in  these  devices  (Ref.  19).  Subsequent 
work,  in  collaboration  with  workers  at  NRL  and  at  Colorado  State  University,  con¬ 
firmed  these  predictions  (Refs.  III-3  and  III-4).  These  results  were  presented  at 
numerous  meetings  and  at  the  NATO  Advanced  Study  Institute  on  the  Physics  of  Non¬ 
linear  Transport  in  Semiconductors.  The  lecture  notes  have  been  published  (Ref. 
III-5)  and  are  part  of  a  major  review  article  to  appear  in  1981  (Ref.  III-6).  In 
addition,  we  have  used  the  two-dimensional  simulation  to  compare  the  switching 
properties  of  GaAs  and  silicon  MESFETs.  We  find  much  higher  switching  speeds  for 
GaAs  at  low  bias  levels,  while  at  high  bias  levels,  due  to  similar  values  of  the 
saturated  drift  velocity,  the  switching  speeds  are  comparable  (Ref.  III-7)  .  There 
appears  to  be  experimental  confirmation  of  these  results  (Ref.  20). 

We  began  an  assessment  of  the  role  of  the  substrate.  Part  of  this  aspect  of 
the  study  was  closely  related  to  the  NRL  collaboration.  In  one  case,  a  two-terminal 
device  whose  channel  height  was  sufficient  to  generate  a  high  field  recycling  domain 
at  a  notched  cathode  was  redesigned  with  a  substrate.  The  substrate  now  occupied 
one-half  of  the  active  region  and  the  oscillation  ceased.  This  result  is  thought 
to  relate  to  a  critical  'thickness-doping  product'  criteria.  In  another  case,  a 
three  terminal  device  with  a  channel  height  sufficient  to  sustain  recycling  high 
field  domains  was  placed  on  a  p-type  substrate  with  varying  doping  levels.  The 
bias  level  sufficient  to  generate  an  instability  in  the  absence  of  a  substrate 
was  insufficient  when  the  substrate  was  present. 

II. 4  Studies  Using  the  Boltzmann  Transport  Equation 

II. 4. i  Transient  Overshoot 

The  second  aspect  of  our  program  involves  solving  the  Boltzmann  transport 
equation  to  examine  the  short  temporal  and  spatial  transients  of  submicron  device 
structures.  At  the  beginning  of  this  program,  the  primary  focus  was  on  evaluating 
the  scattering  integrals.  In  this  evaluation,  we  avoided  most  of  the  approxima¬ 
tions  made  by  others  and  instead  used  the  complete  set  of  integrals  listed  by 


Butcher  (Ref.  18).  The  values  of  the  scattering  integrals  for  the  set  of  parameters 
displayed  in  Table  1  has  been  briefly  mentioned  earlier  and  was  shown  in  Fig.  2. 

The  next  phase  of  the  calculation  involved  computing  the  transient  response 
of  the  carriers  to  a  sudden  change  in  field.  This  is  the  standard  velocity  over¬ 
shoot  calculation  and  was  reviewed  by  the  Principle  Investigator  at  the  NATO 
Advanced  Study  Institute  (summarized  in  Ref.  II1-8).  There  are  several  features 
of  overshoot.  As  discussed  in  an  earlier  paper  (Ref.  III-9),  overshoot  effects 
are  consequences  of  differences  in  momentum  and  energy  relaxation  times.  In  multi¬ 
valley  semiconductors,  the  overshoot  contributions,  therefore,  appear  in  the  momentum, 
as  well  as  the  mean  velocity  computations.  We  review  this  below  for  both  the  cen¬ 
tral  and  satellite  valleys  and  for  the  situations  where  electrons  starting  with  zero 
drift  velocity  are  subjected  to  a  sudden  change  in  electric  field.  In  the  following, 
the  subscripts  "C"  and  "S"  denote  central  and  satellite  valleys. 

For  the  central  valley  and  at  low  values  of  bias,  the  electron  temperature  is 
approximately  equal  to  room  temperature  and  the  ordinary  time  dependent  dynamic 
behavior  occurs.  At  elevated  bias  levels,  the  electron  temperature  is  substantially 
increased  and  the  momentum  relaxation  time,  due  to  strong  intervalley  coupling, 
decreases  with  increasing  temperature  (Fig.  2).  Thus,  we  see  overshoot,  in  that 
the  final  momentum  is  below  the  peak  momentum  (Fig.  5).  (We  point  out  that  over 
moderate  temperature  increases  LO  phonon  intravalley  and  ionized  impurity  scattering 
do  not  provide  a  momentum  relaxation  time  that  decreases  with  energy.  Intervalley 
phonons  are  required.  Indeed,  for  ionized  impurity  scattering,  the  relaxation  time 
increases  with  energy.)  During  this  same  time  interval,  the  increasing  central 
valley  temperature  (Fig.  5b),  results  in  electron  transfer,  and  the  momentum  density, 
Nc  Pc,  shows  an  even  greater  overshoot  (Fig.  5c). 

The  moment  equations  also  contain  a  contribution  from  differential  repopulation. 
This  contribution  enters  as  an  effective  temporal  scattering  term  in  the  momentum 
equation 


d  log  N 
dt 

For  the  central  valley,  where  at  t  =  0,  Ng  ^  0,  and  Pc  =  Ps  =  0,  this  term  is 
approximately  zero.  However,  when  we  consider  the  transient  behavior  of  the 
satellite  valley,  the  time  derivative  of  logNg  is  important  because  the  change  in 
the  satellite  population,  relative  to  the  original  number  present,  is  quite  large. 
In  addition,  we  note  that  the  satellite  valley  momentum  relaxation  contribution  is 
almost  an  order  of  magnitude  larger  than  that  of  the  central  valley.  Thus,  in  a 
time  considerably  shorter  than  that  associated  with  the  central  valley,  the  satel¬ 
lite  valley  momentum  reaches  the  value 
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where,  since  the  satellite  temperature  remains  close  to  room  temperature  for  large 
changes  in  field,  the  scattering  rates  may  often  be  taken  as  approximately  constant. 
Combining  both  scattering  contributions  we  see  some  "overshoot"  due  to  differential 
repopulation,  but  none  as  dramatic  as  that  associated  with  the  central  valley. 

Figure  6  illustrates.  For  part  a,  1  corresponds  to  the  increasing  momentum  prior 
to  any  significant  electron  transfer;  2  represents  the  scattering  rate  due  to 
repopulation;  3  is  the  steady  state  value.  In  Fig.  6  we  show  the  product  of  Ns 
Ps/ (Nc  +  Ng) ,  where  the  effect  of  electron  transfer  is  included.  We  see  that  as 
far  as  the  contribution  to  velocity  overshoot  is  concerned,  the  electron  transfer 
tends  to  wipe  it  out. 


The  complete  calculation  requires  that  we  obtain  the  drift  velocity.  For  the 
calculations  associated  with  Figs.  5  and  6,  the  drift  velocity  is  shown  in  Fig.  7. 
Perhaps  the  most  remarkable  aspect  of  this  result  is  the  very  large  peak  velocity 
prior  to  steady  state  (see,  e.g..  Ref.  3).  This  result  has  been  one  among  many 
that  have  led  some  to  suppose  that  narrow  channel  devices  will  yield  higher  carrier 
velocities.  To  a  large  extent,  very  high  overshoot  velocities  are  illusory  as  they 
are  very  sensitive  to  rise  time.  We  illustrate  this  for  a  sequence  of  trapazoidal 
bias  pulses,  each  with  a  varying  rise  time  (see  Fig.  8). 

II.4.ii  Bias  Dependence  of  Overshoot 

The  first  set  of  results  is  for  a  relatively  slow  rise  time  and  the  dynamic 
curves  come  very  close  to  the  steady  state  curves  (see  Fig.  9).  A  more  significant 
departure  from  steady  state  occurs  for  the  somewhat  steeper  rise  time.  In  Fig.  10, 
we  see  some  asymmetry  in  the  time  dependence  of  the  central  valley  population  and 
temperature,  and  an  increase  in  the  peak  veloc'Ly.  In  the  final  sequence,  we  show 
results  for  a  very  short  rise  time.  We  see  a  dramatic  increase  in  the  peak  velocity 
and  clear  asymmetry  in  the  carrier  dynamics  (Fig.  11).  Indeed,  the  final  point  of 
approximately  zero-field  and  velocity  is  not  an  equilibrium  state.  Rather,  we 
have  a  dramatic  example  of  velocity  undershoot.  A  longer  time  is  needed  for  the 
electron  temperature  to  approach  equilibrium.  There  are  strong  implications  here 
for  upper  frequency  limits  of  device  operation. 

Perhaps  the  earliest  attempt  to  examine  the  upper  frequency  limit  for  large 
signal  oscillation  was  that  of  Butcher  and  Hearn  (Ref.  1).  Here,  using  a  set  of 
displaced  Maxwellian  electron  distributions  for  each  valley,  the  set  of  differential 
equations  for  the  time  dependent  electron  temperatures,  drift  velocities,  and  valley 
populations  was  solved  for  a  dc  bias  field  plus  rf  field. 


II-7 


II.4.iii  Self-Excited  Oscillations 

From  the  point  of  view  of  a  device  physics,  a  driven  oscillator  probably  lies 
somewhere  between  a  small  signal  and  a  large  signal  self-excited  oscillator.  The 
self-excited  oscillator  is  perhaps  the  most  interesting  of  the  three  because  it 
highlights  the  tenacious  balance  between  electron  transfer  and  sustained  oscillators. 
It  is  extremely  sensitive  to  contact,  space  charge  and  circuit  conditions.  We 
have  been  examining  the  upper  freuqency  limit  for  self-excited  oscillations  for 
the  device  in  a  circuit  with  reactive  elements.  Figure  12  shows  the  oscillation. 

Here  the  current  through  the  load  resistor  is  displayed  in  part  b;  the  dynamic 
voltage  and  I  versus  <£,  obtained  by  eliminating  time  between  current  and  device 
is  shown  in  Fig.  12a.  We  also  display  the  mean  velocity  (dynamic  conduction  current). 
The  details  of  the  oscillation  are  as  follows: 

As  the  field  across  the  device  increases  and  exceeds  a  threshold  value  transfer 
begins  to  occur.  But  since  the  field  changes  more  rapidly  than  the  electron  tem¬ 
perature,  more  carriers  are  retained  in  the  central  valley,  and  with  higher  momenta, 
than  steady  state  ^ould  dictate.  This  effect  is  responsible  for  the  higher  peak 
conduction  current.  But  if  the  increasing  electric  field  sustains  high  fields 
for  a  sufficient  duration  enough  carriers  will  transfer  to  result  in  negative 
differential  mobility  (NDM) .  NDM  of  sufficient  magnitude  is  required  for  sustained 
self-excited  oscillations. 

On  the  downswing  the  field  again  changes  more  rapidly  than  the  electron  tem¬ 
perature  and  more  carriers  are  retained  in  the  satellite  valley  -  we  achieve 
transient  undershoot.  Now  while  NDM  is  not  necessary  on  the  downswing  (Ref.  III-l) 
enough  carriers  must  be  returned  to  the  high  mobility  valley  for  transfer  on  the 
upswing  and  NDM  to  occur.  This  means  that  the  field  must  change  slowly  enough  to 
allow  a  dump  of  carriers  from  the  satellite  to  the  central  valley.  If  the  field 
changes  too  rapidly  too  many  carriers  are  retained  in  the  satellite  valley  and 
NDM  is  too  weak  to  sustain  steady  state  oscillations.  In  Fig.  13,  we  plot  the 
maximum  frequency  of  self-excited  oscillations  as  a  function  of  dc  bias. 

We  note  the  results  in  Figs.  5  through  13  are  taken  from  Ref.  III-6. 


11.4. iv  Length  Dependence  of  Negative  Differential  Mobility  Region 


We  have  been  discussing  the  upper  frequency  limit  of  transferred  electron 
devices  from  the  circuit  viewpoint  and  the  transfer  and  return  of  electrons  between 
the  central  and  subsidiary  valleys.  In  this  analysis  we  have  not  explicitly  con¬ 
sidered  the  following  problem:  When  an  electron  enters  the  active  region  of  a 
device,  it  will  accelerate  in  the  presence  of  an  applied  electric  field.  If  the 
initial  drift  velocity  of  the  carrier  is  low,  is  the  transit  length  sufficient  to 
cause  electron  transfer  and  negative  differential  mobility?  The  answer  lies  in 
earlier  calculations.  If  the  carrier  experiences  a  sudden  change  in  field,  the 
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Initial  transient  (t  <<  LO  phonon  scattering  time)  will  be  ballistic,  followed  by 
a  region  where  v  will  approach  -eiF/m*  for  a  single  valley.  If  the  transit  time 
is  short  enough  to  prevent  significant  transfer  NDM  will  be  weak,  if  at  all. 

Figure  14  summarizes  (Ref.  III-10)  where,  for  uniform  fields,  device  length  is  a 
derived  quantity.  Here, 

T 

L  =  J  v ( t )  dt  (19) 

o 

where  T  is  a  transit  time. 

Now  the  velocity  versus  field  curves  and  the  velocity  versus  time  curves  of 
the  type  shown  here  provide  an  indication  of  why  there  is  interest  in  submicron 
devices.  The  possibility  exists  for  achieving  very  high  velocities  over  very 
short  distances.  But  again,  a  word  of  caution.  The  calculations  of  Fig.  14  are 
for  carriers  subjected  to  sudden  changes  in  field.  As  we  have  seen,  finite  rise 
time  dramatically  reduce  this  peak,  and  so  the  results  will  be  somewhat  less 
important . 

II. 4. v  Spatial  Transients  in  Short  Channel  GaAs  Devices 

Velocity  overshoot  and  other  uniform  field  transient  calculations  are  a  very 
common  feature  of  submicron  device  physics  studies.  Spatial  transients,  of  the 
type  discussed  by  Kroemer  (Ref.  21),  where  rapid  changes  in  field,  as  in  a  high 
field  Gunn  domain,  can  cause  velocity  overshoot  transients  has  received  less  ex¬ 
posure.  Spatial  transient  calculations  are,  however,  likely  to  be  essential  for 
realistically  determining  the  effectiveness  of  overshoot  in  improving  the  frequency 
response  of,  e.g.,  submicron  field  effect  transistors.  We  have  begun  to  examine 
some  of  these  problems  by  solving  the  full  spatially  dependent  one-dimensional 
Boltzmann  transport  equations  (Eqs.  (7)  through  (9)).  We  do  this  for  a  one-dimen¬ 
sional  5000  A  GaAs  element  with  a  donor  distribution  as  shown  in  Fig.  15.  The 
element  is  part  of  a  resistive  circuit  which  is  biased  at  a  finite  rate.  The 
carrier  dynamics  is  examined  at  two  instants  of  time.  (Note:  serious  objections 
can  be  raised  to  use  of  "jellium"  distributions.  For  any  combination  of  increased 
donor  density  or  size  reduction  will  necessarily  introduce  effects  due  to  the 
discrete  nature  of  the  donors.  This  is  ignored  below.) 

As  the  bias  is  turned  on,  there  is  an  increase  in  potential  across  the  device, 
and  a  corresponding  increase  in  current  and  field.  The  field  is  computed  self- 
consistentlv  and  its  slopes  reflect  any  incomplete  screening  of  inhomogeneities  by 
the  mobile  carriers.  For  the  device  in  the  schematic  configurations  of  Fig.  15, 
as  the  field  in  rising  energy  relaxation  is  incomplete  and  velocity  overshoot  con¬ 
tributions  are  dramatic  (see  Figs.  16  and  17).  Here,  velocity  is  computed  from  the 
equat ion 
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In  steady  state  the  extent  of  the  spatial  transient  will  depend  on  the  bias.  In 
Figs.  18  and  19,  the  fields  are  high  enough  to  accommodate  almost  complete  trans¬ 
fer.  Here,  the  carrier  temperature  is  reduced,  the  energy  relaxation  rates  are 
shorter,  and  there  is  virtually  no  overshoot.  The  results  here  are  virtually 
the  same  as  we  would  obtain  using  the  steady  state  curves.  The  nonuniform  field 
calculations  are  summarized  in  Ref.  III-6. 
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I.  INTRODUCTION 

The  operating  principles  of  the  field  effect 
transistor  were  set  down  by  Shockley  in  1951. 
The  device  was  a  three-t erminal  device  with  one  of 
its  terminals  connected  to  a  pn  junction  (see  Fig. 
la).  The  pn  junction  provided  a  region  of  local 
charge  depletion  whose  width  W(x)  could  be  varied 
by  changing  the  bia&  on  the  gate  contact.  In  this 
way  the  resistance  to  the  flow  of  carriers  from 
the  source  to  the  drain  contact  could  be  altered 
and  current-voltage  relations  of  the  type  ahovn  in 
Fig.  lb  resulted. 

Since  1952,  as  new  problems  arose  and  a  variety 
of  different  materials  were  used,  FET  behavior 
has  required  models  increasingly  more  sophisti¬ 
cated  than  that  proposed  by  Shockley.  In  general, 
however,  as  new  models  were  proposed;  their 
validity  was  soon  called  into  question.  For 
example,  in  Shockleys  analysis,  which  was  often 
applied  to  the  material  silicon,  the  mobility  was 
assumed  to  be  constant.  While  this  was  expected 
to  be  reasonably  accurate  for  low  bias  levels  it 
was  also  applied  to  those  situations  where  the 
presence  of  high  fields  made  the  assumption 

inadequate.  Modifications  have  since  been  made 

2 

which  include  velocity  saturation  . 

Along  with  velocity  saturation  comes  the  notion 
that  carriers  in  the  FET  behave  like  a  compressi¬ 
ble  fluid,  bunching  together  and  yielding  local 
regions  of  high  electric  field**.  When  this 
occurs  in  other  semiconductors  such  as  gallium 
arsendie,  dramatic  effects  may  occur.  In  the  case 
of  gallium  arsenide,  while  early  interest  was 
based  on  its  high  mobility;  it  was  recognised, 
almost  at  the  onset,  that  nonlinearities  in  the 
carrier  velocity  would  have  to  be  included.  At 
the  beginning,  and  for  FETa,  the  carrier  velocity 
was  modelled  with  two  sections,  with  the  high- 
field  section  accounting  for  velocity  satura¬ 
tion^.  While  this  could  predict  some  of  the 
CaAs  FET  properties,  it  could  not  explain  the 
current  instabilities  that  were  reported  in 


1 9  66  ^  and  1968^.  In  1976-79,  the  first  papers 
to  deal  with  these  dynamic  effects  began  to 
appear  .  These  calculations  were  store  general 
than  the  earlier  studies  in  that  the  nonlineari¬ 
ties  in  the  velocity  electric  field  relation 
included  a  region  of  negative  differential  mo- 
bility,  as  shown  in  Fig.  2a.  Currently,  the  uae 
of  these  curves,  which  represents  steady  state 
equilibrium  situations  for  explaining  the  dynamic 
behavior  of  electrons,  is  being  called  into 
question'®.  In  particular,  ita  uae  for  explain¬ 
ing  highly  doped  submicron  device  behavior  is 
being  treated  gingerly.  The  questions  of  where 
and  when  the  equilibrium  velocity  electric  field 
curve  of  Fig.  2a  can  be  used  is  as  yet  unanswered. 
Meanwhile,  we  will  discuss  what  it  predicts.  We 
will  also  sutmaarize  our  experimental  results  with 

CaAs  Schottky  gate  FETS  which  tend  to  confirm 
11 

these  predict  ions  . 

II.  INSTABILITIES  IN  SCHOTTKY  CATE  GALLIUM 
ARSENIDE  FIELD  EFFECT  TtytNSISTORS 

II. a  Summary 

The  static  and  dynamic  properties  of  a  Schottky 
gate  gallium  arsenide  field  effect  transistor  are 
summarized  in  Fig.  2b.  The  bold  curve  in  the 
figure  is  the  equilibrium  velocity  electric  field 
curve  scaled  to  the  current  and  voltage  parame¬ 
ters.  The  current  is  in  multiples  of  N_ev  A, 
ldiere  NQ  is  the  active  region  doping  level,  v^  the 
peek  carrier  velocity,  and  A  the  cross-section 
area  in  the  direction  of  current  flow.  The 
voltage  ie  in  multiples  of  E  L,  where  is  the 
field  at  peak  velocity.  Vitn  reference  to  Fig. 
2b,  we  have  found  that  whenever  the  drain  current 
density  I  £  Ig,  the  value  associated  with  the 
saturated  drift  velocity  of  the  carriers,  an 
instability  can  occur.  In  two  terminal  device 
language,  I(  is  referred  to  as  the  sustaining 
current*2.  Generally,  beceute  of  the  nonlinear 
distribution  of  potential  under  the  gate  contact 
very  thin  FETa  sustain  lower  current  deoaity 
levels  than  thicker  devices.  This  particular 
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feature  has  led  to  the  folic  wing  device  classifi¬ 
cation  scheme**,  We  have  classified  devices  as 
being  Berbers  of  one  of  two  groups  as  determined 
by  the  rat  10 

^  m  _ _ Cate  voltage  at  cutoff _ 

Drain  voltage  at  the  onset  of  current 
saturation  for  zero  gate  voltage 

Devices  with  K  greater  than  unity  sustain  current 
oscillations,  those  with  K  approximately  equal  or 
less  than  unity  are  electrically  stable.  For 
those  cases  with  K  greater  than  unity,  oscilla¬ 
tions  begin  beyond  the  knee  of  current  saturation 
(denoted  by  dashed  line  in  Fig.  2b).  Oscillations 
cease  at  high  values  of  drain  and  gate  bias  (the 
solid  lines  in  Fig.  2b).  The  instability  is 
initiated  when  the  average  field  under  the  gate 
contact  reaches  the  NDM  threshold  value.  (This 
result  holds  true  at  different  temperatures.)  An 
estimate  of  the  carrier  velocity  in  the  channel 
between  the  gate  and  drain  contacts  indicates  that 
just  prior  to  the  oscillation  the  electrons  are 
usually  drifting  at  speeds  in  excess  of  the  high 
field  saturation  drift  velocity.  When  the  oscil¬ 
lation  ceases  with  increased  gate  bias,  the 
velocity  drops  to  values  significantly  below  the 
saturation  drift  velocity.  The  above  conclusions 
emerge  from  numerical  simulation  and  experiment. 

II. b  Numerical  Results 

The  operation  of  the  CaAs  FET  is  critically 
dependent  on  the  behavior  of  local  regions  of  high 
electric  field  of  'domain*.  And  depending  on 
device  size  and  bias  values,  these  domains  are 
either  trapped  near  the  gate  or  drain  contacts  or 
propagate  between  them**.  The  numerical  simula¬ 
tion  reveals  the  presence  of  two  qualitatively 
different  trapped  and  one  traveling  domain.  The 
first  trapped  domain  forms  within  the  conducting 
channel  and  at  the  drain  edge  of  the  gate  contact. 
Its  formation  is  a  direct  consequence  of  velocity 
saturation,  and  its  presence  is  responsible  for 
current  saturation  at  values  of  drain  voltage 
significantly  below  that  predicted  by  Shockley* 
(thus  values  of  K  greater  than  unity).  The  second 
trapped  domain  is  qua  1 i tat ive 1 y  different  and 
forms  in  FETS  capable  of  sustaining  current 
inst ab i 1 i t les  arising  from  negative  differential 
mobility.  This  domain  forms  near  the  drain 
contact  at  high  drain  and  low  gate  bias  levels 
and  is  coincident  with  the  cessation  of  oscil¬ 
latory  activity  in  GaAs  FETs.  The  domain  is 
qualitatively  similar  to  the  high  anode  field*^ 
conf igurat ion  in  two  terminal  devices,  where  it 
is  responsible  in  many  caaea  for  ceaaation  of 
oscillations.  In  the  case  of  the  traveling 
domain,  it  is  launched  at  the  drain  side  of  the 
gate  contact,  propagatea  to  the  drain  contact,  and 
then  recycles.  The  general  properties  of  thia 


domain  are  qualitatively  similar  to  that  for  two 
terminal  devices;  and  when  the  current  density 

is  somewhat  below  1  ,  all  oscillatory  activity 
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ceases  . 

II.  c  Experimental  Results 

The  oscillatory  activity  of  the  CaAs  FET6  was 
detected  with  spectrum  analyzers  and  sampling 
oscilloscopes.  Spectral  analysis  on  one  device 
showed  oscillations  at  6.5,  13,  23,  30,  36,  and  40 
GHs.  Higher  frequencies  could  not  be  detected. 
The  oscillations  were  strongest  near  the  insta¬ 
bility  threshold  but  ceased  at  high  drain  bias 
levels.  At  high  drain  and  zero  gate  bias  levels, 
white  light  was  detected.  An  initial  reduction  in 
bias  on  the  gate  contact  resulted  in  similar 
behavior,  but  all  oscillatory  activity  ceased  at 
sufficiently  high  gate  bias  levels,  including  the 
appearance  of  light.  The  appearance  of  light  is 
compelling  evidence  for  the  presence  of  regions  of 
high  electric  field  domains.  Others  have  also 
observed  white  light*'*  at  high  drain  bias  levels 
and  attributed  it  to  impact  ionization  within  a 
domain  trapped  at  the  drain  contact.  The  inten¬ 
sity  of  light  was  also  observed  to  increase  with 
drain  bias  and  decrease  with  increased  gate 
bias. 

III.  CONCLUSIONS 

Based  upon  the  results  of  our  numerical  and 
experimental  study  we  have  reason  to  conclude 
that  our  observations  of  oscillations  in  gallium 
arsenide  field  effect  transistors  arise  from 
negative  differential  mobility  contributions. 
Although  the  possibility  exists  chat  these  os¬ 
cillations  are  a  consequence  of  feedback  between 
the  gate  and  drain  loop,  the  critical  parameters 
for  the  initiation  and  cessation  of  oscillations 
appear  to  be  determined  by  the  material  proper¬ 
ties  of  the  device  rather  than  the  circuit.  Of 
course  the  circuit  and  the  device  may  determine 
the  oscillation  frequency  much  in  the  same  way  of 
circuit  controlled  oscillations  occur  in  two 
terminal  device;  but  in  the  study  reported  here 
the  oscillation  conditions  are  primarily  materials 
dependent.  The  details  of  our  result  appear 
elsewhere*  * . 
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la.  Schematic  representation  of  a  junction  field  effect  transistor 


DRIFT  VELOCITY  (107cm/(*d  ”  SOURCE  -  DRAIN  CURRTNT 
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i  source  drain  characteristics  of  a  junction  field  effect  transistor 


ELECTRIC  FIELD  IkV/cm) 


Drift  velocity  versus  electric  field  for  a  semiconductor  with  a  region  of 
negative  differential  mobility 


?b 
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SOURCE  -  DRAIN  VOLTAGE  V/VD 


Expected  current  voltage  characteristics  for  a  negative  differential  mobility 
semiconductor.  Stable  regions  are  denoted  by  solid  curves,  unstable  regions 
by  dashed  curves.  Instabilities  require  Ip>Ig.  Devices  where  Ip(V^*0)  lg 
will  be  electrically  stable. 
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I.  INTRODUCTION 

As  devices  are  made  smaller  and  doping  levels  increased,  transient 
relaxation  effects  begin  to  dominate  the  behavior  of  transferred  electron 
devices.  One  area  that  has  received  most  attention  is  that  of  velocity 
overshoot  in  the  central  valley,  an  effect  due  mainly  to  differences  in 
the  energy  and  momentum  scattering  times  (Ref.  1).  There  are,  however, 
at  least  two  other  relaxation  effects  of  interest.  The  first  is  impor¬ 
tant  mainly  in  the  satellite  valley  and  is  a  direct  consequence  of  its 
short  momentum  scattering  time.  The  second  is  a  consequence  of  the  fact 
that  scattering  between  two  different  quantum  states  takes  place  over  a 
finite  duration  of  time  (Kef.  2).  In  the  discussion  that  follows  we  will 
review  overshoot  contributions  within  the  central  valley  and  then  discuss 
the  remaining  transient  effects.  Our  results  are  presented  within  the 
framework  of  solutions  to  the  moments  of  the  Boltzmann  transport  equa¬ 
tion,  although  many  of  the  results  have  also  been  obtained  using  Monte 
Carlo  methods.  In  these  calculations  we  have  also  included  a  resistive 
circuit.  Its  consequences  will  be  briefly  discussed. 


In  an  actual  device  the  electic  field  will  be  a  function  of  space 
and  time.  The  present  calculations  assumes,  however,  a  spatially  uniform 
electric  field  between  cathode  and  anode  of  a  two  terminal  device.  The 
electrons  are  assumed  to  enter  the  device  with  a  zero  drift  velocity. 
The  response  of  the  average  electron  to  a  step  change  in  bias  is  then 
obtained  by  simultaneously  solving  the  external  circuit  equation  with  the 
zeroth,  first  and  second  moments  of  the  Boltzmann  transport  equation 
assuming  that  the  electron  distribution  may  be  approximated  by  a  shifted 
Maxwe 1 1-Boltzmann  distribution  function  (Ref.  3).  The  zeroth  moment 
yields  the  equation  of  continuity,  and  Che  first  moment  yields  Newton's 
law  of  motion  for  momentum  density: 


i 


-eK 


dP 

dt 
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In  equation  (1)  K  denotes  tut-  momentum  scattering  rate,  its  units  art 
1/sec.  The  seconc  noDent  equation  describes  er.erg;-  balance  and  its  tic.-, 
dependence  is  governea  by  its  scattering'  rates.  Indeea,  tt.c  injjrtor,: 
transient  relaxation  effects  arise  because  the  menu  r.t  ur.  scattering  rate 
is  an  explicit  function  of  electron  tempera  t  ur  e ,  and  because  for  materi¬ 
als  like  GaAs ,  the  time  rate  of  change  of  momer.  turn  and  energy,  which  are 
governea  by  their  scattering  rates,  are  different.  Typical  central 
valley  momentum  and  energy  scattering  rates  for  1  -A  ordering  are  shown  in 
figure  la  (.see  also  Ref.  3) 


11.  CENTRAL  VALLET  RELAXATION 


Central  valley  overshoot  can  be  understood  by  examining  solutions  to 
equation  (1).  For  the  case  of  a  constant,  temperature  independent 
momentum  scattering  time,  the  solution  to  this  equation  subject  to 
P(t=0)=0  is  (ignoring  dlogN/dt): 


P ( t )  =  -eR  1  -  exp  -R  t 

P  P 


(2) 


and  the  time  dependence  is  similar  to  tnat  shown  in  figure  2a,  which  was 
calculated  exactly.  For  the  calculation  of  figure  2a  the  bias  was 
sufficiently  low  so  that  only  a  marginal  increase  in  the  central  valley 
electron  temperature  occurs.  For  sufficiently  high  values  of  bias  the 
electron  temperature  reaches  a  high  enough  value  to  cause  a  sizeable 
increase  in  tne  scattering  rate.  The  result  is  that  shown  in  figure  2b. 
If  we  examine  the  structure  of  the  central  valley  momentum,  we  find  that 
although  the  electron  temperature  is  increasing  the  initial  momentum 
transient  is  governed  primarily  by  the  momentum  balance  equations  (1)  and 
whatever  resistance  is  in  the  external  circuit.  If  the  latter  can  be 

ignored,  then  for  an  average  scattering  rate  of  3.5  x  10  /sec  over  a 

temperature  range  of  300  to  1000°K,  Pc  will  be  within  75%  of  its  peak 
value  -eF/Rp  after  a  time  duration  of  approximately  O.Aps.  The  central 
valley  electron  temperature,  which  has  also  been  increasing  (figure  2c) 
over  this  time  period,  continues  to  increase.  And  its  increase,  which  is 

governed  mainly  by  the  external  circuit  parameters  and  the  energy  relaxa¬ 

tion  time,  governs  the  successive  increase  in  the  momentum  scattering 
rate  and  the  decrease  in  the  steady  state  value  of  the  central  valley 
momentum  (see  figure  2b). 


The  overshoot  displayed  in  figure  2D  is  that  associated  with  a 
single  carrier.  This  is  not  the  quantity  usually  discussed  when  velocity 
overshoot  is  considered.  Rather,  the  average  drift  velocity 

V  =  Vc/mc*  -  NsPfi/ms* 


(3) 


where  Pg  is  the  satellite  valley  moment  urn  ;  mc*,  m,.*,  Nc  and  N£  are 
respectively,  the  central  and  satellite  valley  effective  masses  and 
net  carrier  densities.  In  figure  2d  we  display  the  quantity 


NCPC/  (Nc  *  Ns) 

versus  time.  We  also  superimpose  the  curves  of  figures  2b  and  2c  for 
reference.  We  note  that  the  oversnoot  is  considerably  enhanced  by  elec¬ 
tron  transfer.  Also  note  that  the  relaxation  of  the  central  valley  drift 
velocity  closely  follows  the  relaxation  of  the  electron  temperature. 
This  result  holds  true  even  for  different  resistive  contributions  and 
suggests  that  the  relaxation  to  steady  state  in  an  actual  device/circuit 
configuration  will  be  governed  by  both  energy  relaxation  and  the  external 
c i rcui t . 


III.  SATELLITE  VALLEY  relaxation 
In  the  above  discussion  we  have  ignored  the  quantity 

dlogN 

dt 

appearing  in  equation  1.  For  the  discussion  so  far,  in  particular, 
that  associated  with  figure  1  where  at  t  =  0,  N  =  0,  and  P  *=  P  =  0, 
this  is  approximately  true.  However  when  we  consider  the  transient 
behavior  of  the  eat  ellite  valley,  the  time  derivative  of  logN£  is 
significant . 

To  examine  the  transient  behavior  of  the  satellite  valley  we  first 
estimate  the  extent  to  which  the  satellite  temperature  exceeds  room 
temperature.  Typically,  according  to  the  results  of  reference  A,  the 
satellite  temperature  is  approximately  equal  to  its  room  temperature 
value  for  fields  below  AO  kv/cm.  Over  this  same  temperature  range  the 
central  valley  temperature  can  exceed  200Q°K.  The  origin  of  this  dif¬ 
ference  lies  in  the  scattering  rates  for  the  central  and  satellite 
valleys  which  we  show  below. 

In  steady  state  the  energy  balance  equation  yields  (kef.  3) 

0  =  -NcFPc/mc*  -NckBTcK(5)  +  N£TSK(6)  (4) 

0  =  -n’fPs/*£*  -VbTsR(7)  4  NcTcR(b)  (5) 

for  the  central  and  satellite  valleys  respectively.  lne  prime  on  Ng 
denotes  that  N£  is  the  density  of  only  one  of  the  satellite  valleys. 
The  numbers  on  the  scattering  rates  refer  to  the  central  and  satellite 
valley  rates  of  figure  2.  Using  the  steady  state  values  of  momentum  and 
equations  (4)  and  (5)  we  find  the  following  relation  for  the  central  and 
satellite  temperatures: 


(6) 


NgTs  =  h(3)h(b)  ■*  (ms*/mc*  )k(  -  ;«(  b ) 

NCTC  K(  3 1 kl 6 )  +  (Dis*/mc*)K(4)K(  1) 

wher*  k(3),  k(3)  and  R(fe)  are  central  valley  scattering  rates.  Tne 

others  are  satellite  valley  rates.  ns  an  example  of  the  relative  values 
°-  and  Ts  wc  consider  the  time  rate  of  change  of  carriers  in  the 

central  va 1  ley : 


d_Nf  =  -Nck(1)  +  n’r(2)  (7) 

dt 

From  figure  Id  we  see  that  the  satellite  valley  scattering  rates  are 
relatively  insentive  to  temperature  wnile  the  central  valley  scattering 
rates  are  temperature  sensitive.  We  can  estimate  that  in  steady  state 
Ng  =  would  occur  with  a  central  valley  electron  temperature  of 

approximately  1500°K.  The  corresponding  value  ot  Tg  as  obtained  from 
equation  (b)  is  slightly  larger  than  SOO^k.  The  implication  of  this  is 
that  satellite  valley  momentum  overshoot  is  not  sensitively  dependent  on 
the  satellite  valley  temperature. 

As  mentioned  above,  in  examining  the  transient  behavior  of  the 
central  valley  we  have  ignored  contributions  from  the  time  derivative  of 
log  Nc.  From  equations  (1)  and  (7)  we  see  tnat  the  net  scattering  rate 
i  s 


R(3)  -  K( 1 )  +  R(2)Ng/Nc. 

For  the  initial  transient  where  Ng^^Nc,  and  from  the  results  of  figure 
2  we  see  that  the  dominant  scattering  contribution  is  from  R(3).  The 
situation  is  considerably  different  when  we  consider  the  satellite 
valley.  For  this  case,  with  3  satellite  valleys,  as  for  F-X  ordering  the 
net  momentum  scattering  rate  is 

RU)  +  R(l)N’c/(3Ng)  -  KU)/3 

Now  prior  to  intervalley  transfei  the  sum  of  the  last  two  terms  is  zero. 
But  at  elevated  central  valley  temperature  R(ll  increases  rapidly  with 
electron  temperature.  This  coupled  with  the  fact  that  at  the  initial 
stages  of  electron  transfer  the  ratio  of  Nc/3Ng  is  large  will  lead 
to  the  result  that  the  onset  of  transfer  will  significantly  increase  the 
satellite  valley  scattering  rate. 


In  addition  to  the  importance  of  the  term  dlogNg/dt,  we  also  note 
that  the  satellite  valley  momentum  relaxation  contribution  K(4)  is  almost 
an  order  of  magnitude  large  than  that  of  the  central  valley.  Thus  in  a 
time  considerably  shorter  than  that  associated  with  the  central  valley  we 
expect  the  satellite  valley  momentum  to  reach  its  steady  state  value 
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-eFU) 


K(4)  ♦  aloghs 


-1 


lb) 


dt 


**-*«■»> 


c  r<  Wt-  can  take  K(4)  to  be  constant  and  equal  to  jts  r  oom  tempi  ra:  _rt 
value.  f(t)  will  continue  to  change  in  a  manner  determines  b\  the 
external  circuit  and  the  central  valley  relaxation.  Thus  we  expect  that 
P  will  reach  a  maximum  value  given  approximately  by  -eF(tj/K(4)  and  then 
decrtoSt  as  the  contribution  from  dlogN^/dt  enters.  In  sttad\  state 
tne  final  value  of  P£  will  be  given,  again  by  eF/iUa).  This  result  is 
sn own  in  figure  3a.  In  figure  3a,  (l)  corresponds  to  the  increasing 
momentum  prior  to  any  significant  electron  transfer;  ©  represents  the 
scattering  rate  due  to  repopulation;  ©  is  the  steady  state  value.  lr. 
figure  3b  we  show  tne  product  of  N'P,./(N  +  3N').  We  see  that  as  far 
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as  the  contribution  to  velocity  overshoot  is  concerned  the  electron 
transfer  tends  to  wipe  it  out. 


IV.  FINITE  QUANTUM  TRANSITION  TIMES 


Most  electron  device  calculations  in  which  scattering  events  are 
accounted  for  are  based  upon  the  assumption  that  the-  transitions  between 
two  states  within  a  particular  valley,  or  between  state  within  different 
valleys  takes  place  during  a  zero  time  duration.  This  is  an  approxima¬ 
tion;  the  transition  can  take  place  over  a  finite  time  duration  and  this 
time  duration  may  be  comparable  to  the  transient  response  time.  Ferry 
and  Barker  (Kef.  2)  have  suggested  an  approximation  to  account  for  these 
finite  collision  events.  This  is  to  replace  the  scattering  rates  'R', 
which  are  implicit  functions  of  time  by  the  following  expression 


HU) 


exp(-t'/tc)  R(t  -  t’  )dt 1 


(y> 


Thus  the  instantaneous  scattering  rate  is  averaged  over  an  effective 
collision  duration  t  weighted  by  the  function  exp  (— t/rc ) .  If  as 
in  the  case  of  low  fields  k(t)  does  not  change  much  during  a  collision 
duration,  the  normal  result  is  obtained.  However,  in  a  large  field, 
where  finite  quantum  transition  times  are  important  the  weighting  is 
necessary.  These  effects  will  also  be  important  in  high  frequency 
transport  where  the  collision  duration  becomes  comparable  to  the  relaxa¬ 
tion  times  and  the  period  of  the  wave. 


The  effect  of  including  finite  quantum  transition  times  is  illustrateo  in 
figure  U  where  we  show  the  transient  response  for  a  steady  uniform  rield 
of  5kv/cm  and  25kv/cm.  The  dashed  lines  are  obtained  for  the  case  of 
instantaneous  transition  times.  We  see  that  in  both  cases  the  response 
for  finite  transition  times  is  the  faster  one  (kef.  2).  We  note  that 
results  displayed  in  figure  4  were  obtained  for  a  two  band  F-L  ordering, 
within  the  framework  of  the  shifted  Maxwellian  distribution. 


V.  CONCLUSIONS 


The  purpose  of  the  study  was  to  highlight  several  new  features 
associated  with  the  transient  behavior  of  electrons  in  GaAs.  Particular 
emphasis  was  on  overshoot  in  the  satellite  valley  where  its  origin  was 
due  primarily  to  the  onset  of  electron  transfer.  A  similar  type  of 
overshoot  may  be  expected  to  be  associated  with  the  central  valley  at 
very  high  bias  levels.  In  the  latter  case  the  central  valley  is  almost 
depleted  of  free  carriers  and  a  decrease  in  bias  will  result  in  a  trans¬ 
fer  of  carriers  from  the  satellite  to  central  valley  with  resulting 
repopulation  overshoot  effects.  In  addition  we  also  discussed  the 
consequences  of  finite  quantum  transition  times  on  the  transient  behavior 
of  carriers  in  GaAs.  Noting  that  these  times  may  be  comparable  to  the 
scattering  times  or  the  period  of  a  high  frequency  signal,  the  effect  of 
these  finite  times  on  the  response  of  the  carriers  was  computed.  It  was 
shown  that  the  effect  of  a  finite  transition  time  is  to  increase  the 
scattering  rates  and  hence  reduce  the  response  times  of  the  carriers. 
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Figure  2.  Calculations  for  T-X  ordering,  (a)  Central  valley  momentum 
versus  time  for  a  spatially  uniform  field.  The  average  field  across 
the  device  plus  load  resistor  is  2.2  kV/cm.  T0  is  the  room  temperature 
value,  and  t0  the  central  valley  collision  time  for  zero  drift  energy. 
t0  *  3.17  x  10-33  sec.  (b) ,  (c)  and  (d):  For  an  average  field  across 
device  and  load  resistor  of  10  kV/cm,  the  central  valley  momentum, 
electron  temperature  and  momentum  density,  respectively. 
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INTRODUCTION 

In  the  previous  discussion  (Grubin,  1977)  we  concentrated  on  tran¬ 
sient  hot  electron  effects  in  semiconductor  devices.  We  noted  that 
these  effects  may  be  dominant  for  spatial  dimensions  of  the  order 
of  the  carrier  mean  free  path  and  for  doping  levels  common  to  cur¬ 
rently  fabricated  field  effect  transistors.  Perhaps  one  of  the  most 
significant  consequences  of  these  relaxation  effects  is  that  they 
will  ultimately  serve  to  elimate  the  presence  of  negative  differen¬ 
tial  mobility  from  such  semiconductors  as  gallium  arsenide  and  indium 
phosphide.  Thus,  there  are  at  least  two  questions  that  we  must  con¬ 
sider:  (1)  Under  what  conditions  will  a  semiconductor  device  possess 

a  region  of  negative  differential  mobility,  and  (2)  if  it  does  possess 
this  region,  what  are  its  consequences?  We  will  examine  these  ques¬ 
tions  in  reverse  order. 

In  discussing  the  consequences  of  negative  differential  mobility 
(NDM)  we  are  at  that  point  in  time  where  most  of  our  device  pictures 
are  derived  from  models  in  which  the  equilibrium  velocity  electric 
field  relation  is  taken  to  represent  carrier  drift  in  nonequilibrium 
situations.  While  the  discussion  of  transient  effects  points  to 
severe  limitations  with  this  model  at  extremes  of  size  and  doping 
levels,  dynamic  effects  associated  with  such  things  as  intervally 
scattering  times  will  limit  its  usefulness  to  frequencies  at  and 
below  the  low  millimeter-wave  scale.  Further,  the  equilibrium 
velocity  electric  field  curve  is  dependent  on  impurity  densities  and 
so  may  vary  from  one  region  of  the  device  to  another.  With  these 
qualifications,  it  is  clear  that  equilibrium  velocity  electric  field 
models  cannot  be  used  confidentially  for  providing  quantitative 
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determination  of  device  behvaior.  Rather,  we  may  expect  that  in 
those  devices  whose  dynamic-relaxation  contributions  still  allow  a 
region  of  NDM  to  be  present,  such  equilibrium  velocity  electric  field 
models  will  provide  the  information  necessary  to  present  qualitative 
explanations  of  the  oscillatory  characteristics  of  NDM  semiconductor 
devices.  Indeed  this  is  the  very  least  we  can  expect  of  the  equili¬ 
brium  models;  and  as  we  will  see,  they  have  served  us  well. 

In  discussing  NDM  semiconductor  device  behavior  through  the  equili¬ 
brium  velocity  electric  field  model  we  are  going  to  concentrate  on: 
(1/  the  role  of  the  low  resistance  contacts  on  the  oscillatory  pro¬ 
perties  of  the  NDM  device,  and  (2)  the  influence  of  the  external  cir¬ 
cuit.  We  will  first  examine  the  two-terminal  device,  then  move  to 
the  three-terminal  device. 


TWO-TERMINAL  DEVICES 

To  highlight  the  importance  of  the  cathode  or  anode  region  properties 
on  the  oscillatory  behavior  of  two  terminal  NDM  devices,  we  point  out 
that  while  at  the  early  stages  of  device  development  most  devices 
were  fabricated  with  ostensibly  similar  cathode  and  anode  region  con¬ 
tacts,  it  was  the  usual  case  to  find  that  the  current  voltage  charac¬ 
teristics  of  the  device  were  not  symmetrical  with  respect  to  the 
polarity  of  the  applied  bias  (Shaw  et  aT . ,  1969).  The  importance  of 
the  circuit  emerges  when  we  realize  that  all  current  oscillations  in 
NDM  devices  are  not  characterized  by  isolated  propagating  domains 
whose  transit  time  determines  the  oscillation  frequence  of  the  device 
(Gunn,  1964).  In  some  cases,  the  oscillation  frequence  bears  no 
resemblance  to  the  transit  time  of  the  carriers  (Thim,  1968). 

With  regard  to  the  low-resistance  contact  to  materials  such  as  gallium 
arsenide  and  indium  phosphide,  there  appears  to  be  general  agreement 
that  these  contacts  are  highly  nonlinear  and  have  a  non-negligible 
effect  on  device  behavior  at  all  except  the  lowest  bias  levels. 

There  is  not,  however,  any  theory  that  self-consistently  determines 
the  properties  of  these  contacts  from  first  principles  and  then 
relates  these  properties  to  device  behavior.  Rather,  there  are  a 
variety  of  phenomenological  models.  In  our  case,  we  have  adopted  the 
point  of  view  that  the  alloyed  metal/semiconductor  contact  imitates 
the  properties  of  the  unalloyed  metal/semiconductor  contact.  The 
important  distinction  between  these  being  that  the  alloyed  contact 
has  a  substantially  lower  barrier  height  than  the  unalloyed  contact. 

In  analyzing  the  role  of  the  contacts  to  the  NDM  device,  we  con¬ 
ceptually  divide  the  device  into  three  sections,  as  shown  in  Fig.  1. 

The  center  region  in  Fig.  1  represents  the  NDM  semiconductor.  It  is 
characterized  by  an  equilibrium  velocity  electric  field  relation  v(E) 
and  an  equilibriumdiffusion  electric  field  relation  D(E).  The  time- 
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Fig.  1.  Device  -  circuit  configuration. 


dependent  dynamics  of  the  carriers  are  obtained  by  solving  the  semi¬ 
conductor  equations,  which  for  one  dimension  reduce  to  that  shown 
below,  for  the  total  current  density  through  the  NDM  semiconductor 
(see,  e.g.,  Shaw  et  al . ,  1979), 

yo  ■  v<E)  [eyx)  +  t  - 

82E(x,t)  9N  (x)  9E(x,t) 

-D(E)  e  - r - +  e  — - -  +  e  -  .  (1) 

-  9x  9x  J  9t 

The  above  equation  deals  with  only  one  carrier  specie,  whose  charge 
we  take  conveniently  as  the  magnitude  of  the  electron  charge.  Here, 
e  is  the  permittivity  of  the  semiconductor,  Nq(x)  is  the  background 
doping  density,  which  may  be  spatially  dependent,  and  E(x,t)  is  the 
electric  field. 
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The  carrier  dynamics  in  the  cathode  contact  region  are  taken  to  be 
governed  by  the  differential  equation  (Grubin,  1976a) 

j0(t)  =  W  +  XEEc/dt>  (2> 

where  jc(Ec)  represents  the  cathode  conduction  current  density  and 
has  the  same  general  properties  as  the  'control  characteristic'  dis¬ 
cussed  by  Kroemer  (1968).  Ec  is  the  electric  field  at  the  boundary 
between  the  cathode  region  and  the  semiconductor  region.  Depending 
on  the  contact,  it  may  assume  very  high  values  as  well  as  values  near 
zero.  X  is  a  dimensionless  parameter.  As  indicated  in  Fig.  1,  we 
assume  that  the  cathode  region,  as  represented  by  (2),  separates  the 
NDM  semiconductor  from  the  cathode  proper.  The  rationale  for  this 
separation  lies  in  experiments  which  indicate  that  for  contacts  with 
high  cathode  fields,  the  transition  to  the  high  field  values  may  take 
place  outside  the  NDM  semiconductor.  The  cathode  contact  region  is 
also  taken  to  be  electrically  neutral  and,  for  the  cases  considered 
below,  to  have  a  permittivity  equal  to  that  of  the  NDM  semiconductor 
(i.e.,  A  =  1).  Ec  is  assumed  to  be  continuous  across  the  interface 
between  the  NDM  and  separating  elements. 

As  indicated  above,  the  principle  assumption  we  make  about  the 
alloyed  contact  is  that  it  imitates  the  properties  of  the  unalloyed  metal/ 
semiconductor  contact.  As  used  here,  this  implies  that 

W  =  -jr  t  exp(-eVc/nkT)-exp  [-(n  1  -  l)eVc/kT]  }  ,  (3) 

with  Vc  =  ECLC.  The  form  of  (3)  was  adpated  from  studies  of  the 
unalloyed  metal/semiconductor  contact  (Rideout,  1975).  Its  use  here 
presumes  a  similar  description.  For  the  unalloyed  contact,  'n'  is 
the  'ideality'  factor  and  describes  the  contact  as  dominated  by 
thermionic  emission  (n *  1)  or  by  tunneling  (n>>l);  jr  is  the  reverse 
current  flux  and  may  be  related  to  the  barrier  height  through  the 
Richardson-Dushmann  equation 

jr  =  m*AT2exp[-(B/kT)]  ,  (4) 

where  B  is  the  barrier  height  (in  ev),  A  is  the  Richardson-Dushmann 
constant  (=  120  A/cm2/k2),  and  T  is  the  absolute  temperature,  m*  is 
the  ratio  of  an  appropriate  effective  mass  to  that  of  the  free  elec¬ 
tron  mass.  The  results  depend  significantly  on  the  value  of  jr  which 
we  translate  into  values  for  B.  For  this  we  have  arbitrarily  assigned 
a  value  of  0.072  to  m*  as  representative  of  the  central  valley  of  CaAs. 
The  relation  between  jr  and  B  is  displayed  in  Table  1,  where  jp  = 
N0evp»  with  vp  being  the  peak  drift  velocity  of  the  carriers  as  ob¬ 
tained  from  the  equilibrium  velocity  electric  field  relation.  Typi¬ 
cal  barrier  heights  for  our  discussion  are  about  0.2  ev.  The 


TABLE  1 

BARRIER  HEIGHT 
-B=kT  loge (jr/m*AT2) 


j  /j 

B 

r  p 

1.00 
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parameter  Lc  ,  assoc iated  with  ( 3)  ,  is  specif  ic  to  the  forumulation  we  use 
and  does  not  appear  in  the  formulation  for  metal/semiconductor  contacts. 
While  the  parameter  Lc  is  conceptually  ambiguous,  we  regard  it  as 
representing  the  width  of  the  alloyed  region;  n  and  jr  are  thought  of 
as  more  representative  of  the  properties  of  the  metal/semiconductor 
interface. 


Figure  2  displays  the  dependence  of  the  form  of  jc(Ec)  on  the 
choice  of  parameters.  Note  that  significantly  different  parameters  can 
yield  similar  jc(Ec)  curves.  We  have  also  included  for  reference  the 
neutral  current  density  curve  jc(E)  =  Noev(E)  which  scales  the  velo¬ 
city  electric  field  curve  for  gallium  arsenide.  The  significance  of 
these  curves  for  the  space  charge  distribution  within  the  NDM  semi¬ 
conductor  device  is  considered  below  for  steady-state  time  indepen¬ 
dent  conditions.  For  this  situation,  we  return  to  the  total  current 
density  equation  and  examine  it  for  zero  diffusion  and  for  constant 
doping  density.  At  the  cathode  boundary. 


“.WlW 

9x  ev(E  ) 

c 


(5) 


whereas  within  the  bulk  of  the  semiconductor,  away  from  either  the 
cathode  or  anode  regions 


E  =  j/N^ey  ,  (6) 

U  o 

where  uQ  is  the  low  field  mobility  of  the  semiconductor. 

In  Fig.  3,  we  sketch  the  electric  field  profile  for  x>0  for  the 
situation  where  the  cathode  conduction  current  density  curve  inter¬ 
sects  the  neutral  current  density  curve  at  a  point  somewhere  below 
the  threshold  field  for  negative  differential  mobility.  For  success¬ 
fully  higher  current  density  values,  by  virtue  of  current  continuity, 
the  intersections  of  the  vertical  line  j  =  jj_  with  jc  and  jn  deter¬ 
mines,  respectively,  the  cj  thode  field  Eci  and  the  bulk  field  E^. 

The  value  of  jn  at  Eci  exceeds  jc  at  this  point,  and  there  is  a  par¬ 
tially  depleted  region  of  charge  adjacent  to  the  cathode  region 
boundary.  Increasing  the  current  density  to  the  value  j2»  where  jc 
and  jn  intersect,  yields  a  situation  of  uniform  fields  throughout  the 
cathode  plus  semiconductor  region.  For  the  situation  where  the  current 
density  is  increased  to  j 3 ,  the  field  at  the  cathode  is  less  than  the 
field  in  the  bulk,  jn  is  less  than  jc  and  we  have  a  region  of  charge 
accumulation  adjacent  to  the  cathode  contact  region  boundary.  A  cur¬ 
rent  instability  may  be  expected  to  occur  when  the  field  away  from 
the  cathode  region  enters  the  region  of  negative  differential  mobil¬ 
ity. 
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Fig.  2.  Cathode  conduction  current  density  versus  electric  field 
for  seven  of  parameters,  n,  jr  and  Lc .  Also  shown  is  the 
neutral  conduction  current  density  curve  for  gallium 
arsenide,  denoted  by  Jn(E).  From  Shaw  et  al.  (1979). 


The  next  situation  to  consider  is  that  shown  in  Fig.  4  ,  where  the 
cathode  current  density  curve  intersects  the  neutral  curve  in  the 
region  of  negative  differential  mobility.  At  the  current  density  jj., 
we  have  as  in  the  previous  case,  a  region  of  charge  depletion  occurs, 
adjacent  to  the  cathode  region  boundary.  Increasing  the  current 
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Fig.  3.  Electric  field  versus  distance  profiles  for  the  situation 
where  jc(Ec)  and  Jn(Ec)  intersect  at  values  of  electric 
field  below  the  NDM  threshold  field.  The  shaded  region 
denotes  the  region  of  negative  differential  mobility. 


density  to  J2  results  in  an  increase  in  cathode  field  to  a  value  some¬ 
where  within  the  NDM  region.  The  field  within  the  bulk  is  still 
within  the  low  field  region  as  indicated.  Further  increases  in 
current  density  J3  result  in  the  intersection  of  the  jn  and  jc  curves 
and  charge  neutrality  at  the  cathode  region  boundary.  Note  that  the 
field  within  the  bulk  is  still  significantly  less  than  the  field  at 
the  cathode  boundary  region.  The  situation  at  33  is  unstable  because 
now  further  increases  in  current  density  result  in  the  formation  of 
an  accumulation  layer  adjacent  to  the  cathode  region,  followed  by  a 
downstream  depletion  layer.  This  domain  is  wholly  within  the  NDM 
region  and  may  generate  a  large  signal  instability.  A  point  worth 
drawing  attention  to  at  this  time  is  made  by  comparing  the  current 
densities  and  the  estimates  of  the  average  fields  across  the  semi¬ 
conductor  elements,  represented  by  Figs.  3  and  4.  In  Fig.  3,  the 
current  density  at  the  instability  threshold  is  near  the  peak  value 
for  NDM.  Also,  the  average  field  just  prior  to  the  instability  is 
approximately  equal  the  threshold  field  for  NDM.  For  the  results 
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————— -CURRENT  DENSITY  DISTANCE - - 

Fig.  4.  As  in  Fig.  3,  but  jc(Ec)  and  jn(Ec)  intersect  within  the 
region  of  negative  differential  mobility. 


of  Fig.  4,  both  the  current  densities  and  the  average  fields  prior  to 
the  instability  are  significantly  less  than  the  NDM  threshold  values. 

In  Fig.  5,  we  illustrate  the  situation  where  the  cathode  conduction 
current  density  intersects  the  neutral  curve  somewhere  within  the 
region  of  saturated  drift  velocity  region.  For  this  case,  we  do  not 
expect  any  charge  accumulation  until  the  value  of  current  density 
exceeds  the  current  density  associated  with  the  saturated  drift  velo¬ 
city. 

To  examine  the  dynamic  properties  of  these  space  charge  layers,  we 
resort  to  numerical  simulation  and  solve  the  semiconductor  equations 
simultaneously  with  the  circuit  equations.  For  this  we  need  the  anode 
contact  region  conditions.  This  is  an  easier  contact  to  specify 
phenomenologically,  for  we  have  found  that  if  the  region  is  not  tai¬ 
lored  to  meet  specific  and  unusual  requirements,  its  effect  on  the 
oscillation  and  instability  conditions  can  be  simulated  by  setting 
the  anode  field  value  to  zero.  With  this  in  mind,  we  examine  the 
simulation  results. 
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Fig.  5.  As  in  Fig.  3,-  but  jc(Ec)  and  jn(Ec)  intersect  within  the 
region  of  saturated  drift  velocity. 


We  being  with  the  situation  where  the  cathode  conduction  current 
density  curve  intersects  the  neutral  curve  within  the  region  of  nega¬ 
tive  differential  mobility.  For  a  resistive  circuit,  Z  =  R,we  get  the 
result  in  Fig.  6  (Shaw  jet  £l . ,  1979)  which  shows  the  prethreshold 
electric  field  versus  distance  profile  and  the  postthreshold  electric 
field  versus  distance  profile  (at  one  instant  of  time).  The  nuclea- 
tion  criteria  is  the  same  as  that  discussed  in  Fig.  4,  and  the  re¬ 
sulting  domain  travels  between  the  contacts.  It  drains  at  the  anode 
and  is  renucleated  at  the  cathode  region.  The  dc  current  voltage 
curves  (Grubin,  1976b)  are  shown  in  the  inset.  We  point  out  the 
sublinearity,  just  prior  to  the  instability  (represented  by  the 
’wiggle’  in  the  curve).  The  sublinearity  is  due  to  the  enhanced 
voltage  drop  in  the  vicinity  of  the  cathode  contact.  In  a  circuit 
containing  reactive  components,  e.g.,  serial  inductance  and  parallel 
package  capacitance,  the  dynamic  behavior  of  an  NDM  element  with  a 
similar  contact  may  be  different.  This  is  displayed  in  Fig.  7a  where 
we  plot  current  density  through  a  dc  load  resistor  versus  an  average 
electric  field,  E  (Grubin,  1978).  Time  is  eliminated  between  the  two. 
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Fig.  6.  Resistive  circuit  calculations  for  a  10  micron  long  element 
with  jc(Ec)  and  jn(Ec)  intersecting  within  the  region  of 
NDM.  Curves  1,  2,  and  3  denote  pre-instability  curves  at 
increasing  values  of  bias.  Curve  4  denotes  an  instability 
at  one  instant  of  time.  From  Shaw  et  al.  (1979).  The 
current  density  versus  average  electric  field  curve  for 
this  element  is  shown  in  the  graph  on  the  right.  From 
Grub in  (1976b). 


In  Fig.  7b,  we  show  the  electric  field  versus  distance  profile  at 
four  instants  of  time.  Also  shown  is  the  cathode  conduction  current  den¬ 
sity  curve.  We  note  that  in  the  case  of  Fig.  6,  the  instability  is 
dominated  by  the  cathode  contact  and  the  frequency  is  determined  by 
the  transit  time  of  the  domain.  In  the  circuit  containing  reactive 
components,  the  frequency  here  is  determined  primarily  by  the  circuit 
elements  (although  this  is  not  always  the  case).  We  note  here  that 
the  oscillation  appears  to  consist  of  a  moving  oversized  dipole  layer. 

We  show  a  similar  set  of  calculations  for  a  device  in  which  the 
cathode  conduction  current  density  curve  intersects  the  neutral  curve 
further  within  the  region  of  negative  differential  mobility,  close  to 
the  region  of  velocity  saturation.  The  results  show  a  much  weaker 
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Fig.  7.  Reactive  circuit  computations  for  a  jc  curve  that  inter¬ 
sects  jn  within  the  NDM  region.  From  Grubin  (1976b). 


current  oscillation.  Indeed  with  cathode  curves  of  this  type,  as 
they  intersect  further  out  we  may  expect  the  oscillatory  activity  to 
disappear.  The  current  voltage  curves  for  this  device  are  shown  in 
the  inset  (Grubin,  1976a)  and  we  see  a  situation  where  the  subline¬ 
arity  is  extreme  enough  to  approach  saturation.  This  excessive  sub¬ 
linearity  is  due  to  the  large  voltage  drop  in  the  vicinity  of  the 
cathode  contact,  as  illustrated  in  Fig.  8. 

The  space  charge  prof  lies  discussed  above  can  be  intuitively  under¬ 
stood  on  the  basis  of  some  simple  analytical  concepts  associated  with 
the  intersection  of  the  neutral  current  density  curves  and  the  cathode 
conduction  current  density  curve.  There  is  another  profile  that 
occurs  frequently  in  device  operation  that  is  less  amenable  to  simple 
arguments.  It  occurs  most  often  in  devices  whose  cathode  current 
density  curves  intersect  the  neutral  curves  at  cathode  field  values 
Ec  <  Ep,  and  it  manifests  itself  as  a  large  voltage  drop  in  the  vici¬ 
nity  of  the  anode  contact.  This  domain  requires  a  region  of  NDM  for 
its  existence,  but  does  not  require  the  presence  of  any  special  anode 
property.  The  presence  of  these  anode  domains  was  first  predicted  by 
Shockley  (1954).  In  resistive  circuits  its  presence  usually  occurs 
at  high  bias  levels  and  often  accompanies  the  cessation  of  oscilla¬ 
tions.  In  circuits  containing  reactive  elements  it  usually  manifests 
itself  as  a  permanent  and  residual  voltage  drop  in  the  vicinity  of  the 
anode  as  illustrated  in  Fig.  9.  (In  Fig.  9,  rather  than  deal  with  the 
cathode  conduction  current  curves  of  Fig.  2,  we  have  simply  set  the 
cathode  field  to  zero.)  We  show  the  lissajous  figure  for  the  oscil¬ 
lation  as  well  as  a  multiple  exposure  of  the  electric  field  versus 
distance  profile.  This  profile  represents  an  accumulation  layer  that 
propagates  toward  the  anode  but  never  quenches  at  the  anode.  Instead 
there  is  always  a  residual  space  charge  layer  adjacent  to  the  anode 
contact . 


Each  of  the  above  calculations  was  carried  out  for  a  simulated 
gallium  arsenide  device  that  had  an  active  region  length  of  10  microns 
and  a  doping  density  of  1015/cm3.  The  space  charge  profiles  were 
developed  for  a  rather  sophisticated  model  of  the  cathode.  We  point 
out,  however,  that  the  use  of  this  model  is  not  necessary  for  simu¬ 
lating  the  qualitative  properties  of  the  instability  in  GaAs.  A  far 
simpler  model  can  be  used  to  generate  these  results;  one  in  which  the 
cathode  field,  instead  of  being  dependent  on  current  density,  is  held 
fixed  (Shaw  et  al . ,  1969).  We  have  carried  out  simulations  with  this 
model  for  10  micron  long  devices  as  well  as  100  micron  long  devices. 
The  cathode  dependent  results  are  qualitatively  similar.  We  have  also 
performed  experiments  and  used  the  fixed  cathode  field  model  in  their 
interpretation.  The  experimental  results  are  shown  in  Figs .  10  through 
12;  and  the  summary  of  the  gallium  arsenide  results  is  shown  in  Fig.  13. 
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Fig.  9.  Reactive  circuit  computations  showing  the  presence  of  a 
permanent  residual  voltage  drop  in  the  vicinity  of  the 
anode  contact.  (a)  Current  density  versus  average  elec¬ 
tric  field,  (b)  multiple  exposures  of  the  electric  field 
versus  distance  profiles.  From  Grubin  (1976b). 


The  experimental  results  were  obtained  for  Ga As  samples  cut  from 
bulk,  doped  Monsanto  n-GaAs  wafers  and  homegrown  epitaxial  specimens. 
The  Monsanto  wafers  had  quoted  Hall  mobilities  usually  between  4000  anc 
5000  cnn’/V-sec  and  carrier  concentrations  between  5xl014  and  1016/cm3. 
The  contacts  were  either  tin/nickel/tin  or  gold/germanium/nickel. 

Figure  10  shows  experimental  current  voltage  and  voltage  versus 
distance  profiles  for  a  sample  whose  active  region  boundary  field  was 
below  the  threshold  field  for  negative  differential  mobility  (Solomon 
£jt  al_.  ,  1975).  To  assure  this  we  took  the  expedient  approach  of 
removing  the  active  region  from  the  influence  of  the  metal  contact  by 
substantially  reducing  the  cross  sectional  area  of  a  large  region 
between  the  contacts.  This  sample  and  samples  of  this  type,  when  in  a 
resistive  circuit,  exhibited  oscillatory  characteristics  that  were 
more  closely  related  to  doping  inhomogenieties  than  to  the  transit¬ 
time  between  the  contacts. 

Figure  11  shows  experimental  results  for  a  device  that  sustained  a 
transit-time  oscillation.  We  see  the  enhanced  sublinearity  in  the 
current  versus  voltage  relation.  This  sublinearity,  prior  to  the 
instability  is  due  to  the  increased  voltage  drop  at  the  cathode  as 
the  bias  is  increased.  This  increasing  voltage  drop  is  revealed  in 
the  probe  measurements.  We  estimate  a  cathode  field  prior  to  an 
instability  of  approximately  6.2kV/cm. 

In  Figure  12  we  show  measurements  for  a  device  that  did  not  exhibit 
any  time  dependent  behavior.  Here  one  notices  extreme  saturation  in 
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Fig.  10.  (a)  Current  versus  voltage  for  a  low  cathode  boundary 

field  GaAs  sample.  The  probing  point  is  indicated  by 
the  closed  circle,  (b)  the  probed  voltage  versus  dis¬ 
tance  measurements.  From  Solomon  et  al.  (1975). 


current  at  high  bias  levels,  while  at  low  bias  levels  the  current 
voltage  relation  is  linear.  The  voltage  versus  distance  measurements 
show  a  very  large  cathode  voltage  drop  and  a  cathode  field  which  we 
estimate  to  be  greater  than  60  kV/cn.  Vie  show  measurements  in  both 
polarities  to  offer  evidence  that  the  large  voltage  drop  is  not  a 
consequence  of  a  high  resistance  cathode  region.  In  both  polarities 
we  see  that  there  is  not  a  large  anode  voltage  drop. 

The  results  of  the  cathode  boundary  dependence  of  two  terminal  nega¬ 
tive  differential  mobility  devices  is  summarized  in  Fig.  13.  Here 
for  the  same  device  parameters  we  plot  current  density  versus  average 
electric  field  (sample  voltage/sample  length)  for  different  fixed 
values  of  cathode  field.  We  superimpose  the  equilibrium  drift  velo¬ 
city  versus  electric  field  relation.  The  properties  of  the  insta¬ 
bility  fall  into  three  categories.  At  low  Ec  (c  irve  A),  the  j-  E 
curve  is  essentially  linear  up  to  the  field  at  peak  velocity.  At 
this  point  several  things  may  occur.  The  possibilities  depend  in  a 
detailed  way  on  the  characteristics  of  the  GaAs  and  the  circuit. 
Domains  may  be  nucleated  at  large  doping  fluctuations  in  the  bulk 
and  lead  to  transit-time  oscillations.  A  more  interesting  possi¬ 
bility  is  that  after  a  domain  is  nucleated  and  reaches  the  anode  it 
may  remain  there  in  a  stationary  field  configuration  and  the  current 
will  saturate.  In  a  resonant  circuit  the  field  may  oscillate  at  the 
resonant  frequency.  For  Ec  within  the  shaded  region,  the  current 
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Fig.  11.  Probed  sample  voltage  versus  distance  for  and  inter¬ 
mediate  cathode  boundary  field,  at  the  bias  point 
indicated  in  (b) .  (b)  Current  versus  voltage,  (c) 

instability  current  versus  time  response.  From 
Solomon  et  al.  (1975). 


shows  a  departure  from  linearity  due  to  the  appearance  of  a  large 
cathode  voltage  drop.  At  threshold  the  current  switches  along  the 
load  line.  Threshold  is  controlled  by  the  cathode  field  and  not  by 
the  threshold  field  for  negative  differential  mobility.  Threshold 
occurs  before  the  bulk  enters  the  region  of  negative  differential 
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Fig.  12.  Probed  sample  voltage  versus  distance  for  a  high  cathode 
field  sample  in  (a)  negative  and  (b)  positive  polarities. 
Current  versus  voltage  curves  for  (c)  negative  and  (d) 
positive  polarities.  From  Solomon  et  al.  (1975). 


mobility.  For  high  values  of  cathode  field  the  current  voltage  char¬ 
acteristic  is  linear  only  at  low  bias.  At  high  bias  the  curve  satu¬ 
rates  at  the  current  density  js  =  Noevs.  The  departure  from  line¬ 
arity  corresponds  to  the  appearance  of  a  large  cathode  voltage  drop 
and  transit-time  oscillations  in  a  resistive  circuit  do  not  occur. 

In  a  resonant  circuit  the  sample  may  sustain  some  oscillatory  acti¬ 
vity  for  a  certain  range  of  boundary  field.  We  note  that  in  all 
cases  where  an  oscillation  occurs,  stable  domain  propagation  requires 
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Fig.  13.  The  equilibrium  velocity  electric  field  curve  and  the 

computer  simulated  current  density  versus  average  elec¬ 
tric  field  curve  for  various  fixed  values  of  cathode 
boundary  field  Ec .  For  the  simulation  the  NDM  element 
is  in  the  circuit  shown  in  the  inset.  For  curve  A,  Ec = 
0.0.  The  values  for  curve  Bj  and  B?  are  indicated  by 
the  arrows,  for  curve  C,  Kc  =  24kv/cm.  The  simulated 
element  had  an  active  region  length  of  100  microns  and 
a  doping  level  of  lO'-Vcrr  .  The  right  and  left  hand^ 
ordinates  are  related  by  jn  =  Noev  and  vs  =  0.8  x  10 
cm/sec.  From  Shaw  et  al .  (1969). 


that  the  current  level  exceed  the  saturated  drift  current,  sometimes 
referred  to  as  the  sustaining  current. 

The  above  discussion  appears  to  present  a  rather  complete  phenomeno 
logical  description  of  the  contact  and  circuit  dependence  of  gallium 
arsenide  two  terminal  devices.  The  model  has  also  been  applied  to 
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the  ind ium  phosphide  device.  Here  the  presently  known  contact  behav¬ 
ior  is  somewhat  richer  than  that  for  gallium  arsenide  and  while  (2) 
and  (3)  can  be  used  to  simulate  the  InP  behavior,  the  fixed  cathode 
field  approximation  is  in  some  cases  inadequate.  In  particular,  while 
a  number  of  classes  of  devices  can  be  explained  with  contact  curves 
represented  by  Figs.  3  through  5,  the  appearance  of  very  high  effi¬ 
ciency  oscillations  (Colliver  ej:  al . ,  1973)  with  prethreshold  current 
voltage  curves  similar  to  that  of  Fig.  8,  requires  jc  (Ec) curves  simi¬ 
lar  to  that  of  #4  in  Fig.  2  for  its  simulation.  This  curve  is  quali¬ 
tatively  different  from  the  curves  discussed  earlier,  for  as  the 
range  of  current  oscillations  the  cathode  conduction  current  is 
approximately  constant,  while  the  cathode  field  can  almost  freely 
move  about.  A  plot  of  the  electric  field  profiles  during  the  course 
of  the  oscillation  is  shown  in  Fig.  14.  We  note  that  immediately 
downstream  from  the  cathode  contact  the  electric  field  profile  is 
approximately  uniform.  We  now  turn  to  the  subject  of  three  terminal 
devices . 


THREE  TERMINAL  DEVICES 

We  are  interested  in  the  hot  electron  contributions  to  the  operation 
of  gallium  arsenide  field  effect  transistors  (FET) .  The  FET  in  its 
simplest  form  is  a  semiconductor  slab  with  three  terminals.  Two  of 
these  are  usually  low  resistance  contacts,  while  the  third  is  either 
a  Schottky  contact  or  a  PN  junction  with  an  accompanying  region  of 
charge  depletion.  For  a  unipolar  conduction  device,  operation  is 
based  on  modulation  of  the  depletion  region  which  is  usually  accom¬ 
plished  by  changes  in  the  gate  bias.  Small  and  large  signal  (Colli¬ 
ver  et  al.,  1973)  gain  are  possible.  Figure  15  is  a  sketch  of  the 
device  and  connecting  lumped  elements,  whose  operation  we  have  simu¬ 
lated.  Our  (two  dimensional)  numerical  simulations  and  experiments 
(Grubin  aJL ,  1979)  demonstrate  that  the  GaAs  FET  can  be  placed  in 
one  of  two  groups  as  determined  by  the  ratio 

K _ Gate  voltage  at  cutoff _ 

Drain  voltage  at  the  onset  of  current 
saturation  for  zero  gate  voltage 

Devices  with  K  greater  than  unity  sustain  current  oscillations, 
whose  origin  lies  in  the  presence  of  negative  differential  mobility  in 
the  semiconductor.  Those  with  K  approximately  equal  or  less  than  unity 
are  electrically  stable.  We  point  out  that  in  the  original  FET  dis¬ 
cussion,  as  given  by  Shockley  (1952),  K  was  necessarily  equal  to 
unity.  The  situation  with  K  greater  than  unity  is  summarized  in  Fig. 
16a,  while  that  for  K  approximately  equal  to  unity  is  summarized  in 
Fig.  16b  (Grubin,  1977).  In  Fig.  16  we  sketch  the  current  voltage 
relation  for  a  gallium  arsenide  FET  with  the  parameters  listed  in  the 
caption.  For  reference  we  have  drawn  the  velocity  electric  field 
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Fig.  14.  Reactive  circuit  calculations  for  the  cathode  conduc¬ 
tion  current  curve  #4  of  Fig.  2.  The  inset  shows  the 
pre-instability  current  density  versus  average  elec¬ 
tric  field  curves  for  the  element  in  a  resistive  cir¬ 
cuit.  From  Shaw  et  al.  (1979). 


relation  for  GaAs  scaled  to  the  current  and  voltage  parameters.  The 
first  point  to  note  is  that  the  current  levels  do  not  approach  the 
peak  current  associated  with  GaAs.  This  is  a  consequence  of  the 
additional  resistance  supplied  by  the  gate  region  as  well  as  velocit 
limitations.  We  recall  from  our  earlier  discussion  of  two  terminal 
GaAs  devices  that  sublinearity  in  the  GaAs  I-V  relation  is  often 
accompanied  by  a  current  instability.  For  the  wider  channel  device 
shown  in  Fig.  16  similar  behavior  occurs.  The  instability  is  repre¬ 
sented  by  the  dashed  part-  of  the  I-V  curves  of  Fig.  16a.  The  x's  in 
the  diagram  represent  ave./.ge  current  and  voltage  values  for  the 
instability,  and  the  presence  of  negative  conductance  is  due  to  the 
dynamic  propagating  domain.  The  closed  circles  in  Fig.  16  represent 
stationary,  time-independent  points,  and  we  note  that  when  there  is 
an  instability  it  is  surrounded  by  regions  of  nonzero  gate  and  drain 
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Fig.  15.  Schematic  representation  of  a  Schottky  -  gate  field 
effect  transistor. 

bias  for  which  there  is  no  time-dependent  behavior.  We  now  illus¬ 
trate  the  space  charge  distribution  associated  with  the  x's  and 
closed  circles. 

The  internal  distribution  of  charge  and  current  associated  with  the 
current  and  potential  levels  *  1  *  and  '2'  of  Fig.  16a  are  shown  in 
Fig.  17  (Grubin  and  McHugh,  1978).  Column  a  shows  a  set  of  current 
density  streamlines  through  the  device.  The  length  of  each  stream¬ 
line  is  proportional  to  the  magnitude  of  the  vector  current  density 
at  the  point  in  question.  The  maximum  length  of  the  individual  x 
and  y  components  before  overlap  is  jp  =  Noevp  where  vp  is  the  peak 
carrier  velocity.  We  see  in  both  cases  that  the  current  density  is 
greatest  under  the  gate  contact  as  required  by  current  continuity. 

For  the  higher  bias  case,  the  current  density  under  the  gate  region 
is  at  least  as  great  as  jp  and  velocity  limitation  introduces  carrier 
accumulation.  The  density  of  charge  particles  in  the  FET  is  gener¬ 
ally  nonuniform  and  column  b  of  Fig.  17  is  a  projection  of  this  dis¬ 
tribution  as  it  may  relate  to  the  current  density  profile  of  column 
a.  We  point  out  that  the  particle  density  increases  in  the  downward 
direction.  The  first  frame  shows  a  region  of  charge  depletion 
directly  under  the  gate  contact.  The  second  frame  shows  the  forma¬ 
tion  of  a  weak  stationary  dipole  layer  under  the  gate  contact.  We 
see  from  these  diagrams  that  the  x-component  of  electric  field  has 
reached  the  NDM  threshold  field  value  at  the  drain  edge  of  the  gate 
contact.  We  point  out  that  an  analysis  of  the  potential  distribution 
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Fig.  17.  The  internal  distribution  of  current  and  charge  corre¬ 
sponding  to  the  current  and  voltage  levels  represented 
by  numbers  *1'  and  *2’  of  Fig.  16.  Note  that  the  par¬ 
ticle  density  surrounding  the  nonuniform  distribution 
is  uniformly  distributed  within  the  source-gate  region 
and  the  gate-drain  region.  From  Grubin  and  McHugh 
(1978). 


shows  that  at  the  high  bias  levels  most  of  the  potential  drop  is 
under  the  gate  contact. 

We  now  consider  the  presence  of  an  instability  (Grubin  and  McHugh, 
1978).  We  recall  that  for  two  terminal  devices  the  instability  was 
determined  by  the  value  of  electric  field  at  the  cathode  boundary 
region  and  that  the  threshold  current  density  for  the  instability 
could  be  anywhere  between  js  and  jp,  where  jg  is  the  current  density 
associated  with  the  high-field  saturated  drift  velocity.  In  three 
terminal  devices,  the  initiation  of  a  domain  instability  generally 
occurs  under  the  gate  contact.  The  instability  occurs  at  a  value 
of  current  density  approximately  equal  to  jp.  Figure  18  provides 
a  dramatic  representation  of  the  FET  instability,  where  we  see  the 
nucleation,  propagation  and  recycling  of  a  high  field  domain.  The 
sequence  of  events  associated  with  the  instability  is  as  follows: 
Domain  growth  under  the  gate  is  accompanied  by  an  increase  in  poten¬ 
tial  across  the  device.  A  corresponding  decrease  in  current  occurs 
throughout  the  device  and  circuit,  as  constrained  by  the  dc  load 
line.  As  the  current  decreases,  carriers  with  velocities  below  that 
of  the  peak  velocity  enter  the  accumulation  layer  which  subsequently 


ection  of  the  time  dependent  particle  density  when 
nstability  occurs.  The  parameters  at  which  this 
irs  are  represented  by  the  number  '3'  in  Fig.  16. 
i  Grubin  and  McHugh  (1978). 
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begins  to  detach.  The  domain  spreads  as  it  leaves  the  gate  region 
and  it  settles  into  a  value  of  current  density  somewhat  in  excess  of 
that  associated  with  the  saturated  drift  velocity  of  the  electrons. 
Prior  to  reaching  the  drain  contact  the  domain  dynamics  appear  to  be 
one  dimensional. 

The  numerical  situation  at  high  drain  bias  levels  and  zero  gate 
bias  is  similar  to  that  of  two  terminal  devices.  Namely  the  oscilla¬ 
tions  cease  and  are  accompanied  by  the  presence  of  an  accumulation 
layer  extending  from  the  gate  to  the  drain  contact.  This  is  illustrat¬ 
ed  in  Fig.  19.  The  region  of  space  charge  accumulation  is  accompanied 
by  a  large  nonuniform  potential  drop  in  the  gate  to  drain  region  with 
resulting  high  values  of  electric  field  and  electrons  traveling  at 
their  saturated  drift  values.  This  region  of  space  charge  accumula¬ 
tion  is  qualitatively  similar  to  the  anode-adjacent  domain  discussed 
in  connection  with  Fig.  9. 

Our  discussion  of  two  terminal  devices  indicates  that  large  values 
of  cathode  boundary  fields  limit  the  downstream  carrier  velocity  to 
values  below  the  saturated  drift  value.  Instabilities,  if  present, 
are  damped.  An  analogous  situation  occurs  with  three  terminal  devices 
when  a  large  negative  bias  is  applied  to  the  gate  contact.  For  this 
case  the  depletion  layer  moves  toward  the  bottom  of  the  channel  and 
the  source-drain  current  is  low  (Grubin  and  McHugh,  1978).  A  large 
potential  drop  is  present  under  the  gate  contact  resulting  in  the 
formation  of  a  high  field  domain.  As  seen  in  Fig.  20,  the  domain 
consists  of  an  accumulation  region  surrounded  by  a  depletion  zone. 

The  experimental  situation  with  regard  to  hot  electron  contribu¬ 
tions  to  the  FET  may  be  summarized  by  turning  to  Fig.  21  (Grubin  and 
McHugh, 1978)  where  we  show  the  dc  current  voltage  characteristics  for  a 
device  whose  parameters  are  listed  in  the  caption  and  for  which  K>1. 
The  cross  hatched  region  shows  the  range  over  which  current  oscilla¬ 
tions  were  observed  as  detected  by  a  spectrum  analyzer.  Oscillations 
were  detected  over  a  range  of  6  to  40  GHz.  In  all  of  our  experiments 
in  which  oscillations  were  observed,  the  instability  occurred  beyond 
the  'knee'  of  the  current  voltage  relation.  Prior  to  the  oscillation 
the  drift  velocity  in  the  gate  to  drain  region  was  approximately 
equal  to  or  greater  than  the  saturated  drift  velocity  of  the  carriers. 
Saturation  in  current  occurred  when  the  average  field  within  the  con¬ 
ducting  channel  under  the  gate  contact  was  approximately  equal  to  the 
threshold  field  for  negative  differential  mobility  (Grubin,  1977). 

The  instability  was  suppressed  either  by  going  to  a  sufficiently  high 
drain  bias  or  by  going  to  a  sufficiently  large  and  negative  gate  bias, 
such  that  in  the  latter  case  the  gate  to  drain  channel  current  density 
was  below  js. 

The  above  discussion  demonstrates  that  in  wide  channel  devices  hot 
electron  effects  in  three  terminal  devices  are  similar  to  that  of 
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X  GATE 


Fig.  19.  Distribution  of  current  and  particle  density  for  the  case 
represented  by  '4'  of  Fig.  16.  The  instability  has 
ceased.  Note  that  the  current  density  within  the  gate  to 
drain  region,  where  there  is  particle  accumulation,  is 
accompanied  by  current  density  levels  that  are  approxi¬ 
mately  equal  to  the  saturated  drift  current.  From  Grubin 
and  McHugh  (1978). 


two  terminal  devices.  A  simple  intuitive  picture  of  this  commonality 
arises  after  making  a  one-to-one  correspondence  between  the  potential 
drop  in  that  portion  of  the  conducting  channel  that  is  under  the  gate 
contact  and  a  phenomenological  cathode  boundary  for  two  dimensional 
devices. 

The  situation  with  narrow  channel  devices  is  somewhat  less  dramatic 
from  the  point  of  view  of  domain  instabilities.  Generally  because 
of  the  nonlinear  distribution  of  potential  under  the  gate  contact, 
very  thin  devices  sustain  current  densities  significantly  below  js 
in  the  region  between  the  gate-to-drain  contact.  As  in  two  terminal 
devices,  this  is  sufficient  to  prevent  an  instability. 


DEVICES 


ribution  of  current  and  particle  density  for  the  case 
;sented  by  '5*  in  Fig.  16.  The  instability  has  ceased, 
the  large  dipole  layer  within  the  vicinity  of  the  gate 
>ct,  and  the  fact  that  the  current  density  between  the 
and  drain  region  is  very  small.  From  Grubin  and 
Sh  (1978). 
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Fig.  21.  Measured  dc  current  voltage  data  for  a  gallium  arsenide 

field  effect  transistor.  The  device  had  a  nominal  doping 
level  of  1017/cm3,  and  epitaxial  layer  thickness  of  3000 
±  500°A,  a  source-drain  spacing  of  approximately  8.5 
microns  and  a  gate  length  of  3.0  microns.  The  mobility 
measured  from  adjacent  samples  varied  from  3000  to  4000 
cn7/v-s.  The  cross-hatched  region  of  the  figure  denotes 
the  bias  range  for  the  instability,  as  detected  by  spec¬ 
tral  analysis.  Details  of  the  experimental  results  are 
discussed  in  Grubin  et  al.  (1979),  from  which  this  data 
is  taken. 
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INTRODUCTION 

The  original  impetus  for  much  of  the  hot  electron  work  that  has 
taken  place  in  the  last  decade  was  the  need  to  understand  the  con¬ 
tributions  of  hot  electrons  to  semiconductor  devices.  The  turning 
point  in  this  area  of  study  occurred  when  Gunn  (1964a)  published  his 
results  showing  the  presence  of  spontaneous  and  coherent  oscillations 
in  gallium  arsenide  and  indium  phosphide  two  terminal  devices. 

Today,  of  course,  hot  electron  phenomena  are  being  studied  for  their 
own  interest,  and  with  ends  that  are  ostensibly  independent  of  device 
considerations.  But,  in  fact,  the  two  areas  cannot  be  separated. 
Rather  the  emphasis  is  different,  with  hot-electron-device  studies 
concentrating  on  situations  where  the  electric  field  sustains  large 
and  fast  transients  and  where  the  field  is  highly  nonuniform. 

Typically,  when  it  is  of  interest  to  find  an  experimental  situa¬ 
tion  that  best  illustrates  hot  electron  device  effects,  the  results 
of  Gunn's  early  studies  are  shown.  In  Fig.  1,  we  display  results 
(Gunn,  1964a)  of  a  current  instability  in  gallium  arsenide.  Here  we 
show  sampled  current  -  voltage  characteristics  of  a  device  thirty- 
six  microns  long.  Mathematically,  these  results  are  plots  of  current 
as  a  function  of  voltage  with  time  held  constant.  At  low  voltages 
the  I-V  relationship  is  nearly  linear.  At  higher  voltages  the  dc 
conductance  first  decreases  slightly;  then  at  some  value  of  voltage 
there  occur  wide  differences  in  current  between  successive  samples. 
This  appears  as  a  large  scatter  of  dots  representing  individual 
samples.  In  a  subsequent  series  of  experiments  Gunn  (1964b)  estab¬ 
lished  that  these  fluctuations  in  current  were  consequences  of  co¬ 
herent  oscillations  associated  with  a  moving  distribution  of  poten¬ 
tial  within  the  semiconductor. 
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VOLTAGE  (VOLTS) 


Figure  1.  Sampled  current-voltage  characteristics  of  an  n-type 
GaAs  specimen  The  active  region  lengths  and  scales 
are  shown  in  the  figure.  From  Gunn  (1964a)  with  per¬ 
mission. 


Gunn's  experiments  demonstrated  that  the  semiconductor  gallium 
arsenide  was  capable  of  converting  dc  to  ac  power.  Our  present  ex¬ 
planation  (Kroemer,  1964)  for  this  conversion  is  that,  through  the 
mechanism  of  electron  transfer  (Ridley  and  Watkins,  1961)  from  the 
central  to  satellite  valleys  of  the  conduction  band,  gallium  arsenide 
possesses  a  region  of  negative  differential  mobility  (Hilsum,  1962) 
(see  Fig.  2).  As  discussed  by  Ridley  (1963),  the  presence  of  nega¬ 
tive  differential  mobility  can  lead  to  unstable  localized  space 
charge  regions.  The  observations  of  Gunn  are  then  explained  as 
arising  from  the  periodic  nucleation  and  disappearance  of  traveling 
space  charge  domains  (Kroemer,  1964). 

Figure  2,  which  is  the  drift  velocity  versus  electric  field 
curve  for  gallium  arsenide  (Butcher,  1967),  represents  what  most 
device  physicists  regard  as  the  basis  of  hot-electron  transport 
effects  in  devices.  Indeed,  very  sophisticated  numerical  programs 
have  been  designed  around  this  curve  and  used  to  model  the  spatial 
and  time-dependent  dynamics  of  electrons  in  devices.  Remarkable 
agreement  with  experiment  has  been  obtained  that  has  also  included 
explanations  of  fine  details  of  observations.  I  will  illustrate  some 
of  these  in  a  following  chapter.  Perhaps  the  most  noteworthy  point 
is  that  the  velocity  field  curve,  as  represented  by  Fig.  2,  derived 
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ELECTRIC  FIELD  (kV/cm) 

Figure  2.  Calculation  of  the  steady  state  drift  velocity  versus 
electric  field  curve  for  gallium  arsenide.  From 
Butcher  (1967)  Figure  3,  with  permisssion. 


under  steady-state  conditions  and  spatially-uniform  fields,  is  often 
used  to  explain  transient  phenomena  and  highly  nonuniform  situations. 
It  is  this  latter  use  that  is  currently  the  subject  of  increased 
scrutiny  and  will  occupy  the  rest  of  the  current  discussion.  These 
topics  properly  come  under  the  heading  of  transient  hot-electron 
effects  in  semiconductor  devices. 

Before  beginning  the  discussion  we  must  lay  a  few  ground  rules. 
We  will  provide  a  microscopic  description  of  carriers  in  a  semicon¬ 
ductor  device  in  terms  of  a  displaced  Maxwelliam 


f(p±,  1^)  “  A(ni/Ti3/,2)exp[Hi2(k-p1)2/2mikoT1] 


(1) 


The  dynamics  of  the  particles  are  then  obtained  from  moments  of  the 
Boltzmann  transport  equation  in  the  form  (Bosch  and  Thim,  1974) 
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In  the  above,  a^  represents  the  number  of  equivalent  satellite 
valleys,  n^,  p^  ancj  are  respectively  the  number  density,  the 
momentum  and  temperature  of  the  ith  valley.  Further,  we  have  ig¬ 
nored  the  spatial  contributions  to  the  moment  equations.  The  t's 
in  the  above  equation  are  scattering  times  for  intervalley  and  intra¬ 
valley  transfer  and  significant  effects  arise  because  of  important 
differences  between  the  energy  and  momentum  relaxation  time.  Fig.  3 
shows  a  set  of  scattering  times  for  gallium  arsenide  calculated  in 
1974  by  Bosch  and  Thim  (1974).  These  calculations  are  for  the  T-X 
ordering  in  gallium  arsenide  in  which  the  number  of  equivalent  satel¬ 
lite  valleys  is  equal  to  3.  Calculations  for  the  T-L  ordering 
with  four  equivalent  satellite  valleys  have  not  yet  been  published, 
although  several  groups  are  pursuing  the  problem.  VJe  note  that  at 
room  temperature,  the  momentum  scattering  time  for  the  central  valley 
is  significantly  smaller  than  the  energy  scattering  time.  The  sig- 
lificance  of  this  will  emerge  below  in  connection  with  overshoot 
’ elocity . 


Tl  4NSIENT  HOT  ELECTRON  EFFECTS  IN  SEMICONDUCTOR  DEVICES 
Vc.ocity  Overshoot 

If  a  system  of  electrons  is  subjected  to  the  combined  influence 
of  ai  electric  field  and  scattering  centers,  ignoring  for  the  moment 
elec  r on  transfer  as  a  scattering  mechanism,  then  the  drift  velocity 
of  th>  particles  asymptotically  approaches  the  steady  state  value 


v  =  -e<t>  F/m*. 


(5) 


SCATTERING  RATES  O/sec) 


1/t)  of  GaAs  ve 
The  symbols  a 
d  Thim  (1974),  \ 
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TIME 

Figure  4.  Approach  to  steady  state,  with  and  without  transient 
relaxation  effects.  Part  a  denotes  increasing 
velocity  and  part  b  denotes  decreasing  velocity  values. 


If  we  were  just  to  consider  the  momentum  balance  equation,  (2),  we 
could  be  satisfied  that  the  particles  would  approach  this  value  in  a 
time  approximately  equal  to  2<t>  seconds,  as  illustrated  in  Fig.  4. 
But  considering  the  momentum  balance  equation  alone  would  give  us  an 
incomplete  picture  of  events.  Energy  balance  tells  us  that  the 
electron  temperature  increases  with  increasing  electric  field  and 
departs  significantly  from  room  temperature  when  the  electric  field 
exceeds  a  threshold  value.  The  effect  of  the  increasing  electron 
temperature  is  to  decrease  the  average  collisioi.  time  and  to  decrease 
the  steady-state  velocity,  given  by  (5).  Now,  if  the  momentum  and 
energy  scattering  times  are  similar  in  value  then  b.->th  momentum  and 
energy  will  follow  changes  in  electric  field  at  approximately  the 
same  rate  and  the  solid  curve  of  Fig.  4  describes  the  approach  to 
steady-state.  On  the  other  hand  if  the  energy  scatter iiq  time  is 
significantly  longer  than  the  momentum  scattering  time,  the  average 
velocity  of  the  carriers  will  find  its  value  continually  'corrected' 
until  steady  state  in  the  energy  distribution  is  reached.  The  velo¬ 
city  will  then  relax  in  the  manner  shown  by  the  dashed  cum,  A 
similar  situation  may  be  expected  when  the  electric  field  is  de¬ 
creased,  for  here  it  also  takes  a  finite  time  for  the  electric  'ield 
to  decrease  and  for  the  electron  temperature  to  decrease. 
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We  may  now  ask:  where  do  we  expect  the  overshoot  transient 
effects  to  become  significant?  On  the  basis  of  the  above  discussion, 
substantial  overshoot  effects  occur  when  the  momentum  scattering  time 
is  much  shorter  than  the  energy  scattering  time.  Thus  we  would  ex¬ 
pect  to  detect  these  effects  when  the  time  rate  of  change  of  electric 
field  falls  somewhere  within  the  limits 


F  F 

th  >  d¥  >>  th 


dt 


(6) 


We  are  interested  in  the  upper  bound  of  the  inequality  and  for  this 
case  we  take  Tp  to  be  represented  by  the  momentum  scattering  time 
for  LO  intravailey  phonons.  Intervalley  scattering  is  important  as 
the  electron  temperature  is  elevated,  and  will  reduce  the  momentum 
scattering  time.  For  LO  phonon  intravailey  scattering  (Conwell, 
1967), 


1  FOe 

—  £ — - -  (7) 

x  (2m*k  T  ) 1/2 

P  op 

For  the  central  valley  of  GaAs,  with  F  =  5.6  x  103  V/cm  (Butcher, 
1967), 


—  %  4xl012/sec  (8) 

TP 


which  is  consistent  with  the  more  complete  calculations  of  Fig.  3 

'  3.5  x  1012/sec  for  temperatures  below  1000  K) .  Transient 

TPi  .  .  . 

overshoot  effects  are  then  expected  to  be  significant  when  the 

carriers  experience  changes  in  the  electric  field  of  the  order  of 

Fth  (3kV  at  the  doping  levels  of  1015  / cm3)  occurring  during  a  time 

interval  of  the  order  of;  or  somewhat  larger  than,  xp. 


In  Fig.  5  we  show  (Rees,  1969)  the  response  of  the  central 
valley  velocity,  the  satellite  valley  velocity,  the  satellite  valley 
population  and  the  drift  velocity,  when  the  electric  field  is  stepped 
down  from  6kV/cm  to  5kV/cm.  The  drift  velocity  is  computed  from  the 
equation 


n  v  +  a  n  v 
c  c  s  s  s 

n  +  n 
c  s 


(9) 


For  gallium  arsenide,  the  satellite  valley  scattering  times  are  sig¬ 
nificantly  shorter  than  the  central  valley  scattering  times  and  so 
the  satellite  velocity  readjusts  itself  to  the  change  in  electric 
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Figure  5.  Response  time  of  electrons  to  a  field  stepped  down  from 

6  to  5  kv/cm  at  time  t  =  0.  (a)  and  (b)  are  respectively 

the  central  and  satellite  valley  velocities,  (c)  is  the 
satellite  population  fraction,  (d)  is  the  drift  velocity 
of  the  carriers.  From  Rees  (1969)  with  permission. 
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field  in  less  than  10“ 13  sec.  The  central  valley  velocity  takes 
somewhat  longer  to  reach  its  minimum  value  and  to  then  recover. 

Here  v  falls  for  about  4x10“ 13  sec  and  then  begins  to  recover  over 
a  period  of  approximately  5  picoseconds.  The  transfer  of  electrons 
from  the  central  valley  to  the  satellite  val..-?y  is  also  slow.  The 
net  response  of  the  drift  velocity  shows  an  initial  drop  due  to  the 
central  and  satellite  valley  velocities  followed  by  a  recovery  due 
partly  to  the  central  valley  velocity  and  partly  due  to  differential 
electron  transfer.  The  final  velocity  exceeds  the  initial  velocity 
and  illustrates  the  presence  of  negative  differential  mobility. 

Transient  calculations  of  the  above  type  are  important  for 
determining  the  conditions  under  which  a  material  like  gallium 
arsenide  will  remain  an  active  device.  Rees  (1969),  using  the 
results  of  Fig.  5,  determined  the  frequency  dependence  of  the 
differential  mobility.  These  results  are  shown  in  Fig.  6.  For  this 
calculation,  the  Fourier  transform  of  the  drift  velocity  was  taken, 
yielding  v  (w) .  The  mobility  was  defined  as  the  ratio  of  v(w)/5.5 
kV/cm  (for  a  change  from  6kV/cm  to  5kV/cn);  5 . 5  kV/cra  is  the  average  of 
the  initial  and  final  values  of  the  electric  field.  What  we  notice 
here,  is  that  the  real  part  of  the  mobility  is  negative  at  fre¬ 
quencies  below  80  GHz.  It  is  positive  above  80  GHz,  reflecting  the 
short  time  response  of  the  central  valley  electrons. 


Spatial  Dependence 

The  above  arguments  are  for  uniform  fields.  We  can  use  this 
information  to  estimate  the  extent  to  which  spatial  changes  in  the 
electric  field  will  affect  the  transient  response  of  the  carriers. 
To  carry  out  this  argument,  we  must  first  calculate  the  mean  free 
path  of  electrons  injected,  with  zero  initial  drift  velocity,  into 
a  region  of  high  electric  field.  The  mean  free  path  is  calcualted 
prior  to  thermalization  assuming  a  scattering  time  given  by  (7). 
From  the  momentum  balance  equation,  we  have 


v 


p/m*  =  -(eTp/m*)F 


1  -  exp  (-t/r  ) 

P  J 


(10) 


and  in  a  time  t  =  2tp,  v  will  have  reached  approximately  87%  of  its 
peak  value.  During  this  time  the  drifting  electrons  will  have 
travelled  a  distance 


i  =  (ex2/m*  )  F  (1  -  e~2)  ,  (11) 

o  p 


before  thermalization  occurs.  We  estimate  that  if  the  spatial  varia¬ 
tion  in  the  electric  field  is  such  that  large  changes  in  the  electric 
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Figure  6.  Frequency  dependence  of  the  differential  mobility  for  a 
bias  field  of  5.5  kV/cm.  Pp  and  Pp  denotes  respectively 
the  components  in  phase  and  quadrature  with  the  ac 
field.  From  Rees  (1969)  with  permission. 


field  occur  in  distances  less  than  that  given  by  (11)  ,  then  signi¬ 
ficant  relaxation  effects  may  be  expected  to  take  place.  As  an 
example,  for  an  electric  field  of  5kV/cm  and  a  cold  mobility  of 
6500cm2 /V-sec,  electrons  in  the  central  valley  will  have  attained 
a  velocity  of  2.8xl07cm/sec  in  a  time  2xp,  and  will  have  travelled 
a  distance  of  approximately  0.1  microns  before  the  steady  state 
value  is  reduced  by  increases  in  the  electron  temperature.  Figure 
7  illustrates  this  situation  (Ferry  and  Barker,  1979).  Here  the 
velocity  response  of  electrons  in  GaAs  is  shown  for  a  field  of 
5kV/cm.  These  calculations  include  intervalley  transfer  and  con¬ 
tributions  from  the  satellite  valley,  for  a  T-L  ordering. 

The  above  results  suggest  that  transient  overshoot  effects  will 
be  important  when  the  electric  field  undergoes  a  change  of  the  order 
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Figure  7.  The  velocity  response  of  electrons  in  GaAs  for  an  elec¬ 
tric  field  of  5  kV/cm.  From  Ferry  and  Barker  (1979), 
with  permission. 


of  Fth  (or  greater)  within  a  distance  of  fc0.  In  other  words  over¬ 
shoot  effects  will  be  pronounced  when 


dF 

dx 


> 


(12) 


Using  Poisson's  equation  this  can  be  translated  into  the  condition 


N 


0 


(13) 


where  t  is  the  dielectric  relaxation  time.  At  Nq  =  1015/cm3,  the 
right-hand-side  of  the  above  equation  is  equal  approximately  to  A. 
Thus,  when  the  ratio  of  N/Nq  exceeds  5  we  may  expect  relaxation 
effects  to  enter  in  an  important  way.  In  an  early  discussion. 
Butcher,  et  al.,  (1967)  estimated  the  maximum  value  of  N/No  for  the 
accumulation  portion  of  a  stable  travelling  domain.  As  is  well 
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known  from  domain  theory,  as  the  outside  field  begins  to  decrease 
the  domain  increases  in  size  and  the  accumulation  portion  of  the 
domain  increases  dramatically.  Butcher,  et^  a^l. ,  (1967)  found  that 
when  the  outside  field  is  below  2kV/cm,  the  accumulation  ratio 
exceeds  5,  and  that  for  outside  fields  approaching  the  sustaining 
field  the  ratio  is  substantially  higher.  The  implication  is  that 
a  sizable  portion  of  the  accumulation  layer  is  not  accurately 
modeled  by  steady-state  velocity-electric-field  relationships. 
Rather,  transient  effects  must  be  included  for  a  complete  picture. 
The  degree  to  which  this  particular  transient  contribution  will 
affect  the  output  at  the  terminals  of  device  has  not  yet  been  deter¬ 
mined  . 

As  we  go  to  higher  carrier  concentrations,  x  decreases  and  the 
criteria  expressed  by  (13)  becomes  more  restrictive.  For  example 
at  1017  cm~3,  accumulation  layers  in  excess  of  4%  of  Ng  will  intro¬ 
duce  significant  relaxation  effects. 


Relaxation  Effects  Within  a  Travelling  Accumulation  Layer:  The 
Jones-Rees  Effect 


The  discussion  above  for  nonuniform  fields  was  based  on  con¬ 
cepts  derived  from  a  uniform  field  analysis.  While  this  approach 
is  useful  for  estimating  some  of  the  relaxation  contributions,  it 
is  not  wholly  satisfying.  Jones  and  Rees  (1973)  have  examined  some 
of  these  limitations.  For  this  discussion,  we  refer  to  the  accumu¬ 
lation  layer  profile  sketched  in  Fig.  8. 

The  accumulation  layer  profile  consists  of  two  regions,  a 
high-field  and  a  low-field  region  separated  by  a  region  of  sub¬ 
stantial  charge  accumulation.  The  downstream  portion  of  the  accu¬ 
mulation  layer  consists  mainly  of  satellite  electrons.  In  analyzing 
the  propagation  of  the  accumulation  layer,  Jones  and  Rees  (1973) 
observed  a  phenomena  that  appeared  in  all  of  their  simulations. 

The  fall  of  the  satellite  electron  density  from  the  high-field 
region  to  the  low-field  region  was  quite  sharp.  This  behavior  was 
not  expected.  In  the  spatially  uniform  analysis  in  which  the 
field  decreased  rapidly  to  a  value  below  threshold,  a  large  fraction 
of  the  satellite  electrons  was  initially  retained.  About  5-10  psec 
was  required  for  their  thermalization  (see  Fig.  5).  Since  the  layer 
velocity  in  their  computer  simulations  was  above  the  satellite  elec¬ 
tron  velocity,  a  long  trailing  edge  of  satellite  electrons  was 
expected.  But  in  fact  this  did  not  occur.  Rather  a  time  constant 
of  the  order  of  1  psec  represented  the  rate  of  transfer.  The 
mechanism  proposed  by  Jones  and  Rees  (1973)  to  explain  the  differences 
between  the  relaxation  time  for  uniform  fields  and  nonuniform  fields 
is  described  below. 
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.  Profiles  of  electric  field,  total  electron  density  and 
satellite  valley  density,  and  drift  velocity  for  a 
mature  propagating  domain.  From  Jones  and  Rees  (1973) 
with  permission. 
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Figure  9.  k-space  representation  of  carrier  cooling  behind  an 

accumulation  layer.  'A'  electrons  scattered  to  T  valley- 
backward  momenta  states  lose  energy  and  fall  behind  the 
layer.  'B'  electrons  scattered  to  T  valley-forward 
momenta  states  gain  energy  and  remain  within  the  layer. 


Figure  9  illustrates  possible  k-space  trajectories  for  elec¬ 
trons  initially  within  the  layer.  We  will  assume,  for  illustrative 
purposes,  that  the  electrons  are  initially  occupying  the  satellite 
valley  states.  These  electrons  rapidly  scatter  between  the  satel¬ 
lite  and  central  valley  states  with  high  energy  and  momenta.  Some 
of  these  electrons  will  scatter  into  k-states  whose  direction  is 
opposite  to  that  of  the  moving  distribution.  Those  scattered 
electrons  will  lose  energy  and  will  not  scatter  back  into  the 
satellite  valley.  They  are  then  concentrated  behind  the  accumula¬ 
tion  layer  and  account  for  the  lower  satellite  valley  population. 

In  contrast,  electrons  scattered  to  high  forward  momenta  states 
gain  energy  from  the  field  and  are  concentrated  in  the  front  of  the 
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accumulation  layer.  At  least  part  of  the  propagating  accumulation 
layer  is  due  to  this  scattering  mechanism.  Thus  the  traveling 
accumulation  layer  segregates  the  electrons  into  two  classes  and, 
in  leaving  a  cold  distribution  behind  the  layer,  it  overrides  the 
tendency  to  leave  electrons  heated,  which  characterizes  a  rapid 
fall  of  field  in  a  spatially  uniform  distribution. 


CONCLUSIONS 

In  the  above,  I  have  outlined  several  hot  electron  transient 
effects  in  transferred  electron  device.  As  of  this  writing,  it  is 
not  at  all  clear  what  role  these  specific  effects  will  have  on  the 
operation  and  output  of  small  electron  devices.  What  is  clear  is 
that  it  is  incorrect  to  try  to  discuss  the  behavior  of  very-small, 
highly-doped  devices  without  examining  transient  relaxation  pheno¬ 
mena. 
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I  Introduction 

THE  SEMICONDUCTOR  gallium-arsenide  is  a 
widely  used  field-effect  transistor  material.  It  is  an 
unusual  one  in  that  it  possesses  a  region  of  negative 
differential  mobility  (NDM)  and  two  characteristic  veloci¬ 
ties:  a  peak  velocity  prior  to  the  onset  of  NDM  and  a 
saturated  drift  velocity  following  the  NDM  region.  The 
effects  of  NDM  are  currently  being  mapped  out  and  it  is 
developing  that  NDM  separates  devices  into  two  cate¬ 
gories:  those  that  sustain  spontaneous  oscillations  and 
those  that  do  not  (1).  Devices  sustaining  spontaneous 
oscillations,  whose  origin  is  due  to  transiting  and  recycling 
space  charge  layers,  generally  have  cutoff  voltages  greater 
than  the  drain-voltage  value  at  the  onset  of  current  satura¬ 
tion.  Those  that  do  not,  possess  cutoff  voltages  approxi¬ 
mately  equal  to  saturation  drain  voltage.  The  latter  often 
yield  low  power  and  are  useful  as  switching  devices.  In  the 
last  few  years  there  has  been  an  increasing  interest  in 
GaAs  switches  for  integrated  circuit  applications.  The 
advantages  for  gallium-arsenide  are  thought  to  be  very 
high  speeds  and  low-power-delay  products.  The  purpose 
of  this  study  is  to  examine  these  conclusions,  and  in  doing 
so  we  highlight  the  roles  that  the  low  field  mobility  and 
the  high  field  saturated  drift  velocity  have  on  the  switch¬ 
ing  times.  To  place  the  results  in  perspective  we  do  similar 
calculations  for  the  semiconductor  silicon. 


II.  Numerical  Simulation  Modeling 

The  switching  calculations  discussed  below  are  for  four 
different  elements,  including  silicon.  For  the  GaAs  FET 
we  examine  only  those  devices  that  do  not  sustain  sponta¬ 
neous  oscillations.  To  highlight  the  importance  of  the  high 
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Fig.  1.  Drift  velocity  versus  electric -field  relauon  used  in  the  calcula¬ 
tions.  The  galhum-arscnide  results  are  taken  from  [7,  fig.  3). 


field  drift-velocity  values  we  also  present  results  of  calcu¬ 
lations  wilb  nonlinear  elements  whose  low  field  mobility  is 
equal  to  that  of  gallium-arsenide  but  whose  high  field 
velocity  saturates  without  NDM.  In  one  case  the  drift 
velocity  saturates  to  the  saturated  drift  velocity  of 
gallium-arsenide.  In  the  second  case  it  saturates  to  the 
peak  velocity.  The  velocity  electric-field  characteristics  for 
these  calculations  are  displayed  in  Fig.  1,  where  we  desig¬ 
nate  the  high-saturated  drift-velocity  element  HiV,  and 
the  lower  saturated  drift-velocity  element  Lo  V. 

The  calculations  with  the  Lo  V  and  HiF  elements  are 
more  than  academic.  It  is  now  generally  agreed  that  these 
nonlinear  elements  lead  to  dipole  formation  within  the 
gate  to  drain  region  [2],  [3]  and  that  at  large  values  of  bias 
the  electrons  within  this  region  are  travehng  at  their 
saturated  drift-velocity  values.  If  this  is  the  case  then  it 
has  been  argued  that  the  saturated  drift  velocity  is  the 
primary  factor  at  high-bias  levels  for  determining  the 
switching  speeds. 

In  the  following  calculation  we  use  an  already  devel¬ 
oped  [4]  large-signal  numerical  simulation  designed  to 
determine  the  transient  behavior  of  nonlinear  semicon¬ 
ducting  FETs.  The  simulation  results  in  self-consistent 
solutions  to  Poisson's  equation,  the  equation  of  continuity 
and  the  external  circuit  equations.  The  solutions  are  for 
two  dimensions  plus  time  with  the  electrons  described  by 
a  specific  drift  velocity  versus  electric-field  curve  shown  in 
Fig.  1  and  a  diffusion  electric-field  curve  shown  in  Fig.  2. 
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Fig.  2.  Diffusiviiy  versus  electric  field  used  in  the  calculations.  The 
gallium-arsenide  results  are  taken  from  |7.  fig.  4). 


placement  of  the  contacts.  The  length  of  the  device  chosen 
for  the  GaAs,  LoT,  and  Hi  t  calculations  was  10  pm 
The  silicon  parameters,  which  were  different  from  those 
of  the  other  three  elements,  were  chosen  in  a  wa>  to  make 
the  environment  as  represented  by  the  semiconductor 
equations,  similar  for  electrons  in  gallium-arsenide  and 
silicon.  For  this  we  concentrated  on  the  equation  for  total 
current  density  for  the  carriers  in  the  FET: 

J  =  NcME  +  cl)  grad  V  + 1  ( | ) 

*  <17 

where  M  is  the  field  dependent  mobility  (\V\/\E\)  and  D 
is  the  field  dependent  dtffusivity.  N{X,  Y,T )  is  the  mobile 
carrier  density  and  E(X,Y,Y)  the  electric  field.  <  is  the 
permittivity.  We  normalize  |6]  (! ).  writing  it  as 


i 


electric  field  t(> 


Fig.  3  Silicon  parameters  used  in  the  calculation.  The  parameters  are 
adapted  from  [5], 


Fig.  4.  Schematic  representation  of  a  Schottky  gate  field-effect 
transistor. 


For  the  silicon  calculations  the  velocity  and  diffusion  field 
curves  [5]  are  shown  in  Fig.  3. 

The  FET  device  circuit  configuration  is  shown  in  Fig.  4 
and  the  bulk  of  our  calculations  are  for  /V0—  10l5/cm3, 
with  the  exception  of  silicon,  for  which  the  calculations 
were  done  at  5x  IOl5/cm3.  The  reason  for  this  difference 
is  discussed  below.  With  regard  to  Fig.  4  we  note  that 
placement  of  the  source  and  drain  regions.  In  this  config¬ 
uration  we  concentrate  on  the  physics  associated  with  the 
space  charge  and  avoid  effects  due  to  the  geometrical 


where  f=  T/t,  (=X/ Z-0.  T  =  E/Ep,  p  =  M  / Mu.  /.0=  Vpl. 
n  =  .V/  N0,  and 

Ld  =  (Dt),/2  (3) 

which  is  a  field  dependent  Debye  length  .  r  is  the 
dielectric  relaxation  time.  At  low  values  of  electric  field 
and  for  a  doping  level  of  10'Vcm3.  LD  is  approximately 
equal  to  0.14  pm  for  GaAs.  Also,  in  (2).  Vp  is  the  peak 
electron  velocity,  which  for  GaAs  is  approximately  2.2  x 
107  cm/s  at  10l5/cm3,  M0  is  the  low  field  mobility,  and  Ep 
is  the  field  at  Vp.  In  our  simulations  we  chose  the  ratio 
Z-p/Eo  to  be  the  same  for  GaAs  and  silicon.  For  this  case, 
at  least  at  low  fields,  the  carriers  would  find  the  different 
semiconductor  equations  indistinguishable.  Differences 
will  occur  at  high  bias  levels  where  the  effects  of  NDM 
and  velocity  saturation  come  into  play. 

If  we  use  the  Einstein  relation  £>=  M0k0T0/e  then  the 
ratio  (Lp/  L0)2  is 


(4) 


and  for  GaAs  at  A/0=  I015,  Ep  =  3.25  x  I03  V/cm.  silicon 
with  A0=5x  10''/cnv\  and  Ep  =  7.14x  |07  V/cm  we  find 

/  Lp  \2  <  _  f  0.95  x  1 08/ V2-  cm  for  GaAs 

\L0j  k0T0  \  o.98x  108/V2  cm  for  Si 

We  use  these  parameters  in  the  simulation. 

The  normalization  associated  with  (2)  is  an  intrinsic 
normalization.  But  the  equation  for  current  flow  is  also 
subject  to  boundary  condition  at  the  source  and  the  dram 
contact.  To  give  the  carriers  the  impression  that  the  nor¬ 
malized  sample  lengths  are  the  same  it  is  necessary  to 
require  that  the  ratio  Lso/  be  the  same  for  the  silicon 
and  gallium-arsenide  device.  ( LSD  is  the  source-to-drain 
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Fig-  5  Tunc  evolution  of  the  current  and  potential  at  the  three  contact 
for  the  GaAs  FET.  The  normalization  scheme  is  indicated  in  the  test. 
Time  is  in  multiples  of  the  dielectric  relaxation  time 


separation.)  The  ratio  Lsn/ La  is  equal  to 


1-so  _  Lsn  _  1-spN o  /  e  i 
Lo  V,r  Ep  U 


(6) 


and  in  the  simulation  we  have  chosen  LSD  =  4.4  /cm  for 
silicon. 


III.  Illustrative  Calculations 

We  begin  the  discussion  with  the  gallium-arsenide  ele¬ 
ment.  The  output  of  the  simulation  is  illustrated  in  the 
next  few  figures.  In  Tig.  5  we  display  the  evolution  of  the 
current  at  the  three  contacts  and  the  potential  at  the  gate 
and  drain  contacts  for  the  situation  where,  at  first  the 
drain  bias  is  increased  to  0.5  tp  V(tp  =  EpLs[i)  and  the 
gate  bias  .  s  zero.  For  the  second  stage  the  gate  bias  is 
decreased  relative  to  ground  to  the  value  -0.6  vp.  The 
bias  changes  at  a  finite  rate  to  the  values  indicated  in  the 
figures.  The  sign  conventions  are  as  follows:  positive 
source  current  denotes  current  flow  into  the  device 
through  the  source  contact:  positive  gate  current  denotes 
current  flow  out  of  the  device  through  the  gate  contact; 
positive  drain  current  denotes  current  flows  out  of  the 
device  through  the  drain  contact.  Thus  from  Fig.  5  an 
increase  in  drain  bias  results  in  an  increase  in  drain  and 
source  current,  with  more  current  leaving  the  device 
through  the  drain  current  than  entering  through  the 
source  contact.  There  is  a  large  displacement  current  with 
carriers  building  up  on  the  gate  contact.  After  the  initial 
transient  there  is  a  significant  conduction  current  contri¬ 
bution  associated  with  the  physical  movement  of  the 
gate-depletion  region  until  the  steady  state  is  reached.  The 
subsequent  decrease  in  gate  bias  from  0  to  -0.6  vp  results 
in  an  increase  in  potential  on  the  drain  contact.  There  is  a 
further  buildup  of  charge  on  the  gate  contact  and  a 


movement  of  the  depletion  layer  to  the  bottom  of  the 
device.  This  effectively  cuts  off  the  current  flow  out  of  the 
drain  contact.  I  he  switching  times  are  estimated  from  Tig. 
5.  and  in  going  to  the  off-state  are  equal  to  the  magnitude 
of  the  difference  between  the  time  of  initial  relaxation  and 
the  time  at  which  the  drain  current  reaches  a  negligible 
value. 

The  internal  distribution  of  charge  and  current 
associated  with  the  current  and  potential  levels  of  Tig  5  is 
shown  in  Fig  6.  In  Fig  6(a)  we  display  a  set  of  current 
density  streamlines  through  the  device.  The  length  of  each 
streamline  is  proportional  to  the  magnitude  of  the  vector 
current  density  at  the  point  in  question.  The  maximum 
length  of  the  individual  v  and  y  component  before  overlap 
is  Jp  =  N0el/p.  For  the  stationary  stale  at 
60  r  s  in  Fig.  5  we  see  the  current  density  to  be  greatest 
under  the  gate  contact  as  required  by  current  continuity. 
In  Fig.  6(b)  we  show  the  current  density  for  the  situation 
where  the  dram  current  has  been  reduced  to  zero  by 
virtue  of  the  gate  bias  being  set  equal  to  -0.6  r  For 
these  GaAs  calculations  tp  =  3.2  V  The  density  of  charged 
particles  within  the  device  is  generally  nonuniform  and 
Figs.  6(c)  and  (d)  are  projections  of  this  distribution.  We 
point  out  that  the  particle  density  increases  in  the  down¬ 
ward  direction,  and  we  see  the  presence  of  the  depletion 
layer.  The  position  of  the  gate  for  these  calculations  are 
indicated. 

Quantitative  estimates  of  the  particle  density  and  poten¬ 
tial  for  the  steady-state  case  at  60r  s  are  illustrated  in  Fig. 
7.  where  we  display  contour  plots  of  these  quantities.  We 
note  that  most  of  the  potential  drop  is  under  the  gate 
contact.  It  is  clear  that  the  potential  difference  between 
the  contours  0.105  vp  and  0.210  r  are  insufficient  for  the 
field  to  reach  the  threshold  for  NI)M 
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Ftg  6  The  internal  distribution  of  current  and  charge  for  the  parame¬ 
ters  of  Fig.  5.  The  time  is  keyed  to  the  computations  of  Fig  5 
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Fig.  7  Contours  of  charge  density  and  potential  for  the  60r  s  Calcula¬ 
tion  of  Fig.  5 


Fig.  8.  Internal  distribution  of  charge  and  current  density  for  the 
situation  where  a  dipole  layer  forms  under  the  gate  contact  Here 
m  1  On,  and  oc —0.0. 

An  increase  in  drain  bias  to  vf  results  in  an  increase  in 
potential  on  the  drain  contact  and  a  moderate  increase  in 
drain  current.  We  are  into  current  saturation.  A  dipole 
forms  and  grows -under  the  gate  contact,  but  there  is  no 
current  instability.  For  a  lO-jrm  long  device  there  is  about 
2  pm  of  the  device  within  the  NDM  region.  Fig.  8 


illustrates  with  the  carrier  density  projection  and 
current-density  vector  display.  Fig.  9  illustrates  with  con¬ 
tour  plots  of  potential  and  earner  density.  We  have  also 
estimated  the  average  field  between  the  contours  of  con¬ 
stant  potential  and  find  that  at  the  drain  side  of  the  gate 
contact  the  field  exceeds  the  threshold  field  for  negative 
differential  mobility.  The  values  are  indicated  in  the  fig¬ 
ure. 

As  indicated  in  the  introduction,  calculations  with 
matenals  with  high  field  saturated  dnft  velocities  also 
show  dipole  formation.  Fig.  10  illustrates  this  for  silicon, 
where  we  show  the  contour  plots  of  earner  density  and 
potential.  The  average  field  under  the  gate  is  seen  to 
exceed  the  cntical  field  of  7  kV/cm.  We  also  show  the 
vector  current  density  display.  Here  the  lines  begin  to 
overlap  when  the  current  density  exceeds  N0eyp,  with 
equal  to  1 .0 x  107  cm/ s. 

IV.  Switching  Cai  citations 

The  above  discussion  presented  the  distribution  of 
charge,  potential,  and  current  within  a  device.  We  next 
compute  the  switching  properties  of  the  device.  We  begin 
by  presenting  the  current  voltage  charactenstics  under 
zero  gate  bias  conditions  for  the  four  simulated  elements. 
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Fig.  9,  Carrier  density  and  potential  roniours  fir  the  calculation  of 
Fig.  8.  Also  shown  are  estimate  of  the  field  between  the  contours. 


Fig.  10.  Contours  of  carrier  density,  electric  field,  and  current  vector 
for  a  three-terminal  silicon  element  showing  the  presence  of  a  dipole 
layer  under  the  gate.  Here  og  —  3.01),,,  uc«0.0.  The  peak  field  at  the 
bottom  of  the  channel  occurs  between  the  7.89 Ef  contours. 


This  is  displayed  in  Fig.  II.  The  current-voltage  relation 
for  each  element  is  normalized  to  lp  and  vp,  where  lp  = 
N0eVpA  and  vp~  EpL.  (A  is  the  cross-sectional  area  in  the 
Y-Z  plane.)  We  note  that  the  simulated  points  for  GaAs 
are  bracketed  by  those  of  Hi  V  and  Lo  V,  while  those  for 
silicon  are  approximately  equal  to  that  for  Hi  V. 

In  our  calculations  switching  takes  place  between  u 
point  lD.  i p  at  zero  gate  bias  and  a  second  point  at  a  level 
>  f  zero  drain  current.  Both  points  lie  on  a  dc  load  line.  In 
the  switching  process  an  amount  of  energy 

m =  fT‘vD(T)Ip(T)dT  (7) 

JT, 

•  unsferred  We  estimate  fifT  by  the  power-delay  prod- 

,  .  i  r  i /„(  i  ,m  Tj  r,) 

'  •» r  -esuits  in  Fig  12  show  the  power-delay  product 
‘  n  ru-menis  with  a  load  /!«  R0%  where  R0  is  the 
.  ,  .  - r  of  r.iih  of  the  elements.  We  also  show 
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ig.  II.  On-stale  current-voltage  relation  for  the  four  elements.  Note 
vp  —  3  2  V  for  GaAs.  Lo  y,  and  HI  V.  vp  m  3. 1  V  for  silicon. 
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Fig.  12.  Normalized  power-delay  product  and  normalized  switching 
times  versus  normalized  drain  voltage  for  the  four  elements  used  in  the 
calculation. 


the  switching  times  normalized  to  the  dielectric  relaxation 
time.  Perhaps  the  most  significant  feature  of  this  calcula¬ 
tion  is  that  GaAs,  Lo  V,  Hi  V,  and  the  scaled  silicon  device 
all  appear  to  have  similar  normalized  characteristics  as  far 
as  (7)  is  concerned.  The  differences  appear  in  the 
normalized  switching  times,  with  the  HiF  being  the 
fastest. 

In  examin'  these  switching  results  two  extreme  situa¬ 
tions  should  be  considered:  low-  and  high-bias  levels.  At 
the  high-bias  end  the  GaAs  and  Lo  V  element  show  similar 
behavior,  with  both  elements  displaying  longer  delay 
times  than  that  of  Hi  V.  This  result  provides  strong  evi¬ 
dence  that  saturation  in  the  drift  velocity  of  GaAs  is  the 
principle  determinant  of  its  high-bias  switching  speed. 
With  regard  to  the  scaled  silicon  device,  its  switching 
times  and  power  delay  product  are  similar  to  that  of  Hi  V, 
an  element  whose  mobility  is  better  than  four  times  hig¬ 
her,  and  whose  velocity  at  saturation  is  2.2  times  higher. 
In  our  calculations  the  high-bias  switching  times  appear  to 
be  determined  by  the  ratio  L /  VMV  where  L  represents  the 
length  of  the  high  field  region.  The  silicon  device  has  a 
gate-to-drain  spacing  that  is  approximately  1/2  that  of 
Hi  V ,  which  may  account  for  the  similar  switching  times. 

At  low-bias  levels  we  may  expect  the  switching  proper¬ 
ties  to  be  adequately  represented  by  the  Shockley  analysis 
[8]  in  which  case  the  switching  times  are  proportional  to 
[8],  [9]  r  ( L/af  where  a  is  the  channel  height  and  r  is  the 
dielectric  relaxation  time.  This  result  teaches  that  materi¬ 
als  with  different  mobilities  can  yield  similar  delay  times, 
provided  their  “aspect”  ratios  and  their  dielectric  relaxa¬ 
tion  times  are  similar.  The  scaling  associated  with  (l)-(6) 
resulted  in  comparable  values  of  r  and  L/a  for  all  four 
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nonlinear  elements.  Thus  the  low-bias  switching  speeds  of 
all  elements  should  be  similar. 

Generally,  the  net  switching  time  of  any  bias  level  is 
composed  of  several  contributions  including  transit  and 
dielectric  relaxation  times,  parasitic  and  dissipative  load¬ 
time  constants.  It  would  appear  that  if  the  high-bias 
switching  times  are  determined  by  velocity  saturation  and 
transit  through  high  field  regions,  than  so  long  as  t  is  less 
than  the  transit  time,  the  value  of  the  doping  level  will  be 
important  only  insofar  as  it  influences  the  field  distribu¬ 
tion  and  breakdown  conditions.  On  the  other  hand,  at 
low-bias  levels  the  switching  times  will  be  more  sensitive 
to  doping  as  indicated  in  the  above  paragraph. 

The  above  discussion  clearly  suggests  that  at  high-bias 
levels  similarly  designed  silicon  and  gallium  FET  switches 
should  yield  comparable  delay  times.  There  are,  however, 
limits  to  which  this  result  is  valid.  Continual  reduction  in 
the  active  region  length  to  achieve  shorter  switching  times, 
will  begin  to  lead  to  nonequilibrium  transient  behavior. 
These  transient  effects  lead  to  higher  transient  velocities 
in  gallium-arsenide  than  in  silicon,  which  should  result  in 
shorter  switching  times  for  the  former.  At  low-bias  levels, 
but  for  similar  doping,  silicon  should  yield  a  delay  time 
longer  than  that  of  GaAs  by  the  ratio  of  their  mobilities. 
At  intermediate  bias  levels  GaAs  should  generally  be 
faster  than  silicon. 

With  regard  to  the  power  -delay  product  in  Fig.  12,  it  is 
normalized  to  the  quantity 

Now  for  gallium-arsenide  with  L,  =  10  gm  and  l^  =  1.21 
gm 

/plpr=  1.28<L.  (GaAs).  (9) 

For  the  scaling  scheme  we  have  used  for  silicon,  with 
Lz  =  4.4  gm  and  =0.536  gm.  the  normalized  products 
are  approximately  the  same: 

lpvpT=\.20  tL.  (silicon).  (10) 

Thus  for  a  device  whose  width  is  L.  =  10  gm.  there  is  a 
power-delay  product,  at  a  normalized  drain  potential  of 
vp  that  for  both  is  approximately  equal  to  0.012  pj  as 
indicated  in  Fig.  12.  We  note,  however  that  for  a  silicon 
and  GaAs  switch  of  similar  doping  and  dimensions,  low- 
bias  level  switching  should  lead  to  higher  power-delay 
products  for  silicon.  At  bias  levels  well  into  saturation 
they  should  be  more  nearly  similar. 

A  recent  compilation  of  data  by  Greiling  [10]  on  simi¬ 
larly  designed  GaAs  and  Si  switches  shows  silicon  to  have 
(low-bias)  delay  times  six  times  longer  than  that  for  GaAs. 
At  high-bias  levels  the  switching  times  of  both  reflected 
the  fact  that  the  high  field  carriers  were  into  velocity 
saturation.  Further,  at  the  high-bias  levels  GaAs  was  only 
within  a  factor  of  two  faster  than  silicon. 

V.  Conclusions 

The  purpose  of  this  study  was  to  examine  some  of  the 
switching  properties  of  GaAs  FETs  and  to  compare  them 
to  silicon  FET’s.  The  results  demonstrate  that  at  high-bias 
levels  the  delay  time  of  GaAs  is  determined  by  velocity 
saturation.  Further  so  long  as  nonequilibnum  effects  may 
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be  ignored,  switching  speeds  for  GaAs  and  silicon  will  be 
similar  for  comparably  designed  structures  operated  at 
high-bias  levels.  (This  latter  conclusion  ignores  differences 
in  parasitic  contributions.)  At  low-bias  levels,  scaling 
principles  indicate  the  silicon  delay  times  will  be  longer 
than  those  of  gallium-arsenide  by  the  ratio  of  their  mobili¬ 
ties.  Intermediate  fields  will  result  in  shorter  delay  times 
for  gallium-arsenide. 

Rhirinos 

|l|  H  l.  Grubin.  D.  K  Kerry,  and  K  R  Gleason,  Solid-State  Elec¬ 
tron..  vol.  23.  p  157.  1980. 


[2)  D.  P.  Kennedy  and  R  R  O’Brien,  IBM  J  Ha  Develop  ,  vol  14, 
p  95.  1970 

[3)  H.  L.  Grubin  in  Proc  Sixth  Biennial  Cornell  Elect  Eng.  ConJ.. 
(Cornell  University,  llhaca,  NY),  p.  399,  1977. 

[4)  H.  L.  Grubin  and  T.  M.  McHugh,  Solid-Stale  Electron.,  vol.  21, 
p.  69,  1978 

[5)  C.  Jacaboiu,  C.  Canaii.  G.  Ottaviam,  and  A.  Albergigi  Quaranta, 
Solid-State  Electron.,  vol.  20,  p  77,  1977. 

[6]  B.  W.  Knight  and  G  A  Peterson.  Phys.  Rev.,  vol.  155.  p.  393, 
1967. 

[7]  P.  N.  Butcher,  Rep  Progr  Phys.,  vol.  30,  p  97,  1967. 

[8]  W.  Shockley.  Proc  IRE,  vol.  40,  p  1365,  1952 

[9]  G.  C.  Dacey  and  I.  M  Ross,  Bell  Sy.it  Tech  J  ,  vol.  34,  p.  1149, 
1955. 

(10)  P.  Greiling,  pnvale  communication. 


St>iui  Slait  Cfa-fANUt  s  Vol  2J.  pp  !'M72 
©  Pef|«mon  Press  1  id  IHM  Pnnied  in  (<«jl  BnUtfl 


00  M-  I I0I/KM0I0I  -01 57/S02  «V0 


SPONTANEOUS  OSCILLATIONS  IN  GALLIUM  ARSENIDE 
FIELD  EFFECT  TRANSISTORS 

H.  L.  Grubin 

United  Technologies  Research  Center.  East  Hartford.  CT  06108,  U  S  A 
D.  K.  Ferry 

Colorado  State  University.  Fort  Collins,  CO  80521.  U  S  A 
and 

K.  R.  GLEASONt 

Naval  Research  Laboratory,  Washington.  DC  20375.  U  S  A 

(Received  19  February  1979;  in  revised  form  19  May  19791 

Abstract— We  present  results  of  experiments  and  numerical  simulations  designed  lo  reveal  the  presence  of 
spontaneous  oscillations  arising  from  negative  differential  mobility  effects  in  gallium  arsenide  field  effect  tran¬ 
sistors  The  measurements  include  d  c.  and  pulsed  current/voltage  vs  temperature  characterization,  sampling  scope 
measurements,  spectral  analysis  to  40  GHz  and  observation  of  light  emission  The  simulation  is  a  time  dependent 
large  signal  transient  analysis  arising  from  a  fully  two-dimensional  solution  of  the  self-consistent  potential  and 
charge  within  the  device 


1  INTRODUCTION 

The  gallium  arsenide  field  effect  transistor  has  been 
around  a  long  time  [1.21  showing  promise  of  delivering 
high  power  while  sustaining  low  noise  levels.  GaAs  is. 
however,  subject  to  hot  electron  effects  that  manifest 
themselves  in  the  presence  of  highly  nonuniform  charge 
layers  or  "domains".  And  depending  on  device  size  and 
bias  values,  these  domains  are  either  trapped  near  the 
gate  or  drain  contacts(3],  or  propagate  between  them. 

The  presence  of  trapped  high  field  domains  is  expected 
to  profoundly  affect  FET  performance  and  is  currently 
subject  to  much  study [4-6],  However,  with  the  excep¬ 
tion  of  early  empirical  studies [2, 7)  and  recent  numerical 
studies[3. 8.  10).  the  dynamic  properties  of  the  propagat¬ 
ing  domain,  the  conditions  for  their  existence,  and  their 
influence  on  device  performance  has  until  recently ! II- 
14]  been  largely  ignored.  The  renewed  interest  arises 
from  a  need  to  determine  the  extent  to  which  the  occur¬ 
rence  of  spontaneous  oscillations  in  GaAs  FETsTl  1-14) 
are  a  consequence  of  negative  differential  mobility 
(NDM).  In  the  following  we  discuss  the  influence  of 
trapped  and  propagating  domains  on  FET  operation.  We 
also  present  evidence  for  the  highly  probable  conclusion 
that  the  propagating  domain  is  the  origin  of  our  obser¬ 
vation  of  spontaneous  oscillations  in  GaAs  FETs. 

The  study  reported  here  is  in  two  parts:  (i)  numerical 
simulation  and  (ii)  experiment.  The  numerical  studies 
reveal  the  presence  of  two  qualitatively  different  trapped 
domains  and  one  traveling  domain.  The  first  trapped 
domain  forms  within  the  conducting  channel  and  at  the 
drain  edge  of  the  gate  contact.  (4-35).  Its  formation  is  a 
direct  consequence  of  velocity  limitation  in  GaAs  and  its 
presence  is  responsible  for  current  saturation  at  values 
of  drain  voltage  significantly  below  that  predicted  by 
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Shor<:/ey(15].  The  second  trapped  domain  is  qualita¬ 
tively  different  and  forms  in  FETs  capable  of 
sustaining  a  current  instability  arising  from  (NDM).  This 
domain  forms  near  the  drain  contact,  at  high  drain  and 
low  gate  bias  levels  and  is  coincident  with  the  cessation 
of  oscillatory  activity  in  GaAs  FETs.  (The  extent  to 
which  it  is  removed  from  the  gate  region  depends  on  the 
gate  to  drain  spacing  and  the  doping  concentration  of  the 
epitaxial  region.)  The  domain  is  qualitatively  similar  to 
the  high  anode  field  configuration  in  two  terminal 
devices,  where  it  is  also  responsible,  in  many  cases,  for 
the  cessation  of  oscillation!  16).  Each  of  the  low  and  high 
drain  bias  trapped  domains  may  pulsate  in  time  and  so 
enhance  the  noise  existent  in  the  FET.  In  the  case  of  the 
traveling  domain,  it.  as  for  some  transferred  electron 
logic  devices  (TELD)  is  launched  at  the  drain  aide  of  the 
gate  contact,  propagates  to  the  drain  contact  and  then 
recycles  successively.  While  the  traveling  domain  need 
not  be  "stable"  in  the  sense  'of  two  terminal  NDM 
devices[16).  its  threshold  and  quenching  conditions  are 
similar  to  that  for  two  terminal  devices  [3, 8). 

The  experimental  studies  reported  below  consist  of 
pulsed  and  d.c.  measurements  performed  as  a  function  of 
temperature  and  under  different  circuit  conditions.  We 
have  determined  the  spectral  content  of  the  oscillations 
when  they  occurred,  and  for  pulsed  measurements  have 
sampled  the  current  and  voltage  at  different  points  within 
the  pulse.  The  devices  studied  were  empirically 
separated  into  two  groups  as  determined  by  the  ratio 

K  = 

_ Gate  voltage  at  cutoff  ( »  Vco) _ 

Drain  voltage  at  the  onset  of  current  saturation  ( ■  V*)’ 

(I) 

Devices  with  K  greater  than  unity  sustained  current 


157 

III-F 


H  I  (.kihiw/jj / 


I'S 

ouillntumy .  those  with  k  approximately  equal  1 o  or  lets 
thiin  unity  Mere  eleelrietillx  quiet  I  or  those  cases  with  k 
greater  than  unitx  oscillations  began  bexond  ihc  knee  of 
current  saturation  In  each  case  the  oscillations  ceased  .it 
high  x. lilies  of  dram  and  gate  bias  In  most  of  the 
electric. illx  un-t.ible  ilex  ices  the  onset  of  the  oscillation 
is  marked  bx  a  drop  in  current  In  till  uno  the  mslabihlx 
.i/’/’iiiri  jo  be  imliiileil  alien  tin  aieraite  held  under  the 
Hate  eontaet  hr\t  reaehex  the  \I>M  threshold  held  i.ilie 
This  result  holds  true  under  different  ambient  tem¬ 
perature  conditions  An  estimate  ot  the  carrier  xelocitx  in 
the  channel  hetxxeen  the  fate  and  drain  contact'  in¬ 
dicates  that  /'nor  ro  the  o»r  illation  the  eleetron\  are 
U  slut  II  x  dntline  at  speedy  in  exet'w  or  the  saturation  drift 
teloeitx  When  the  oscillations  ceased  xxtth  increasing 
eate  bias  the  carrier  xelocitx  dropped  to  xalues 
sifniffcantlx  beloxx  the  saturation  drift  xelocitx  l-.ich  of 
the  aboxe  results  are  consistent  xxith  an  instahililx  arising 
from  the  presence  of  nefatixe  differential  mobilitx 

’  Dt-Xlft  Cl  \ss| HI  xrtllN 

The  classification  parameter  k  isee  eqn  li  emphasizes 
differences  hetxxeen  the  (la  As  and  classical  I  H.  xxhere 
for  the  latter 1 l<]  l..  and  \ ...  are  equal  For  (ia  As.  the 
xerx  earliest  experiments 1 1.  ’.  '|  shoxxed  k  x nixing  from 
unitx  to  xalues  significanllx  hif her  than  umtx  iCurrent 
dexices|l'|  haxe  not  altered  this  situation'  In  those 
experiment'll.?!  xxhen  K  -  1.  no  instabiltix  xx.ts  repor¬ 
ted  W  hen  k  xx.is  at  the  high  end  of  the  scale,  (he  studies 
in  Refs  (?|  and  (’I  showed  that  the  dexice  no  longer 
exhibited  the  classical  characteristics,  but  instead  sus¬ 
tained  spontaneous  oscillations  for  a  ran.ee  of  bias 
xalues 

The  results  of  the  computer  simulation  corroborate  the 
empirical  classification  But.  in  addition  relate  stabilitx  to 
the  material  parameter  r„  the  electron  saturated  drift 
xelocitx  Specifically,  the  simulations  shoxx  that  a  neces- 
s.irx  condition  for  a  current  instahililx  is  that  the  current  at 
saturation  either  exceed  or  be  xxithm  a  narrow  band  about 
the  xalue 

/.  =  t?i 

In  eqn  t?i.  is  the  open  channel  conductance 
i  =  Area  ( L  '  Resistix  itx  n.  E.  =  i .  u.  I  is  the  source  to 
drain  separation,  and  a  the  loxx  field  mobilitx  of  the 
semiconductor  In  txxo  terminal  dexice  lanfua.ee.  E.  is  the 
sustaining  field { 16].  Summarizing  If  IDDS  denotes  the 
zero  gate  bias  current  at  saturation  then 

IDSSs-  /,  ■  Stable current-xc'ltagecharacter'stics 
IDSS  >  I. :  Unstable  current-soilage  characteristics.  FM 

The  aboxe  conclusion  i»  essentiallx  similar  "o  that  ixf 
Vamagucht  et  et!  (8) 

The  results  discussed  aboxe  summarize  the  general 
features  of  the  GaAs  FFT  The  space  charge  distnbution 
xxithm  the  FFT  is  responsible  for  these  properties  and 
six  ill  be  discussed  in  the  next  section 
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The  folloxxing  analxsis  is  for  dcxiccs  in  the  contigura 
lion  of  lie  ll.il  xxhere  the  source  and  drain  contacts  are  at 
parallel  ends  ol  the  dexice  This  configuration  eliminates 
geometrical  effects  due  to  coplan.ir  placement  of  the 
source  and  dram  contacts  (  oplan.u  contact  calculations 
haxe  been  performed  for  the  configuration  shoxx n  in  Tie 
Itbi.  xx Inch  also  include  the  presence  of  a  substrate  Amt 
the  results  of  these  studies  slum  injection  into  the  sub 
strate.  xxith  ihe  instabilities  still  remaining  Ml  of  the 
i  ilculutions  reported  beloxx  are  foi  elements  xxith  doping 
lexels  nominally  equal  to  Iff1'  cn.‘.  a  xalue  approx  I  < 
to  ?  orders  of  magnitude  beloxx  that  of  presenllx  f.ihri 
cated  devices  This  choice  xxas  made  because  the  cost  of 
computing  the  contributions  of  transiting  domains  in 
creases  monotomcallx  xxith  increased  doping  lex  el 1 1 S | 
Hoxxexer.  a  number  of  representative  calculations  haxe 
been  performed  at  higher  doping  lexels  xxith  similar 
results  We  note  that  the  results  of  Yamaguchi  el  al .  |X| 
xxeie  at  )  '•  lff"'icm‘  and  xxhere  there  is  overlap  xxith  our 
studies  there  i'  agreement  Also,  the  main  results  of  the 
experiment  confirm  the  conclusion  developed  at  these 
loxxer  doping  lexels 

bl  Current  eoltaee  <  haraetenxlii  x  and  \puee  elnirue 
distribution  for  FFT  with  K  --  I 

The  space  charge  distributions  for  a  k  I  dexice  are 
displaced  in  Figs  ?-4  The  inset  of  each  figure  shoxxs  its 
current-collage  characteristic  xxith  the  loxxer  case  letters 
kexed  to  the  space  charge  distribution  The  current  is 
normalized  to  the  xalue 

I,  =  (i„F.  I  1 4 r 

xxhere  Er  o  the  electric  field  it  peak  xelocitx  And  the 
xoltage  is  normalized  to 

l  ,  =  E,  I  (5 1 

The  bold  curve  is  the  electron  drift  xelocitx  electric  field 
relation  scaled  to  current  and  xoltage  The  closed  circles 
denote  computed  point'  and  each  point  is  a  stable  time 
independent  point 

For  the  space  charge  distribution  the  t  axis  represent' 
the  longitudinal  spatial  dependence  The  vertical  charge 
distribution  i'  along  the  v-axis  The  source  contact  is  in 
the  x-:  plane  at  r  =  0.  the  drain  contact  is  at  t  =  L  The 
gate  contact  is  in  the  plane  at  the  indicated  points 

At  loxx  values  of  drain  bias  and  zero  gate  bias  (Fig  ?al. 
there  is  a  region  of  charge  depletion  under  the  gate 
contact  The  charge  region  spreads  somewhat  beyond  the 
gate  contact  boundaries  (enclosed  rectangular  parallele¬ 
piped*.  but  pretty  much  the  space  charge  distribution  is 
what  one  would  expect  from  classical  arguments]  15)  An 
increase  in  drain  xoltage  to  the  knee  of  the  current 
xoltage  relation  iFig  ?bi  leads  to  a  domain,  ^presenting 
the  first  significant  departure  from  the  classical  analysis 
The  domain  consists  of  an  accumulation  layer,  which 
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I  ig  I  IX’sice  circuit  configuration  used  in  the  FFT  simulations,  (al  Parallel  source  and  drain  contacts  (hi  C'oplanar 
sourci  and  drain  contacts  Substrate  is  included 


forms  within  the  boundaries  of  the  gate  region  and  is 
followed  downstream  by  a  depletion  layer.  Domain  for¬ 
mation  is  a  consequence  of  current  continuity  and  velo¬ 
city  limitation  and  occurs  w  hen  the  current  density  under 
the  gate  exceeds 

Jr  -  Sn(lr  (6) 

wheret,,  is  the  maximum  carrier  velocity.  Further  in¬ 
creases  in  drain  bias  result  in  an  increase  in  the  amount 

of  charge  accumulated  and  a  broadening  of  the  domain 
(Fig.  2c I.  We  point  out  that  while  the  presence  of  a 
domain  is  necessary  for  a  current  instability,  it  is  not 
sufficient.  By  analogy  with  two  terminal  devices  the 
current  density  downstream  from  the  dipole  must  exceed 
the  sustaining  current,  J„  where 

./,  =  /Voff.  <71 

(see  also  eqn  2).  This  criteria  is  approximately  satisfied 
for  the  FET  hut  is  not  as  stringent.  It  suffices  to  state 
that  the  current  density  for  the  K  =  I  devices  is  generally 
too  low  to  sustain  an  instability. 

The  calculations  of  Fig.  .1  are  for  finite  values  of  gate 
bias.  At  low  drain  bias  levels  the  depletion  layer  extends 
well  beyond  the  gate  boundaries.  We  see  a  region  of 
local  charge  accumulation,  although  clearly  there  is  a  net 
depletion  of  carriers.  However,  increasing  the  drain  bias 
results  in  a  strong  region  of  charge  accumulation,  fol¬ 


lowed  by  significant  depletion  until  the  drain  contact  is 
reached  where  the  boundary  conditions  require  a 
significant  increase  in  the  charge  density. 

Two  extreme  conditions  are  displayed  in  Fig.  4  Here 
at  low  values  of  drain  bias  and  a  value  of  gate  bias  near 
cutoff  we  see  the  presence  of  a  depletion  layer  broadly 
surrounding  the  gate  region.  At  zero  gate  bias  levels  but 
very  high  drain  bias  levels  we  see  the  presence  of  a 
region  of  charge  accumulation  extending  from  the  gate 
region  to  the  drain  contact. 

We  summarize  the  above  results  hy  stating  that  for  a 
collection  of  ostensibly  classical  current-voltage  curves 
the  space  charge  distribution  is  far  different  in  many 
circumstances  from  that  developed  originally  for  the 
FET(|5],  Thus  while  the  current-voltage  relation  does 
not  provide  us  with  detailed  information  about  the  dis¬ 
tribution  of  space  charge  for  the  K  =*  1  device,  we  expect 
that  the  form  of  the  distribution  will  affect  the  equivalent 
circuit  element  representations  of  the  device  [4.51. 

?(c).  Current  voltage  characteristics  and  space  charge 
distributions  for  FETs  with  K  ~>  \ 

For  a  given  doping  level  and  source  drain  spacing,  the 
principle  difference  between  devices  with  K  >  I  and 
K  *  1  is  that  the  former  have  wider  channels,  therefore, 
lower  open  channel  resistances,  and  for  a  given  bias 
draw  more  current.  For  these  devices  the  current  density 
is  usually  high  enough  to  sustain  an  instability.  In  our 
simulations  the  onset  of  the  instability  marked  the  onset 
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Fig.  2.  Space  charge  configuration  as  a  function  of  bias  Profiles  correspond  to  indicated  points  on  the  current  voltage 
characteristics.  Solid  line  represents  GaAs  velocity  vs  electric  field  relation  scaled  to  current  and  voltage  Device 
configuration  is  also  shown,  with  dimensions  in  microns 


of  current  saturation.  In  our  experiments,  the  instability 
occurs  after  the  device  enters  saturation.  The  simula¬ 
tions  for  the  wider  channel  device  are  displayed  in  Figs. 

S  and  6. 

Figure  5(a)  shows  the  charge  distribution  at  a  bias  level 
just  below  that  necessary  to  launch  a  traveling  domain. 
Within  the  gate  boundaries  there  is  charge  depletion.  But 
towards  the  bottom  of  the  channel  and  just  before  the 
boundary  there  is  a  trace  of  an  accumulation  layer- 
implying  that  the  current  density  is  highest  there.  Coup¬ 
ling  this  result  tc  the  current  value  in  the  inset  to  Fig.  5, 
we  have  the  result  that  within  the  gate  boundary  the 


current  density  is  approximately  equal  to  Jr.  while  within 
the  gate-to-drain  region  the  current  density  exceeds  J,. 

The  above  conditions  are  sufficient  to  launch  a  domain 
in  two-terminal  devices,  and  indeed  here  a  slight  increase 
in  drain  bias  results  in  the  launching  of  a  domain  (Fig. 
5b).  In  this  case  further  increases  in  bias  result  in  the 
charge  accumulation  layer  spreading  toward  the  bottom 
of  the  channel  (within  the  gate  boundaries)  and  also 
toward  the  drain  contact.  The  domain  detaches  from  the 
gate  region  and  periodic  transit  time  oscillations  occur. 
Figure  5(b)  shows  a  "snapshot"  of  the  domain  at  one 
instant  of  time  during  propagation  (Recycling  domains 
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are  illustrated  in  Fig.  5.  Ref.[J|.i  The  propagating  and 
recycling  domain  can  persist  for  a  rather  wide  range  of 
drain  bias  values.  The  amplitude  of  these  resulting  transit 
time  oscillations  are  bias  dependent  as  illustrated  in  Fig. 
7. 

Figure  7  contains  plots  of  drain  and  source  current  vs 
drain  potential  at  different  times,  with  time  eliminated 
between  these  In  Fig,  7(ai  we  plot  drain  current  vs  drain 
voltage  for  four  different  values  of  dram  bias.  The  bold 
lines  denote  calculated  results,  the  dashed  lines  are  an 
estimate  of  the  current  \s  voltage:  the  skew  lines  are  the 
drain  d  c  load  lines  In  each  calculation  displacement 
current  contributions  result  in  transient  drain  current 
values  in  excess  of  their  steady  state  values.  This  is 
apparent  for  the  bias  levels  of  14  and  1.6  V„,  where  after 
steady  state  has  been  reached  the  displacement  current 


goes  to  zero  and  the  system  relaxes  to  the  value  denoted 
by  the  closed  circle.  When  the  domain  is  propagating  the 
current  oscillates  along  the  load  line.  At  a  bias  of  1.6  V',, 
the  current  oscillation  amplitude  is  somewhere  around 
20W  of  its  d.c.  value.  Further  increases  in  drain  bias 
result  in  an  increase  in  the  amplitude  of  the  oscillation 
The  average  current  associated  with  these  oscillations  is 
denoted  by  the  ,x's.  We  note  that  the  average-current- 
voltage  relation  for  the  range  in  which  there  are  transit 
time  oscillations .  exhibits  a  region  of  negative  differential 
conductivity  (NDC).  The  NDC  is.  for  these  calculations 
dynamic  in  origin.  The  second  part  of  Fig.  7  shows  the 
source  current  vs  drain  voltage  for  the  same  calculation 
as  that  of  7(a).  Here,  in  Fig.  7(b)  we  see  the  presence  of 
superlinear  current-voltage  relations  at  low  drain  bias 
levels,  followed  by  looping  when  domains  are  present. 
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Both  of  these  effects  occur  because  of  large  transients 
through  the  gate  contact. 

With  respect  to  transit  time  oscillations  we  point  out 
that  no  experiment  has  yet  been  designed  that  pros  ides 
unequivocal  evidence  for  their  existence  in  commonly 
fabricated  GaAs  FF.T'  In  most  common  designs  the  gate 
to  drain  separation  is  somewhere  between  one-half  and 
five  microns,  and  for  these  values  it  is  not  dear  that 
conditions  exist  for  the  propagation  of  stable" 
domains! I6|.  The  result  is  that  the  properties  of  these 
propagating  domains  are  readily  influenced  by  the 
external  circuit  and  can  in  some  cases  be  controlled  bv  it 
An  example  of  this  complex  dipole-circuit  interaction 
exists  in  two  terminal  devices  (e  g  the  quenched  multiple 
dipole  mode)  16])  There  is  also  the  possibility  that  in 
place  of  the  negative  differential  mobility  induced  oscil¬ 


lations  in  GaAs  FETs  there  are  feedback  circuit  oscil¬ 
lations  that  are  independent  of  any  hot  electron  effects 
However,  in  the  experiments  reported  below  we  'how 
that  the  evidence  points  to  the  oscillations  as  being 
negative  differential  mobility  induced  Two  additional 
space  charge  profiles  are  important  in  this  regard  and  are 
discussed  below 

The  two  profiles  of  mitres!  arc  coincident  with  the 
suppression  of  transit  lime  oscillations.  In  one  case  a 
significant  increase  in  drain  bias  alone  suppresses  the 
instability  In  the  second  case  an  increase  in  the  mag¬ 
nitude  of  the  gate  bias  is  sufficient  to  suppress  it  The 
space  charge  profile  associated  with  the  high  dram  bias 
suppression  of  the  instability  is  displayed  in  Fig  6(a)  We 
see  that  an  accumulation  of  charge  forms  downstream 
from  the  gate  boundaries  and  extends  towards  the  drain 
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Hi;  '  \s  in  Fig  2.  hut  for  i  different  m/c  device  \  s  denote  transit  time  oscillations 


contact.  White  this  layer  generally  is  not  a  static  one  hut 
instead  pulsates  in  time  there  are  some  hroad  features 
associated  with  it  For  example.  Ihe  profile  forms  tit 
tallies  of  current  density  in  excess  of  1 ,  and  at  high 
tallies  i>f  drain  potential  Most  of  the  potential  drop  is 
across  the  gale-to-drain  region  l Ref  |3|.  Fig  71.  where  we 
may  expect  the  electric  fields  to  he  well  within  the 
saturated  drift  velocilv  region  Thus  the  carriers  within 
the  gate-to-dr.iin  region  are  traveling  at  their  saturated 
drift  velocity  values  This  fact  and  current  contimiitv 
with  the  current  densitv  exceeding  f.  are  responsible  for 
the  region  of  charge  accumulation  Bv  wav  of  analogs 
this  solution  is  qualitatively  similar  to  two  terminal 
device  solutions  which  show  the  presence  of  anode 
adjacent  domains|l6|  Further,  in  three  terminal  devices 
with  long  gate-to-drain  spaemgs  and  high  doping  con¬ 


centrations  ihe  anode  adjacent  domain  has  been 
experimentally  probed  I  lyl 

The  space  charge  profile  associated  with  the  high  gate 
bias  suppression  of  (he  instability  is  displaced  in  Fig. 
Mb)  This  shows  a  snapshot  of  the  space  charge  layer 
while  undergoing  a  damped  oscillation  This  laver  of 
charge  also  pulsates  m  time,  but  in  contrast  to  the  space 
charge  laver  of  Fig  Ma).  this  one  forme  at  current 
density  i  allies  that  tire  less  than  J.  Indeed  for  this  case 
most  of  the  potential  drop  falls  under  the  gate  contact 
and  extends  some  distance  downstream  from  the  edge  of 
the  gate  boundaries  (M  (Fig  Si  \s  long  as  the  space 
charge  laver  does  not  reach  the  drain  contact  the  electric 
field  downstream  from  the  edge  of  the  nonimiform  space 
charge  region  will  be  less  than  the  "two-terminal  sus¬ 
taining  field" |  lb)  Fe  i.  a  This  situation  is  analogous 


f  lu  fi  As  in  Kn>.  5.  hut  for  different  bia-  value' 


lo  the  Iwo  terminal  high  cathode  boundary  field 
profiles  (16). 

3(d).  Summary 

To  summarize,  there  are  several  key  features  asso¬ 
ciated  with  negative  differential  mobility  induced  in¬ 
stabilities  in  field  effect  transistors.  These  are: 

(I)  The  instability  occurs  in  devices  with  K  greater 
than  unity. 

(2t  The  instability  occurs  at  or  beyond  the  "knee"  of 
the  current-voltage  relation. 

(3)  A  necessary  condition  for  an  instability  is  that  the 
current  density  downstream  from  the  gate  boundary  be 
greater  than  or  approximately  equal  to  7,.  or 
equivalently,  that  the  drift  velocity  cither  exceeds  r, ;  or 
is  approximately  equal  to  r, 

(4)  Saturation  in  currert  occurs  when  the  average  field 
within  the  conducting  channel  under  the  gate  contact  is 


approximately  equal  to  the  threshold  field  for  negative 
differential  mobility. 

(5)  The  instability  can  be  suppressed  by  imposing  a 
gate  bias  sufficiently  high  to  reduce  the  current  density 
downstream  from  the  gate  boundary  to  a  value 
significantly  less  than  J,. 

(fi)  The  instability  can  be  suppressed  by  imposing  a 
sufficiently  high  value  of  drain  bias. 

I.ct  us  now  turn  our  attention  to  how  these  instabilities 
manifest  themselves  in  actual  devices. 

5  FXPFRIMFNTS  WITH  I.M.MI  M  AUSFNIDE 
FIFI  D  BTK  T  TRANSISTORS 

3(a).  Experiments  with  K  >  I 

Figures  R  and  9  display,  respectively,  temperature 
dependent,  d  c.  and  pulsed  drain  current  vs  drain  voltage 
characteristics  for  a  GaAs  FF.T  The  source,  drain  and 


Spontaneous  oscillations  in  gallium  arsenide  field  effect  transistors 


16' 


LOAD  LACS 


vOSAtN/  VP 


^3  18MCXONS 
LG-OMMOiONS 


!2  kV/cnt.  anJ  we  have  reduced  these  results  hy  approx 
1055  for  higher  fields  The  closed  circles  in  big  10(a)  are 
for  the  d.c  measurements,  the  t's  represent  pulsed 
measurements  1 1  n  sec  pulse  lengths.  Vi  dulv  eyelet  The 
pulsed  results  generally  surround  the  solid  line  and  arc 
higher  than  the  d  c  results  In  general  user  the  temperature 
range  exceeding  1^0  K  the  pulsed  measurements  yield  a 
yelocity  exceeding  i , 

( 4 )  The  average  held  under  the  gate  boundary  is  c-ti 
mated  from  the  relation 

/.  I  (4 1 

where  l’,.  is  the  longitudinal  potential  drop  under  the 
gate  contact  and  is  obtained  from  an  equation  that 
ignores  contact  resistance  contributions 
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big  ’  Instantaneous  current  vs  voltage  tor  the  device  of  lags  s 
and  6  tat  Drain  current  vs  drain  voltage  (hi  Source  current  vs 
drain  voltage 

gate  contacts  are  coplanar  with  the  source  and  drain 
contacts  separated  hy  approx  8.5  pm.  The  gate  is  more 
or  less  centrally  placed  with  a  length.  /....  equal  to 
3.0pm.  The  gate  width.  /  is  250  pm  The  epitaxial  layer 
thickness,  H.  was  3000  -  500°  A  with  a  nominal  doping 
level  of  I0l7/em\  The  mobility  measured  from  adjacent 
samples  varied  from  3000  to  4000  cm:/V-s.  The  cutoff 
voltage  for  this  device  was  84  V  This  device  sustained 
current  oscillations  as  detected  hs  a  spectrum  analyzer. 
All  six  criteria  listed  in  the  last  section  for  NDM  induced 
instabilities  were  satisfied  We  discuss  each  separately 
below.  Oscillation  characteristics 

( 1 )  K  =  V,  J  V'j,  =  18.9  +  0.8  V  1/2.5  =  3.9.  where  the 
cutoff  voltage  includes  as  assumed  built  in  potential  of 
approx.  0.8  V. 

(2)  At  each  value  of  ambient  temperature,  the  in¬ 
stability  occurred  beyond  the  knee  of  the  current  voltage 
characteristic. 

(3)  The  velocity  was  computed  from  the  relation 

t  IlK.cHl  I  (81 
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where  V is  the  drain  voltage  and  /  the  drain  current 
R„ I"  1 10..)  is  the  open  channel  resistance  and  may  he 
obtained  from  the  measured  drain  current  vol’age 
characteristics  of  the  device  via  the  relation 

R„  Rji  l’..  It  I  -  |l„,  •  V  ..1  V  :i/(  I  -  [l  -*  |(V„, 

♦  i..t  )  (ih 

where  we  have  adopted  the  viewpoint  of  eqns  (2l  and  <3j 
of  (irehene[22|  V...  in  eqn  Mil  includes  the  built  in 
potential  The  above  expression  for  R„  is  expected  to  he 
valid  at  low  values  of  drain  bias.  Assuming  a  built  in 
potential  of  0.8  V  and  reading  R,,  at  V ,  =  0  from  the  d.c. 
characteristics  we  compute  R„  The  values  are  listed  in 
Table  I.  where  we  haxe  ignored  the  temperature  varia¬ 
tions  in  l'*,  and  V',...  We  see  that  the  average  field  under 
the  gate  is  relatively  insensitive  to  temperature,  as  is  the 
threshold  field  for  negative  differential  mobility (20).  An 
average  field  of  4.2  kV/cm  is  consistent  with  a  mobility 
for  iC’/cm'  material  of  4200cmJ/V-s. 

(5 1  As  shown  in  Fig  10(b)  the  carrier  drift  velocity  at 
ihe  cessation  of  oscillations  is  below  t  „ 

(6)  All  oscillations  cease  at  high  drain  bias  levels. 

The  next  set  of  data  are  displayed  m  Fig.  1 1  where  we 
show  d.c.  and  sampling  scope  data  (taken  25  nsec  into  a 
50  nsec  pulse,  with  low  duly  cycle).  The  GaAs  FET  had  a 
source  to  drain  spacing  of  3.5  pm  and  a  gate  length  of 
1  pm  The  cutoff  voltage  for  this  device  was  3.0  V  and 
the  mobility  was  quoted  at  around  4500crrr/V-s.  With 
regard  to  the  first  two  oscillatory  criteria.  K  =  4  8  and  the 
instability  occurs  just  beyond  the  knee  of  the  current  vs 
voltage  characteristic.  The  peak  velocity  just  prior  to  fhe 
instability  at  ^7  ma  is  computed  from  the  equation 


And  the  values  prior  to  an  instability  are  displayed  as  a 
function  of  umbient  temperature  in  Fig,  10(a)  where  the 
error  bars  represent  uncertainties  in  the  value  of  H.  In 
Fig.  10(a)  the  solid  line  corresponds  to  .n  estimate  of  the 
saturated  drift  velocity  from  Ruch  and  Fawcett  |20|  and 
Freeman  and  Hobson |2 1 1.  This  data  is  available  only  to 


i  -  ulRol  (12) 

with  R«  ~  9  7  U  (from  pulsed  data)  and  p  =  4500cm J/V- 
s.  The  result  is  r  -068'  10’ cm/sec  which  falls  ap¬ 
proximately  within  the  range  of  values  of  Fig.  10  An 
,tv  ■’•■ace  threshold  field  of  4  kV/cm  for  a  one  micron  gate 
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Table  I  d  i  IXila 


T  CASE 

VDS  (vl 

10  X  10  3(a) 

Rd  (OHMS) 

Ro  (OHMS) 

Eav  (Kv/cm) 

383°K 

2.5 

82 

26.3 

23.0 

4  26 

299°K 

25 

90 

24.2 

21.' 

4.22 

166°K 

2.5 

105 

21.1 

18.4 

4  17 

77°K 

2  5 

115 

20.0 

17  5 

4.00 

Fig  8  feniperature  dependent  d  i  current  voltage  data  for  a  long  gate  IT- 1  IVmvc  r  t;.«metei-  oe  »n  text 


is  reached  at  a  drain  voltage  of  approx  O  K  V .  and  it  is 
here  that  we  see  substantial  preinstahilitv  current 
broadening.  The  oscillations  which  are  initiated  at  about 
1.2  V  are  not  suppressed  until  a  sufficiently  low  value  of 
current  density  is  attained  The  characteristic  of  the 
instability  is  suggestive  of  hot  electron  contributions 
richer  than  that  discussed  in  the  last  section  (see  also 
Ref  (I3|)  However,  consistent  with  (he  sixth  criteria  the 
oscillations  are  suppressed  at  high  values  of  dram  bias 
This  is  displayed  in  Fig  12  where  we  show  sampling 
scope  measurements  at'  increasing  drain  bias  for  the 
device  subjected  a  I  ^  jisec  pulse  We  see  (he  initiation 


of  oscillations  at  low  dram  bias  levels  with  subsequent 
suppression  at  high  drain  bias  levels 

4(c).  Experiment  »ilh  K  I 
The  final  set  of  measurements  is  for  a  device  with  fx 
approximately  equal  to  unitv  These  are  pulsed 
measurements.  The  vies  ice  has  a  cutoff  voltage  of  ap¬ 
prox  I  K  V  and  no  oscillations  were  detected  on  a  spec¬ 
trum  unulv/er  over  a  range  of  40 tilt/  The  device  had  a 
quoted  doping  concentration  of  *  ■  I0""cm  and  an  epi 
taxial  thickness  of  nm  Ihexe  values  were  not 


Spontaneous  oscillations  in  gallium  arsenide  held  effect  transistor- 


l  ie  s  I  emperature  rlependenl  pul'Cvl  rl.tl.i  f«*r  the  rles  ii.e  r’f  lag  *  I  ot  these  nic.isutetacnt- '  f 

soutse  from  ground 


consistent  with  pmchoff  sallies  using  the  equation 

V.  <A„//?i2f>.  (Hi 

With  a  quoted  salue  of  I  ViOqm.  eqn  (8)  yields  a  room 
temperature  range  of  salues  of  i.  associated  with  the 
product  V  ,H 

i  -I,'  •  If)'"  S„H.  P4i 

From  eqn  (I4l  with  \„ H  varying  from  IO':/cm‘  to  15  ■ 
10' ’/cm',  t  varies  from  about  0  5  ■  10' cm  sec  to  0  '  • 
10’ cm  sec  Ml  values  are  too  low  for  an  instability  to 
occur 

4IC)  /Jtu/mion 

The  experiments  discussed  above  provide  evidence 
that  the  GaAs  (  FT  t'  .apable  sif  sustaining  negative 
differential  mohilitv  induced  current  instabilities  These 
instabilities  occur  when  the  average  field  under  the  gate 
contact  exceeds  the  NDM  threshold  field  and  when  the 
current  densitv  at  the  onset  of  saturation  is  up- 
proximutelv  equal  to  or  exceeds  the  saturation  drift 
current  density  for  GaAs  There  are  several  additional 
experimental  aspects  that  should  he  highlighted  The  d  c 
measurements  of  fig  8  were  characterized  by  a  drop  in 
current  at  the  instahilitv  threshold:  the  pulsed  measure¬ 
ments  also  show  a  drop  in  current  The  d  c  measurements 
of  Fig.  1 1  do  not  show  a  dropbuck  in  current;  the  sampling 
scope  measurements  do  There  is  hvsteresis  in  the  d.c 
measurements  of  Fig.  8  with  the  instability  persisting  for 
increasing  and  decreasing  values  of  dram  bias.  Flowevcr 
on  the  decrease  the  instability  persisted  at  drain  bias 
levels  somewhat  below  the  instahilitv  threshold  This  is  a 


characteristic  of  two  tctmuiu!  tt.i'"i!  1  m 
With  regard  to  the  sampling  s.opc  m»  ■.:•  t  I 

1 1 .  we  see  Mime  spreading  of  the  curren'  v-  -a  ,  ■  p. -. .< 
to  the  instahilitv  Ibis  m.iv  .orre'p.  •  p  t 
pulsating  trapped  domain  nc.u  the  ca'i  r.fut  -  ,i  - 
cussed  in  the  earlier  sections 

The  oscillation  properties  reported  let  .  ,1. 
using  sampling  oscilloscope'  and  spcC'i.m.  .  mI./ci- 
Spectral  analvsis  »,i>  used  on  the  Jr'.kcs  o'  tic-  X.  a 
and  P  Sampling  scope  measurement'  wc.e  tin  t he 
device  shown  in  Figs  II  and  I'  and  other  device'  no! 
reported  here  At  first  the  spectrum  an.itv/er  u.o  set  up 
to  detect  frequencies  to  ISfiHz  In  tho  set  up.  frequeu 
cies  of  I  2.  I  b.  2  '  and  4  2  GHz  wete  detected  These 
oscillations  were  bias  dependent  but  ce.i'ed  at  drain  bias 
levels  in  excess  of  h\  The  oscillations  i cased  .t  sub 
st.inti.illv  higher  gate  bias  levels  f  nd  e  assumption 
that  there  were  additional  frequencies  m  excess  of 
18GHz  fha!  were  going  undetected,  and  that  he  inalogv 
to  two  terminal  devices  where  the  introduction  of  ,t 
inductor  reduces  the  circuit  frequence,  we  introduced  a 
•  nH  inductor  in  scries  with  the  source  loop  The  low 
frequency  oscillations  disappeared  and  in  their  place 
were  instabilities  at  IV  P  h  and  '.4  'GHz  More  detailed 
spectral  analysis  measurement'  were  suhscqucnilv  made 
with  the  spectrum  analvzer  open  to  frequencies  to 
4f)GHz.  The  pulsed  measurement'  were  then  analvzcd. 
with  the  following  details.  Oscillations  began  .it  A  x  \  iruj 
persisted  to  just  below  4.2'  \  Oscillations  were  obser¬ 
ved  at  f>  5.  IV  2V  20.  2(s  and  40GHz  Higher  frequencies 
could  not  be  detected  The  os. illations  were  strongest  at 
4  1  V  where  the  ‘.,5  and  4f)(iHz  signals  dominated  Al 
4.2'  V  on  the  drain  the  oscillations  ceased,  at  least  up  to 
40 (i Hz  In  the  same  measuring  circuit  and  with  the  gate 


.si  i  ;>  \  ■  i 


168 


HI  Grubin  m  at 


TEMPERATURE  (CK) 

Fig.  10.  Solid  curve  is  an  estimate  of  high  field  saturated  velocity  vs  ambient  temperature  (al  Velocity,  from  eqn  8.  ss 
ambient  temperature,  prior  to  oscillations,  (hi  Velocity,  from  eqn  8.  vs  ambient  temperature,  at  cessation  of  oscillations 

with  application  of  a  high  cate  bias 


wide  open  visible  light  was  detected  at  6.4  V  on  the 
drain.  This  was  not  accompanied  hy  any  oscillations 
whose  frequency  was  below  40GHz.  An  initial  reduction 
in  potential  on  the  gate  contact  resulted  in  similar  behavior 
but  at  -0.8  V  on  the  gate, all  oscillatory  activity  ceased, 
including  the  appearance  of  light.  The  electrical  behavior 
of  this  device  was  measured  at  a  variety  of  temperatures, 
with  qualitatively  similar  behavior. 

The  oscillation  frequencies  reported  above  bear  no 
clear  relation  to  the  transit-time  oscillation  frequency.  As 
indicated  in  Section  3.  the  gate  drain  separation  may  be 
too  short  to  guarantee  stable  domain  propagation  and  so 
the  oscillatory  properties,  including  the  frequency  are 
likely  to  be  determined  by  a  complex  device-circuit 
interaction.  There  is  also  the  possibility  that  the  oscil¬ 
lations  we  are  observing  are  a  consequence  of  feedback 


between  the  gate  and  drain  loop.  While  this  situation 
cannot  be  ignored  the  evidence  in  terms  of  the  six 
instability  criteria  given  in  the  last  section,  including  the 
suppression  of  oscillations  at  high  drain  bias  levels,  make 
this  possibility  small. 

In  the  above  paragraphs  we  reported  the  observation 
of  light.  Light  was  also  observed  under  d.c.  conditions 
The  details  are  as  follows.  Under  d.c.  conditions  light 
was  seen  at  the  beginning  of  the  current  instability, 
became  dimmer  as  the  drain  bias  increased  but  remained 
while  the  device  oscillated.  The  light  was  stronger  at 
drain  bias  levels  sufficient  to  eliminate  the  oscillations 
At  moderate  gate  bias  levels  faint  light  remained;  hut  at 
sufficiently  high  gate  bias  levels  large  enough  to  eliminate 
the  oscillations,  all  light  activity  ceased.  The  light,  for  all 
piactical  purposes  was  while  The  observations  of  light  is 
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PULSED 


Vn  (VOLTS) 

Fig.  1 1 .  (a)  d.c.  and.  (b)  pulsed  sampling  scope  measurements  of  the  current  voltage  characteristics  for  a  short  gate  FET 
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Fig  13.  Pulsed  current  vs  voltage  data  for  a  device  that  did  not  exhibit  any  oscillations 


compelling  evidence  for  the  presence  of  regions  of  high 
electric  field  and  excess  carriers.  At  low  drain  bias  levels, 
the  dimness  of  the  light  may  be  accounted  for  by  the  fact 
that  the  time  averaged  electric  field  at  a  given  point 
within  the  sample  is  less  for  propagating  domains  then 
for  the  trapped  domains.  The  latter  occurs  at  high  drain 
bias  levels.  The  absence  of  light  at  high  gate  bias  levels 
may  possibly  be  due  to  the  broad  region  of  charge 
depletion  in  the  vicinity  of  the  gate  boundary.  The 
arguments  for  white  light  are  somewhat  uncertain.  We 
recall,  however,  that  the  separation  between  the  central 
portion  of  the  conduction  band  and  the  valence  band  is 
within  the  infrared,  while  the  separation  between  the 
satellite  conduction  and  valence  bands  is  within  the 
ultraviolet.  The  entire  visible  spectra  is  within  these  two 
bounds. 

We  note  that  Mimura  et  al  [23]  on  IO'’/cm'  GaAs 
FETs  also  observed  visible  white  light  at  high  drain  bias 
levels  For  devices  with  highly  doped  n  *  drain  contacts 
light  was  observed  at  the  n‘n  interface;  for  devices 
without  an  n‘  drain  contact  radiation  was  observed  at 
the  drain  metalization  edge  Mimura  et  al. [23],  attribute 
the  radiation  to  impact  ionization  within  a  domain  trap¬ 
ped  at  the  drain  contact.  The  intensity  of  this  light  was 
ooserved  to  increase  with  drain  bias  and  decrease  as  the 
gate  voltage  was  made  negative.  These  results  are  con¬ 
sistent  with  our  findings. 

S.  SOME  MODELLING  CONSIDERATIONS 

The  results  of  our  analysis  are  consistent  with  the 
identity 

IDSS  y  R„  =  E„J.  (15) 


where  is  the  longitudinal  electric  field  in  the  source 
to  gate  region.  In  devices  with  K  greater  than  unity. 
will  generally  be  at  least  equal  to  F.,  \  and  so  for  this  case 
a  lower  bound  for  IDSS  exists 

IDSS  x  R„  S’  F.J..  (16) 

For  example,  with  the  device  reported  in  Figs.  8  and  9. 

the  room  temperature  values  of  the  l.h.s.  of  eqn  (16) 

(from  Table  1)  is  equal  to  1.9  V  The  r  h.s.  side  with 
E,  =  2.2  kV/cm  is  equal  to  1 .87  V  For  the  device  whose 
properties  are  displayed  in  Figs  1 1  and  12.  IDSS  x  R0  = 
0.52  V.  while  E,L-  0.77  V.  By  way  of  comparison 
Shockley's  study  [24]  teaches  that 

IDSS  x  R„  =  V,J1  (17) 

while  more  recently  Shur|25]  have  concluded  that 

IDSS  x  R„  =  (18) 

Equation  (18)  yields  results  significantly  lower  than 
observations,  and  less  than  eqn  (16). 

*.  CONCLUSIONS 

In  conclusion,  through  analysis  and  experiment  we 
have  presented  evidence  for  the  presence  of  negative 
differential  mobility  induced  current  instabilities  in  GaAs 
FETs.  We  have  listed  six  criteria  that  identify  these 
NDM  instabilities,  paramount  among  these  is  the  fact 
that  the  longitudinal  electric  field  under  the  gate  must  be 
greater  than  the  threshold  field  for  negative  differential 
mobility  while  the  current  density  must  exceed  the  sus- 
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taining  current.  The  fact  that  the  instabilities  in  some 
cases  are  delayed  until  higher  values  of  bias  are  reached 
implies  that  the  field  must  exceed  the  NDM  threshold 
field  for  the  instability.  We  have  also  presented  cor¬ 
roborating  experimental  and  theoretical  evidence,  con¬ 
sistent  with  the  conclusions  of  Yamaguchi  el  a/.  [8],  that 
devices  whose  current  densities,  within  the  source  to  gate 
region,  are  less  than  J,  will  not  go  unstable.  Those  whose 
current  densities  exceed  J,  will  sustain  instabilities. 
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ABSTRACT 

We  compute  the  transient  response  of 
gallium  arsenide  submicron  two  terminal  devices 
and  examine  their  load  line  dependence. 
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and  is  referred  to  as  a  'displaced'  or  'shifted' 
Maxwellian  (see,  e.g.,  Butcher,  reference  3). 


For  the  small  device  (SD)  whose  active 
region  length  is  typically  between  a  mean  free 
path  and  two  microns,  the  average  electron 
drift  velocity,  V,  cannot  be  represented  as 
an  explicit  function  of  electric  field.  Rather, 
since  the  important  time  processes  such  as 
transit-time,  momentum  and  energy  relaxation 
time,  are  comparable,  the  carrier  velocity  must 
be  calculated  as  a  function  of  space  and  time 
(see,  e.g.,  Ruch,  reference  I).  In  a  typical 
situation  where  the  SD  is  connected  to  a  power 
supply  through  a  load,  the  velocity  response 
will  depend  on  the  properties  of  the  load.  We 
demonstrate  this  dependence  for  a  purely 
resistive  load. 

The  study  is  numerical  and  involves  solving 
the  blocks  of  equations  in  figure  1.  Blocks 
I  and  IT  are  apparent.  Block  III  contains  the 
device  physics.  For  example  in  the  case  of  long 
GaAs  devices  the  carrier  dynamics  have  been 
successfully  represented  by  explicit  relations 
between  field  and  velocity,  and  field  and  dif¬ 
fusion  (see,  e.g.,  reference  2).  For  the  SD 
a  more  accurate  description  of  the  carrier 
dynamics  is  required  and  we  solve  the  Boltzmann 
transport  equation  (3TE). 


For  multivalley  semiconductors,  such  as  GaAs, 
the  subscript  'i'  =  1 ,2 , . . . identif ies  the  specific 
valley.  And  for  each  valley,  the  carrier  density, 
n^,  momenta,  p^,  and  temperature  T^,  must  be 
determined.  Each  of  these  quantities  have  a 
space  and  time  dependence,  and  vary  as  the  field, 
F,  varies.  In  the  calculations  below,  we  ignore 
spatial  dependence  and  solve  the  homogeneous  zero, 
first  and  second  moments  of  the  BTE  (see,  e.g., 
references  4  and  5).  These  are,  respectively 


(particle  conservation) 

V"i  -  ”ni  +  "j 

3 1  t  r 

ni  "j 

(momentum  balance) 
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nipi 
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(energy  balance) 

+  2"ikoTl\ 

at  \  2  2  / 


(2) 
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II 

The  solutions  to  the  BTE  are  approximate.  We 
assume  that  the  electron  distribution  is 
Maxwellian,  and  peaks  at  a  finite  value  of 
momenta.  The  distribution  function  has  the 
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IV 


where  a j  represents  the  number  of  equivalent 
valleys,  m- ,  the  effective  mass  of  the  ith  valley 
ami  *  l/v*  represents  density,  momenta,  and  energy 
scattering  rates.  Table  1  identifies  these  rates. 
Table  2  lists  the  values  of  the  relevant  para¬ 
meters  used  in  this  calculation.  The  conse¬ 
quences  of  varying  these  parameters  will  be 
discussed  at  a  later  time. 


Ill 

The  calculation  is  performed  for  the  situa¬ 
tion  where  the  SI)  is  part  of  a  resistive  circuit. 
Hn*  SI)  is  uniformly  doped  to  Nq  part  ic les/cc . 
Initially  all  carriers  have  zero  velocity.  A 
bias  0  across  the  load  and  SI)  is  then  stepped 
from  zero  to  a  finite  value.  The  time  depend¬ 
ence  of  V  and  the  potential  across  the  SI)  art* 
constrained  bv  the  load  line  equation 

0O  =  0c;n  +  RNQeV(t)A  (h) 

where  A  is  the  cross  sectional  area  of  the  SI) 
and 

N0V(t)  =  n1p1/m1  +  a2n2p2/m2  ,  (7) 

is  the  net  drift  velocity.  The  subscripts  '1' 
and  '2'  denote  respectively,  the  central  and 
satellite  valleys.  The  solutions  to  the  above 
equations  along  with  those  represented  by  blocks 
t  and  IT  of  figure  1  result  in  a  uniform  time 
dependent  change  in  field  and  velocity  across 
ttie  SI).  Because  of  the  assumption  of  uniformity, 
carriers  at  the  anode,  as  well  as  at  other  points 
in  the  SD,  have  the  same  velocity  as  carriers  at 
the  cathode.  Apart  from  the  artificial  nature 
of  this  assumption,  the  calculated  velocity  can 
be  used  to  determine  the  time  T,  for  a  carrier 
with  a  zero  initial  drift  velocity  to  travel  a 
distance  L,  in  a  time  varying  field: 

T 

L  -  f  V( t)dt  (8) 

The  unifotm  field  calculations  are  relevent  to 
transit-time  considerations. 

Examples  of  the  calculations  are  shown  in 
figures  2  through  4  for  different  values  of  the 
ratio  R/Rq.  Here  Rg  -  L/(NgevigA),  where  pg  is 
the  low  field  mobility  of  the  central  valley 
corresponding  only  to  LO  phonon  intravalley 
scattering.  Also,  recalling  for  GaAs  with  doping 
levels  not  exceeding  10^/cc,  that  the  equili¬ 
brium  velocity  electric  field  relation  differs 
only  marginally  from  that  for  zero  impurity 
scattering,  we  perform  the  calculation  for 
zero  impurity  scattering. 


Figure  2  shows  V(t)  for  two  values  of  K/Kg. 
Hie  value  of  0  exceeds  the  threshold  field  for 
negative  differential  mobility  (NDM)  and  we  see 
the  familiar  overshoot  effects.  With  respect  to 
the  influence  of  the  load  we  see  that  the  peak 
velocity  is  higher  and  the  relaxation  time 
shorter  for  the  R/Kfj  -  0.01  circuit.  The  major 
contribution  to  this  difference  is  the  larger 
potential  drop  across  the  SI)  for  the  lightly 
loaded  circuit.  We  point  out  that  it  is  possible 
to  find  a  particular  pair  of  values  for  0  and 
R/R{)  for  which  the  response  times  for  light  and 
heavily  load  circuits  are  similar.  From  another 
point  of  view  the  load  may  play  a  more  signifi¬ 
cant  role.  Til  is  is  shown  in  figure  3  and  4  where 
we  display  the  velocity  versus  electric  field 
relation  as  a  function  of  SD  length. 

The  calculations  show  the  length  and  load 
dependence  of  the  NDM  slope  for  the  two  terminal 
SD.  In  our  calculations  NDM  arises  only  from 
electron  transfer.  The  shortest  device  exhibits 
the  highest  velocities  and  the  weakest  NDM  region 
(off  scale),  both  features  reflecting  the  fact  that 
the  transit  times  are  considerably  shorter  than 
the  thermal izat ion  times  (see  also  reference  1). 
With  regard  to  the  load  line  dependence  we  see, 
for  a  SD  of  given  length,  the  presence  of  a 
’softer'  NDM  region  and  somewhat  higher  velocities 
when  situated  in  the  more  heavily  Loaded  circuit. 
The  consequences  of  this  result  are  somewhat 
difficult  to  predict  at  this  time,  insotar  as  the 
onlv  detailed  device  analysis  concerned  with  the 
influence  of  NDM  slopes  and  velocities  comes 
from  studies  on  longer  devices.  In  these  earlier 
studies  we  find  that  shallow  NDM  regions  lead  to 
broad  domains,  a  feature  important  to  device 
stability  as  well  as  gain  when  operated  as  a  two 
terminal  amplifier.  On  the  other  hand,  the  NDM 
slope  is  less  significant  than  the  ratio  of  the 
'peak'  and  'valley*  drift  velocities  when  the 
device  is  operated  as  a  uniform  field  oscillator 
(see  reference  2).  For  two  terminal  SD  the 
consequences  of  the  load  line  dependence  must 
await  further  studies  which  include  both  reactive 
circuit  elements  and  nonuniform  fields. 

In  conclusion,  through  solutions  of  the 
Boltzmann  transport  equation,  the  circuit  and 
Poissons  equations,  we  have  demonstrated  that 
the  transient  characteristics,  including  the 
drift  velocity  versus  electric  field  relation  are 
dependent  on  the  properties  of  the  external 
circuit . 
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TABLE  L 

SCATTERING  PROCESSES  : 

(CV)  CENTRAL 

VALLEY 

(SV)  SATELLITE 

VALLEY 

1/t 
,  nl 

l'-n2 

1  /  '  K1 

1  /  x  1 2 

CV  to  all  SV 

All  SV  to  CV 

\h  s2 

Intravalley  CV  and  CV  to 
Intravalley  SV  and  SV  to 
SV  and  one  SV  to  CV 

all  SV 
remaining 

CV  to  one  SV 

TABLE  2 

PARAMETERS  USED  IN  CALCULATION 


PARAMETERS 

r(ooo) 

X(100) 

COMMON 

Number  of  equivalent 

l 

3 

valleys 

Effective  mass  (me) 

0.06  7 

0.40 

T-X  Separation  (ev) 

0.36 

Lattice  Const,  (o^) 

3.64 

Density  (gm/cnw) 

5.37 

POLAR  OPTICAL 

SCATTERING 

Static  Diel.  Const. 

12.53 

High  Frcq.  Diel.  C 

10.82 

L0  phonon  (ev) 

0.0354 

f-X  Scattering 

Coupl .  Cons .  (ev/cm) 

0.62U109 

Phonon  energy  (ev) 

0.0300 

X-X  Xcatterlng 

Coupl.  Const,  (ev/cm) 

1.064xl39 

Phonon  energy  (ev) 

0.0300 

Acoustic  Scattering 

Deform.  Pot  (ev) 

7.0 

7.0 

Acoustic  Vel.(c/s) 

5.22x10 

Figure  1.  Summary  of  device  simulation. 


TIME  (T/Tq) 


Figure  2.  Magnitude  of  drift  velocity  versus 
time  for  different  resistive  loads.  Here  the 
bias  field  is  25  kv/cm  Vq  *  3.68  x  10^  cra/sec 
and  Tq  *  .315  psec. 
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GAAS  (1): 
GAAS  (2): 
GAAS  (3) : 

GAAS  (3)  1 
GAAS  (4) : 

GAASFET: 


IV.  NUMERICAL  CODES,  PARTIALLY  SPONSORED  BY 
ONR,  CONTRACT  N00014-78-C-0269 


BOLTZMANN  TRANSPORT  EQUATION  (BTE)  SCATTERING  INTEGRALS. 
UNIFORM  FIELD  BTE  TRANSIENT  EQUATIONS. 

ONE-DIMENSIONAL  NONUNIFORM  FIELD  BTE  SPATIAL  AND  TEMPORAL 
TRANSIENTS. 

BTE  INTEGRALS  FOR  THREE-LEVEL  TRANSFER. 

TWO-DIMENSIONAL  NONUNIFORM  FIELD  BTE  SPATIAL  AND  TEMPORAL 
TRANSIENTS. 

FIELD  EFFECT  TRANSISTOR  SIMULATIONS 
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VI.  FIGURE  CAPTIONS 


1.  Flow  diagram  representing  motivating  disciplines  in  device  physics. 

2.  Scattering  rates  for  the  parameters  of  Table  1.  Integrals  are  from  Ref.  18. 

CV  refers  to  the  central  valley  and  SV  to  the  staellite  valley. 

3.  Field  profiles  for  a  high  field  propagating  domain.  Dotted  curve  is  obtained 
from  the  mobility  equations.  Solid  curve  is  obtained  from  moment  equations. 
From  Ref.  12. 

4.  Schematic  representation  of  a  coplanar  two-dimensional  device  with  a  cathode 
notch. 

5.  For  No=10^',/am^  and  xQ=.3ps,  transient  central  valley  behavior.  (a)  Momentum 
versus  time,  (b)  Temperature  versus  time,  (c)  momentum  density  versus  time. 

b.  As  in  Fig.  5,  but  for  satellite  valley. 

7.  Mean  velocity,  from  Figs.  3  and  6.  Here  »2mck0T0  =  3.7  x  10^  cm/sec. 

8.  Time  dependent  bias  field  used  to  show  rise  time  dependence  of  transient 

velocity.  For  all  of  these  calculations  FqM  =  21.3  kV/cm  and  F^  =  .43  kV/cm. 

9.  For  this  pulse  t^  =  500/xo  t2  =  600/xQ,  t^  =  1200/x  .  x0  =  .3ps.  (a)  Mean 

velocity  versus  field  (Fq  =  4.3  kV/cm),  (b)  central  valley  population  versus 
time,  (c)  central  valley  temperature  versus  time. 

10.  As  in  Fig.  9.  Here  t^  =  50/x^,  t2  =  60/xo,  t3  ^  110/ T' 

11.  As  in  Fig.  9.  Here  t^  =  5/xq,  t2  =  10/tq,  t3  =  15 hQ. 

15  3 

12.  Self-excited  oscillation,  for  the  circuit  displayed  in  Fig.  13.  Nq  =  10  /cm  . 
Oscillation  frequency  is  78  GHz. 

13.  Maximum  frequency  for  self-excited  oscillations. 

14.  Length  dependence  of  negative  differentail  mobility. 

15.  Schematic  representation  of  device  and  circuit  configuration  for  submicron 
homogeneous  field  profiles.  Inhomogeneous  doping  profile  is  treated  as  a 
"jellium"  distribution. 

16.  (a)  Fractional  central  valley  population  versus  distance,  (b)  Central  valley 
temperature  versus  distance,  (c)  Staellite  valley  temperature  versus  distance. 
Computation  occurs  at  time  t fx  =4.  x0  is  the  LO  phonon  internally  scattering 
time  and  equals  0.32  ps.  T  in  this  figure  is  room  temperature.  The  gradient 
of  N  is  zero  at  the  boundaries.  All  nonuniformities  are  due  to  the  notch. 
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VI.  FIGURE  CAPTIONS  (Cont'd) 


(a)  Mean  velocity  versus  distance,  (b)  Field  versus  distance,  (c)  Free 
carrier  density  versus  distance  at  time  t/i  =  4. 

As  in  III  13,  but  at  t/x  =  16. 

As  in  III  14,  but  at  t/r  =  16. 

Boundary  Condition  approximations  to  one-  and  two-dimensional  problems. 


TABLE  1 


PARAMETERS  USED  IN  CALCULATION 


Parameters 

T(000) 

X(100) 

Common 

Number  of  Equivalent  valleys 
Effective  mass  (me) 

T-X  Separation  (ev) 

Lattice  Const.  (oA) 

Density  (gm/cm^) 

1 

0.067 

3 

0.40 

0.36 

5.64 

5.37 

POLAR  OPTICAL  SCATTERING 

Static  Diel.  Const. 

High  Freq.  Diel.  C. 

LO  phonon  ( ev) 

12.53 

10.82 

0.0354 

T-X  Scattering 

Coupl.  Cons,  (ev/cm) 

Phonon  energy  (ev) 

0. 621xl09 
0.0300 

X-X  Scattering 

Coupl.  Const,  (ev/cm) 

Phonon  energy  (ev) 

1.064xl09 

0.0300 

Acoustic  Scattering 


Deform.  Pot  (ev) 
Acoustic  Vel.  (c/s) 


7.0  7.0 


5.22xl05 
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